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PREFACE 


A sound mathematical training for the modern engineer 
is a sine qua non for attaining new heights in all aspects 
of engineering practice. One of the areas of mathematics 
that plays a big role in the mathematical education of 
the engineer is vector analysis, which is now invariably 
included in the curriculum of higher mathematics in 
engineering colleges. 

The present collection of problems in vector analysis 
contains the required minimum of problems and exercises 
for the course of vector analysis of engineering colleges. 

Each section starts with a brief review of theory and 
detailed solutions of a sufficient number of typical prob- 
lems. The text contains 100 worked problems and there 
are 314 problems left to the student. There are also a 
certain number of problems of an applied nature that 
have been chosen so that their analysis does not require 
supplementary information in specialized fields. The 
material of the sixth chapter is devoted to curvilinear 
coordinates and the hasic operations of vector analysis 
in curvilinear coordinates. Its purpose is to give the reader 
at least a few problems to develop the necessary skills. 

The exposition in this text follows closely the lines 
currently employed at the chair of higher mathematics 
of the Moscow Power Institute. 

The present text may be regarded as a short course in 
vector analysis in which the basic facts are given without 
proof but with illustrative examples of a practical nature. 
Hence this problem book may be used in a recapitulation 
of the essentials of vector analysis or as a text for readers 
who wish merely to master the techniques of vector 
analysis, while dispensing with the proofs of propositions 
and theorems. 

In compiling this problem book, the authors made 
extensive use of material in published courses of vector 
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calculus and collections of problems. Many problems 
were made up by the authors themselves. 

This collection of problems is designed for students of 
day and evening departments at engineering colleges and 
also for correspondence students with a background of 
vector algebra and calculus as given in the first two years 
of college study. 

We would like to express our sincere gratitude to 
Professor V. P. Gromov (the Moscow Krupskaya Pedagog- 
ical Institute), Professor A. V. Efimov and Associate 
Professors I. M. Petrov, B. I. Fridlender, and V. N. Zem- 
skov (the Institute of Electronics) for their thorough 
scrutiny of the manuscript of the book and for valuable 
tuggestions and remarks that were made full use of in 
she final editing. 


M. L. Krasnov 


A. I. Kiselev 
G. I. Makarenko 


Moscow-Dubna, 1977. 


CHAPTER І 


THE VECTOR FUNCTION 
OF A SCALAR ARGUMENT 


Sec. 1. The hodograph of a vector 
function 


Definition 1. A vector r is said to be the vector function 
of a scalar argument t if each value of the scalar taken 
from the domain of admissible values is associated with 
a definite value of the vector r. This can be written as 
follows: 

r =r (t). 


If the vector r is a function of the scalar argument t, 
г =r (t), 


then the coordinates z, y, z of the vector r are also func- 
tions of t: 


z=z(t), y-y(t 2=2 (0. 


Conversely, if the coordinates of the vector r are func- 
tions of t, then the vector r itself is also a function of t: 


r — z(t)i + y(t) j +2 (0) К. 


Thus, specifying a vector function r (t) is the same as 
specifying three scalar functions z (t), y (t), z (t). 

Definition 2. The hodograph of the vector function 
r (t) of a scalar argument is the locus described by the 
terminus of the vector r (t), as the scalar t varies, when 
(i origin of the vector г (t) is fixed at a point O in space 
Fig. 1). 

The hodograph of a radius vector г = г (t) of a moving 
point is the trajectory L of that point. Some other line L, 


10 The Vector Function of a Scalar Argument (Ch. I 


(Fig. 2) is the hodograph of the velocity у = v (t) of that 
point. Thus, if a material point (particle) is in motion 
around a circle with constant ve- 
4 locity, | у | = constant, then its 
hodograph of velocities is like- 
S wise a circle with centre at O, 

4 and with radius equal to | v |. 
Example 1. Construct the ho- 


0 8 dograph of the vector r = ti+ 
tj + Ck. 
Fig. 1 Solution. 1°. This construc- 


tion may be carried out by using 
points and setting up a table: 


ci i ol 4 


t 0 
o| i+j+k 31 + 3j + 9k | 4i + 4j + 16k 


r 2i-+ 2j+-4k 


2°. Alternative solution. Denote by z, y, 2 the coordi- 
nates of vector r; we have 


T4, узі, 2 == Й, 


Eliminating the parameter ¢ from these equations, we 
get equations of the surfaces у = =, 2 = 2’, the line L 


Fig. 2 


of intersection of which is what defines the hodograph of 
the vector г (t) (Fig. 3). 
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1. Construct the hodographs of the following vectors: 
(a) r= 21+ t?j — ?k. 


о. 
(b) r= vy it wey 


(с) r=cost-i+sint-j+k. 


Fig. 3 


(d) к= ti - 3 tj E tk. 


.. 2t 4-2tj + ((1—2) К 
(e) r= 2 о . 


Sec. 2. The limit and continuity 
of a vector function of a scalar argument 


Suppose a vector function г = г (t) of a scalar argument 
t is defined in some neighbourhood of the value ¢, of the 
argument t, except perhaps for the value t, itself. 

Definition 1. A constant vector A is said to be the 
limit of the vector г (t), as t — to, if for any e > 0 there 
is аб >0 such that for all t 5& t, that satisfy the con- 
dition | — ty | < ô the following inequality holds true: 


Ir) = А |<. 
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As in the case of ordinary calculus, we write lim r(t) = 


t-^to 

= А. 

Geometrically (Fig. 4), this means that the vector 
г (t) tends, as і — £to, to the vector A both in length and 
in direction. 

Definition 2. A vector 
æ (t) is said to be infini- 
tesimal, as t— to, if 
a(t) has a limit, as t — to, 
and that limit is zero: 


lim @ (t) = 0, 
tto 
or, what is the same 
thing, if for any е 0 
Fig. 4 there exists a 6 >0 such 
that for all t = t, that 
satisfy the condition |t—t, |< ô, the inequality 
1 a(t) | < e holds true. 
Example 1. Show that the vector a (t) = ti + sin tj is 
infinitesimal when t > 0. 
Solution. We have 


ja (2) | = |ti -- sin tj] < |t| + |sin t| <2 |t]. 


From this it is evident that if for every e — 0 we take 
$ = e/2, then for | t — 0 | < 6 = e/2 we have | а (t) | < 
< e. By the definition, this means that æ (t) is ап infini- 
tesimal vector when #— 0. 

2. Show that the limit of the modulus of a vector is 
equal to the modulus of its limit (if the limit exists). 
3. Demonstrate that for a vector function г (t) to have 
a limit A, as tt, it is necessary and sufficient that 
г (t) be representable in the form 


г (t) = А {а (t), 
where c (t) is an infinitesimal vector when t — fp. 


4. Show that if the vector functions a (t) and b (t) have 
limits as £ — t,, 


lima(t)=A,  limb(t) = B, 
tto 


tto 
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then their sum a (t) + b (t) and their difference a (t) — 
— b (t) also have limits as t — і, and 

lim [а (t) +b (t)) = А + B. 

tto 


5. Let 
lima(t) 2A, limb(t) = B. 
t-to 


t+to 
Prove that 
lim (a (t), b (t)) = (А, B), 
t-to 
where (a (t), b (t)) із a scalar product of the vector func- 


tions a (t) and b (t). 
6. Let 


r(t)—z()i-c-y(0i)-c-z()k, A= ai + aj + 
+ a,k. 
Show that if lim r(t)=A, then 
+70 


limz(t) 2a, limy(t)=a,, limz(t)=as. 
1-40 toto 1-to 


Find the following limits: 
7. lim (= 14 282—1 jre k). 


: iri t 


. sint 
9. lim (= i+cost- dV). 


10. lim ( sine ap ERIT yt), 


11. ЕЕ те 3-2). 


Definition 3. A vector function r = г (t) defined in 
some neighbourhood of the value t = f, is said to be 
continuous when t — t, if 


lim r (t) =r (to). 
tte 
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In other words, г = г (t) is continuous for ¢ = t, if for 
every e >O there is a 6 >O such that for all t that 
satisfy the condition |1 — 4, | < $ the inequality 
|r (t) — г (to) | < е holds true. 

The hodograph of a continuous vector function of a 
scalar argument is a continuous curve. 
12. Start with the familiar inequality |а — b | > 
> ||а| — |Ь | | and demonstrate that the continuity 
of a vector function implies the continuity of its modulus. 
Is the converse true? 
13. Show that if a (t) and b (t) are continuous for t = В, 
then the vector function a (t) + b (t) is also continuous 
for t = ty. 
14. A vector function а (t) + b (t) is continuous for 
t = 1,. Does it follow from this that the vectors a (t) 
and b (t) are also continuous when ? = tọ? 
15. Prove that if a(t) and b (f) are continuous vector 
functions, then their scalar product (а (t), b (t) and 
vector product [a (2), b (t)] are also continuous. 


Sec. 3. The derivative of a vector 
function with respect to 
a scalar argument 


Suppose a vector function г = г (t) is defined for all t 
on the interval (to, ї,). Take some value ż € (to, &), then 
give t an increment А? such thet t + At € (to, В) and find 
the corresponding increment Лг = г (t + At) — r (t) in 
the vector function г (t). Now consider the ratio Ar/At. 

Definition. If, as At — 0, the ratio Ar/At has a limit, 
then that limit is called the derivative of the vector 
function r = r (t) with respect to the scalar argument t 
for a given value ¢ of the argument and is denoted as 


dr (t)/dt or r' (t) or r (t). Thus, 


zo =lim E um lim 
t м-0 А-0 


r(t+At)—r(d) 
At Ы 


In this case the vector function г = г (t) is said to be 
differentiable. 


Sec. 3} The derivative of a vector function 15 


16. Show that if the vector function r — r (t) has а 
derivative for some value ¢ of the argument, then it is 
continuous for that value t. 

The derivative of a vector function г (1) of a scalar 
argument ¢ is a vector directed along the tangent to the 


ar 
Mat 


Fig. 5 


hodograph of the original vector at the point under con- 
sideration (Fig. 5). The vector dr/dt is in the direction 
of the terminus of the vector r (t) as it moves along the 
hodograph when the parameter ¢ increases. 

Suppose r — r (t) is the radius vector of a moving 
point. Then the vector v = dr/dt is the velocity vector 
of that point. 

Suppose 


r (t) = z (t)i + y (t) j + z (t)k, 


where the functions z (t), y (t), z (t) are differentiable at 
the point t. Then there exists dr/dt for that value of t and 


dr _ dz dy ., dz 
ae ar шт! К. (1) 


Example 1. Find dr/dt if r = ia cos t + jb sin t (the 
point is moving in an ellipse). 
Solution. From formula (1), 
dr 


== — ia sin t + jb cost. 
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By analogy with the differential of a scalar function, 
the differential of a vector function г = г (t) is a vector dr 
defined by the equality 


where df = At is the increment in the scalar argument t. 
As in the case of scalar functions, 


Ar = dr + o: At, 
where а = a (t, At) — 0 as At — 0. 


Basic rules for differentiating a vector function 


Assume that all functions being considered (both scalar 
and vector) are continuous and differentiable. 
1°. If c is a constant vector, then de/dt = 0. 
2°. The derivative of a sum of vector functions is equal 
to the sum of the derivatives jd the summands: 
d(a(t)+b(t)) _ 4 db 
dt = +r 
3°. Suppose a vector function a(t) is multiplied by 
a scalar function m (t) of the same scalar argument. Then 
(та) dà | dm 
d ae a 


(а 


5°. srl = [3 , ь]+[а, 4] 


(In this formula, the order of the factors a and b in the 
right-hand member must be the same as that in the left- 
hand member.) 

Let us prove formula 4°. We set q (t) = (a (t), b (t)). 
Give t an increment At; then, by the distributive prop- 
erty, we have for the scalar product 


Аф = 9 (t + At) — Ф (t) = (a + Aa, b + Ab) — (a, b) 
= (Aa, b) + (a, Ab) + (Aa, Ab), 


SpA, b)+ (ade) (A a). © 
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It is given that the functions a (t) and b (2) have deriva- 
tives for the value Е of the argument and, hence, are 
continuous for that value of t. Therefore 


A» да ë jim Sb. № and lim Ab—0. 


мо dt? Ano At dt м-0 

Passing to the limit in (2) as At > 0, we obtain 
d(a,b) da db 
dt oc a b)t (nx) 
17. Given r = г (t). Find the derivatives: 
d d d d d 

(а) 6. (кт). 00 ғ [535]: 
18. Prove that if the modulus | r | of the vector function 
r — r(t) remains constant for all values of t, then dr/dt | r. 
What is the geometrical meaning of this fact? 
„19: Prove that if e is a unit vector in the direction of 
the vector E, then 


[E, dE) 


[e, de] = ТЕГ 


20. Зиррозе 
u = и, (2, у, 2, t) i + из (2, у, 2, t) E + из (2, у, 2, t) К, 
where иу, из, из are continuously differentiable functions 
of their arguments, and z, y, z are continuously differentia- 
ble functions of t. Show that 

du _ да | Ju dz ou dy дий 
de atts dt + ay dt ^ às dt^ 
21. Find the trajectory of motion for which the radius 
vector г (t) of a moving point satisfies the condition 
dr/dt = [a, r], where a is a constant vector. 

The derivative dr/dt of the vector function r (t) of a 
scalar argument is a vector function of the same argument. 
If there exists a derivative of dr/dt, then it is called the 
second derivative and is denoted d*r/dt*. Generally, 

am d анылан 
dt^ dt \ ат у’ E EE 
22. Given the radius vector of a point moving in space: 
r{asint, —acost, bt?) 


(t is time, and a and b are constants). Find the hodographs 
of velocity and acceleration. 


2-910 
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23. Given: г = a cos ot + b зіп оѓ, where о, a, b are 
constants. Prove that 


(1) [r. | = tos, by, 
(2) r+ ot =0. 


24. Show that if г = аео! + Бего‘, where a and b are 
constant vectors, then d?r/dt — wr = 0, 

25. Show that the modulus of the differential of the 
radius vector of a point is equal to the differential of 
the length of the arc described by the point. 

26. Suppose a — a (u) is a vector function of a scalar u, 
where u in turn is a certain scalar function of the basic 
scalar t. Assuming a (u) and u — u (t) to be differentiable 
the necessary number of times, find an expression for the 
derivatives of the composite function da/dt, d*a/di?. 


Sec. 4. Integrating a vector function 
of a scalar argument 


Definition 1. We will say that the vector function А (t) 
is the primitive of the vector function a (t) when t, < t < 
< t, if A (t) is differentiable and 


E =a (t), tC(to,t,). 


Definition 2. The collection of all primitive functions 
of a(t) is termed the indefinite integral of the vector 
function of a scalar argument a = a (t). As in integral 
calculus, the indefinite integral of a vector function is 


denoted by the symbol |, and we have 
| a (t) dt =A (t)+C, 


where A (t) is one of the primitive functions of a (t), and 
C is an arbitrary constant vector. 

The following properties hold true for integrals of 
vector functions: 


1°. faa (t) dt =a аб dt (aisa numerical constant). 
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2°, | (a(t)+b(t)) dt = fae dt 4- | b(t) dt. 


27. Show that if е is a constant vector, and а (t) is a 
variable vector, then 


| (с, а (£)) dt = (e. | a (t) dt) А 


j le, а (2)] dt = [ e, fa а). 
If 
а (t) = a, (t) i + a, (t) j фа, (t) k, 
then 


face dt =i | a, (t) 4+} j a, (t) dt +k | а (1) 4. (1) 


That is, the integration of а vector function reduces to 
three ordinary integrations. 

Example 1. Find the indefinite integral for the vector 
function a (t) = i cos t + je™ + К. 

Solution. According to formula (1), 


[adii j cos tdt i | etat+k | dt = 
=isint—je*+kt+e, 


where с is an arbitrary constant vector. 
Find the integrals of the following vector functions: 
i k 
28. a (t) = te'i + sin? tj — TFF: 
ti РА 
Ей T te*j 4- cos t- k. 


30. а (t) = соз teint .§ —tcost?-j + К. 


29. a (t) = 


31. a (t) = £1 —t sin t-j 4 2'k. 


Let a vector function a (t) be defined and continuous 
over a certain interval [t,, 7], which is the range of the 
argument t. 

Definition 3. We define the definite integral of a vector 
function a (t) on the interval [t,, T] as the limit of the 


2* 
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vector integral sums 
n-1 


= 2 a (тһ) At,, Th Е UT e] 


as the length At of the largest of the subintervals [t,, t, ‘al 
(К =0, 1, ..., n — 1) into which the interval [t,, Т 
is partitioned tends to zero: 
T п-1 
faw dt=lim > a (Ta) Ath. 
At+0 
to k=0 
The following formula holds true: 
T 
} a (t) dt=A (T)—A (to), (2) 
to 
where A (t) is some primitive for the function a (t) on the 
interval [t,, T 


a (t) = a, ( i + a4 (t) j + as (t) k, 
then 
T 


T T T 
faar-i | a, (t) dt 4-j | a, (t) dt +k | аз(0) 4. (3) 
10 to 


to to 


л/2 
Example 2. Compute | а (t) dt, where а (t) = i созі — 
0 


—jsin?t. 
Solution. By virtue of formula (3), 
n/2 л/2 л/2 
| а(00=1 | cost at—i sin? t dt 
0 0 0 
. [3/2 ,/t — sin 2e \ [2/2 
=isine|" —] (-— 252) =1— 4: 


Compute the following integrals: 


л 
32. | 804, where а = sin?t cos ¢-i+cos?¢sint-j+k. 
0 
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33. | a (t) dt, where a— —5,— -— + ke’. 


34. Fi 


35. 


Ote y Ot ot—- 


а (t) dt, where а = Зл cos ла 2tk. 


а (2) 41, where а = (2 - л) 1-31.) + nk. 
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Example 3. An electric current J flows upwards along 
an infinite wire that coincides with the z-axis. Find the 


vector H of the magnetic Coen ear te 
field intensity set up by 4 
this current at an arbitrary 
point М (т, у, 2) of space 
(Fig. 6). 

Solution. We consider a 
sufficiently small element 
PP, = dt ofthez-axis. By 
the Biot-Savart law, the 
intensity dH of the mag- 
netic field set up at point 
M by the current flowing 
through element 45 of the 
wire coincides in direction 
with the vector product 
(dg, r,], where 4 = PP,, 
| dẹ |=%, т, = РМ (see 
Fig. 6). By this same law, the 
modulus of the vector dH is 


I YN 
IaH] = sin (46, г.) dt, 


Fig. 6 


рм(2,0,2) 


РАХ 
where (4$, г,) is the angle formed by the vectors d$ and r,. 


Since 
А /N 
1145, г1] | = 706 sin (46, г), 
we can write 
1 
dH = +14, F4). 


(4) 
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In order to obtain the desired vector H at the point M 
we have to sum all vectors dH pertaining to distinct 
elements PP, of the wire, that is, we have to integrate 
the expression (4) over the whole z-axis: 


+ < 


Н= | = (a, ri. (5) 


We have 
r, — OM — OP. 
But 
OM = zi + yj + 2k, OP = Lk, 
and therefore 
n = zi + yj + (2 — $) k 
so that 


ry = [ry] =V 22+ y? + (2— 0)? = Ир (2—0), 


where р = V T + у? is the distance of point M from the 
axis of the wire. 

For the vector product [45, rj] we have 
ijk 
00% 
zy 2—5 


196, г,| = = —iy db 4- jz dt, 


and formula (5) takes the form [the point М (z, y, 2) is 
fixed, J = constant] 


+o 


To compute the integral on the right-hand member of (6), 
make the substitution 
PA Shep 79%" ae 


& _ pdt 
C—z=ptant, d= =. 
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We then have 
+оо 1/2 


а m | pat 
| тесир = „ STOR OF таптр® 
— ео -л 
4 2/2 
= 7+ \ cos t dt = — 
-3/2 


Thus, the intensity vector H of the magnetic field is 
in our case given by 


I А А 
Ни 
or 
2 
H=- Il, г], 


where I = Г.К is the current vector, г is the radius vector 
of point M (z, y, 2) of the field, and p is the distance of 
M to the axis of the wire. 

Example 4. The motion of an electron in a homogeneous 
magnetic field. 

1°. Suppose a magnetic field H is set up in some region 
of space; let it be constant in magnitude and direction 
(a homogeneous field). Suppose at time t = t,, an elec- 
tron enters the field with an initial velocity vy. Determine 
the path the electron will take. 

Solution. First suppose the vector v, is perpendicular 
to H and that the initial position of the electron is at 
point Мо. Choose the origin О at an arbitrary point of the 
plane P passing through M, at right angles to the vector Н 
(Fig. 7). Let the initial radius vector OM, be ro, let r be 
the radius vector of the electron at the current instant of 
time ft, and let v be the instantaneous velocity at that 
instant. The basic differential equation of motion is 

dir 
It will be recalled that the force F acting at time t оп the 
electron by the magnetic field is 


Е = —e, [Н, vl, 
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where е„ is the absolute value of the electron charge. 
Thus, 


dtr 
m r= elv, Н]. (7) 


At every instant t, the force Е is perpendicular to the 
direction of velocity and to the direction of the field H; 


Fig. 7 


at every instant, it will force the electron to deviate from 
a rectilinear path and to describe a certain curvilinear 
trajectory. 

Let us rewrite (7) as 


dv dr 
m e |, Н | 


and integrate from t, to t with respect to t. This yields 


mv — mv, = eo lr, Н) — e, Ir, Н) 
or 
ту = elr, Н] + (nv, — ео Го, HI). (8) 


Now choose a coordinate origin O' such that the term in 
parentheses in the right-hand member of (8) vanishes, that 
is, so that 

ео го, Hl = ту. (9) 
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From (9) it follows that the initial vector го must be 
perpendicular to the vector у, and must lie on the straight 
line М,К, which is perpendicular to the plane of the 
vectors v, and Н. By virtue of (9), the modulus of the 
vector го must satisfy the relation 

eo | Го || H | = т | у |, 
whence 


_ m |vol 
[го] EM eo | HI М (10) 


This determines the position of the new origin O’. Relative 
to this origin, equation (8) is rewritten thus: 


mv = e, lr, H] (11) 
or 
m S. — e, [r, Н]. 42) 


From equation (11) it follows that the trajectory of the 
electron is a plane curve lying in the plane P because at 
every instant the vector v is perpendicular to H. Now take 
the scalar product of both sides of (12) by r: 


m (т. Gp) = 6209.17. 80). (13) 
The mixed product in the right-hand member of (13) 


is zero, so that 
dr 
(r. A) =0, 


$ (r2) =0 or 4r (72) = 0, that is, г? = constant. 
This is the equation of a circle lying in the plane P with 
centre at the chosen point O'. The radius of the circle is 
found from formula (10) since the initial point М, must 
also lie on that eircle. Thus, we finally have 
r= T= (10) 
Thus, if an electron enters a homogeneous magnetic 
field H with an initial velocity v, at right angles to H, 
then it will describe, in that field, a circular trajectory 
lying in the plane P perpendicular to H and passing 
through the initial point. The radius of the circle is 


whence 
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given by formula (10) and its centre O’ lies on the straight 
line perpendicular to the plane of the vectors v, and Н; 
note that a rotation from v, to H must be seen from point 
О' as a counterclockwise rotation. 

From (10) it is evident that the radius r, of the circle 
is inversely proportional to | H |. Thus, the greater the 
intensity of the magnetic field, the greater the curvature 
of the trajectory. 

From formula (11), 


ту = e, ir, H), 


it is clear that if r is constant in modulus and is all the 
itme perpendicular to H, then also the velocity v of the 


W 


Fig. 8 


point will be constant in magnitude, 
|v | = v = constant, 


so that the electron is in uniform motion in the orbit. 
The period of revolution T is 


T = 270 On (14) 
Vo 


cU. 
eo | H] 7 
This formula does not involve the initial velocity vo. 
Thus, irrespective of the initial velocity у, which is 
perpendicular to H with which the electron enters the 
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homogeneous magnetic field H, it will perform a single 
orbital revolution and always in the same time 7. 

2°. Now suppose an electron enters a homogeneous 
magnetic field H with some initial velocity V that is 
not perpendicular to the vector H. This velocity may 
then be resolved into two components: the vector v, 
at right angles to the field, and a vector v, parallel to the 
magnetic field. 

From the formula 


Е —e[V, Н] = e Ivo, Н) 


it is evident that the "twisting" force F is given only by 
the perpendicular component v, and that it imparts to 
the electron a rotational motion about the circle (centred 
at О’) discussed above. As for the other component v,, 
the electron will retain it by inertia and, besides having 
a uniform circular motion, it will have a rectilinear and 
uniform motion in the direction of H with a velocity 
v = |У | соѕ а. The combination of these motions 
yields a helical curve with axis parallel to the vector H 
and passing through the point O’ (Fig. 8). 


Sec. 5. The first and second derivatives 
of a vector with respect to the arc length 
of a curve. The curvature of a curve. 
The principal normal 


Consider a curve L in space. On it, choose a point M, 
as the origin and also choose a direction along L that will 
be regarded as positive. For a parameter, take the arc 
length s reckoned from M, of the curve (Fig. 9). Then 
the radius vector of a point M of the curve is 


r — r (s). 


With that choice of parameter, 


where t? is a unit vector directed along the tangent to 
the curve L in the direction of increasing values of the 
parameter s. 
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If the vector r is given by the coordinates 
r = zi + yj + zk, 
then 
аг. dy,.,d 
v-z it iti К 
and 
dz \2 dy \2 dz \2 — 
(Ж) +(2) (6) =. 
Since | 1° | = 1, the vector dx?/ds is orthogonal to the 


vector 70. 
The modulus of the vector dr°/ds is 


dv? 
| |- к. 
Неге, К is the curvature of the curve L at the point M. 


The straight line having the direction of the vector 
dv°/ds and passing through the point М of the curve is 


Fig. 9 


termed the principal normal of the curve at the point M. 
Denoting the unit vector of that direction by по, we have 
£t Rn (4) 


ds 
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The inverse of the curvature of a curve at a given point 
is called the radius of curvature of the curve at that point 
and is denoted by R: 


1 
В=-. 


Thus, formula (1) may be rewritten as 


dr dy Jm 
ds? ds R^ 
From this, 
1 d?r 
K= R |as 
or 


к=+=ү (9) (m) (Rm). © 


Using (2), we can compute the curvature of a curve at any 
point if the curve is specified by parametric equations in 
which the parameter is the arc length s. 

In the particular case of a plane curve, a circle of 
radius a, 


$ 

z=acos—, 

а 

P . $ 

уёаз =, 

we have 

dz 1 $ азу __ 1. s 
qi et at ae aua. 


and formula (2) yields 


ride ndo ic rec d 
К=-к= ar cos ata sin ата: 


This means that the curvature of a circle of radius а is 
constant and is equal to the inverse of the radius of the 
circle. 
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If the curve Z is given by the vector-parametric equa- 
tion r = г (t), where the parameter ¢ is arbitrary, then 


lar ae | 

4 dt ^ dt? 

R dr |3 2 (3) 
à 


Formula (3) permits computing the curvature of the 
curve at any point provided we have an arbitrary para- 


metric specification of that curve. 
Example 1. Compute the curvature of the helical curve 


r=acost-i+ asin t-j + ktk. 
Solution. Since 


K= 


Ca —asint-i+acost-j+hk, 
a 
ore —acost-i—asint-j, 
the vector product 
P i j k 
(+. ав |= —asint acosth 
—acost —asintO 
=sahsint-i—ahcost-j+ ak. 
Consequently, 
d d? —UTI d = Le 
[i э |а Иа, || Иа. 
By virtue of (3), 
1 a 
K=p = aq 
or 
3 
RO — constant. 


Thus, a helical curve has a constant radius of curvature. 
Find the radius of curvature of each of the given curves: 

36. r = In cos t-1 + In sin t-j + V 2t-k. 

97. r = ťi + 28}. 
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38. r = 3t*i + (3t — 13) j + 2k for t = 1. 

39. r = а (cost + t sin t) i + а (sin t — t cos t) j 

for t = n/2. 

40. г = a cosh t-i + a sinh t-j + atk at any point t. 


Sec. 6. Osculating plane. Binormal. 
Torsion. The Frenet formulas 


The plane passing through the tangent line and principal 
normal to a given curve L at a point M is termed the 
osculating plane at the point M. 

For a plane curve, the osculating plane coincides with 
the plane of the curve. 

If the vector г = г (t) has a continuous derivative 
dr/dt in the neighbourhood of а point & and, besides, a 
Second derivative d?r (t,)/dt? such that 


dr (to)  d'r(to) 
[ ae cu | #9, 


then at the point ¢ = t, there is an osculating plane ‘o 
the curve г = г (1) whose vector equation is 


(oro, [ т || 


where p — p (t) is the radius vector of the current point 
of the plane. 

The normal to the curve at the point M, which normal 
is perpendicular to the osculating plane of the curve at 
that point, is called the binormal of the curve at the given 
point M. 

Denote by b? the unit vector of the binormal oriented 
so that the vectors то, n°, b? form a right-handed trihedral 
(Fig. 10). Then 


Ьо = 1, b9— [10, по]. 
For the derivative db°/ds we get 


The vector db/ds is perpendicular both to the vector 1? 
and to the vector - b^, that is, it is collinear with the vec- 
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tor n°. Set 
del oi 
ds | T° 
We then have 
ab? _ 1 
d = 


The quantity 1/T is termed the torsion of the given 
curve, and Т is called the radius of torsion of the curve. 

The torsion of a curve 
is given by the formula 


1 dr d'r dr 
pu (agra): 
where the symbol (a, b, c) 
denotes a mixed product of 
the vectors a, b, c, that is, 
(a, b, c) = (a, (b, cl). 

For the case where the 
curve is given by the vec- 
tor-parametric equation r — 
—r(t) we have 


(4,5 a 
dt’ dt?’ às) 


Те, ==]. (0 


La 
D 


"dt ^ at? 


Example 1. Find the torsion of the helical curve 
г = a соз t-i + a sin t-j + МЕ. 


Solution. We find the derivatives of the given vector: 


“= —asint-i+acost-j+hk, 
d*r 5 i F 
ча = —acost-i—asint-}, 


=, =asint.i—acost-j. 


ат 
dt? 
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The mixed product of these vectors is 
—asin і acost h 
)= —acost —asint 0|=а?һ. 
asint —acost 0 


(2. da dr 
dt ’ data’ а 


In example 1, Sec. 5, we found that 
d dir 3|? 
[a> ar] = ee. 


Using (1), we obtain for the torsion 


LE 
Т ah 
Thus, the torsion of a helical curve is the same at all 


its points. 
Example 2. Write the equation of the osculating plane 
at the point £ = 0 of the helical curve 
r = a cos t-i + asin t-j + ktk. 


Solution. We find the values of the derivatives of the 
given vector and its derivatives dr/dt and d?r/di? at the 


point t = 0: 
3 
r(0)=ai, "f capra, SD 
Consequently (see example 1, Sec. 5), 


dr (0 dêr (0 : 
(SP. Set |= -aj ak. 


The vector equation of the osculating plane is 


dr (0 d?r (0 
(p—r(0), au , 20.) =0 


= — аі. 


ог 
(p — аі, —ahj + аж) = 0. 

Since the radius vector of the current point of the osculat- 
ing plane p = zi + yj + zk, it follows that by passing 
to coordinate notation we obtain an equation of the 
desired plane in the form hy — az = 0. 

Formulas expressing the derivatives of the vectors 
10, b°, n? are called Frenet formulas: 


dt _ o 4% „ 1 о dn — іо í po 
—— = п", — — n = g'-T,9M 


1 
ds R ds T° ds 


3-910 
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41. Write down the equation of the osculating plane at 
the point Е = 2 of the curve 

тей + ИК. 


42. Write down the equation of the osculating plane at 
the point t = 0 of the curve 


r — eli 4 ej + V 2 tk. 
43. Find the torsion at the point і = 0 of the curve 
r = е! cos t-i + e! sin t-j + e'k. 
44. Find the torsion at any point ¢ of the curve 
г = a cosh t-i + a sinh t-j + atk. 


CHAPTER II 
SCALAR FIELDS 


Sec. 7. Examples of scalar fields. 
Level surfaces and level lines 


Definition. If a value of a certain quantity is defined 
at every point of space or a portion of space, then we say 
that the field of the given quantity has been specified. ~ 

The field is termed a scalar field if the quantity in 
question is a scalar quantity, that is, if it is fully de- 
scribed by its numerical value.-: 

Examples of scalar fields are: a temperature field, an 
electrostatic field. 

Specifying a scalar field is accomplished by specifying 
the scalar function of a point M: 


и = f (M). 
If a Cartesian coordinate system zyz is introduced in 
space, we have 
u = f (x, у, 2). 
Geometrically, a scalar field is characterized by a 
level surface; this is a locus of points at which the scalar 


function of the field assumes the same value. The level 
surface of a given field is defined by the equation 


f(z, у, 2) = C, where С = constant. 


In the case of a temperature field set up in a homo- 
geneous and isotropic medium by a point source of heat, 
the level surfaces are spheres centred at the source (this 
is a central-symmetric field). 


3* 
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In the case of an infinite uniformly heated wire, the 
level surfaces (isothermic surfaces) are circular cylinders 
whose axes coincide with that of the wire. 

Example 1. Construct the level surfaces of the scalar 


field 
и = х + 2y + 32. 
Solution. The level surfaces are given by the equation 


z + 2y + 32 = C, where С = constant. 
This is a one-parameter family of parallel planes, 
Example 2. Find the level surfaces of the scalar field 
и = r + у? — 22. 
Solution. The level surfaces are given by the equation 
z?4-y?— 22 = C, where С = constant. 
For C = 0, we obtain a circular cone. For any C >0, we 
obtain a hyperboloid of revolution of one sheet with the 
axis coincident with the z-axis. For С < 0, we obtain 


a hyperboloid of revolution of two sheets. 
Example 3. Find the level surfaces of the scalar field 


2 
и = arcsin Per . 


Solution. The domain of definition of the given scalar 
field is found from the inequality 


Б that is, Oct —— <1, 


worst, 


| vate , 


whence 0 < 22 < z? + у. This double inequality shows 
that the field is defined outside a circular cone z? — 
= z? + у? and on it, with the exception of its vertex 
O (0, 0, 0). 

The level surfaces are found from the equation 


arcsin =C, where — «06x. 


2 
y zs 
That is, z/V z? + y? = sin C or z* = (2° + y?) sin? C. 
This is a family of circular cones located outside the cone 
z? = 1? + y? with a common axis of symmetry, Oz, and 
a common vertex, O (0, 0, 0), at which the given field 
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is not defined; note that the cone itself, 22 = z? + у, 
is also included in the family. 
Example 4. Find the level surfaces of the scalar field 


и = ее, г), 
where а is а constant vector and г is the radius vector 


of a point. 
Solution. Here, 


г = (z, у, 2} = zi + yj + 2k 
and let 
а = {a;, аз, аз} = ай + аз + ask 

Then the scalar product 

(а, г) = az + asy + azz. 
The equation of the level surfaces is 

ee. n—C, C0, 
whence 
(a, г) = In С 

or 

аут + ayy + аз: = 1а С. 
This is a family of parallel planes. 

Find the level surfaces of the following scalar fields: 


23 yt zg 
45. и= 7+5 +46. 
46. и= 224+ y?— 2. 
47. ponte. 


2 


48. и=2у? + 922. 
49. и = 3**2v-*, 
(а, г) 

50. и= (Б, г) (а, b are constant vectors). 
51. и = In |г!. 
52. и = ea. ^n (a, b are constant vectors). 

A scalar field is said to be piane if there is a plane 
such that in all planes parallel to the given plane the 
scalar field is the same. 
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If we take this plane as the ry-plane, then the scalar 
field is given by the scalar function 


и = f (z, у), 


which is to say it is not dependent оп 2. 
An example of a plane scalar field is the temperature 
field of an infinite, uniformly heated wire. 


ух у VI 
(«0 
6>0 6>0 
I 
C«0 
Fig. 11 


Geometrically, plane scalar fields are characterized 
by level lines; these are loci in which the scalar function 
has one and the same value. 

Example 5. Find the level lines of the scalar field 


и = 23 — y}. 
Solution. The level lines are given by the equations 
a? — y? = С, where С = constant. 
When C = 0, we obtain a pair of straight lines: 
у = =, = —т. 


For С 0, we obtain a family of hyperbolas (Fig. 11). 
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Find the level lines of the following plane fields: 
53. и= 2z — y. 


54. ив И 2. 
2 


55. u= 2., 
r 

56. u= ех'–и. 

57. Find the level lines of the scalar field u given im- 
plicitly by the equation 


и + тши- у=0. 


Sec. 8. Directional derivative 


Suppose we have a scalar field defined by a scalar 
function 
и =f (M). 


In the field, take a point M, and choose a direction 
indicated by the vector 1. Then in the field take another 
point M so that the vector M,M is parallel to 1. Denote 
by Au the difference 


Ди = f (M) — f (Mo) 


and by Al the length of the vector М.М. The ratio Au/Al 
defines the average rate of change of the scalar field per. 
unit of length' in the given direction. Allow the point M 
to move towards the point M, so that the vector М.М 
is always collinear with the vector 1. Then Al — 0. 

Definition. If, as Al — 0, there is a limit to the ratio 
Au/Al, then it is called the derivative of the function 
и = f (M) at the given point M, in the direction of 1 
and is denoted by the symbol 8u/ól so that, by definition, 
we have 


2 — lim Ar = lim LORE) мом |. 


м-0 M м-0 

This definition of а directional derivative is invariant, 
that is, it is not connected with any choice of coordinate 
system. 
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Suppose a Cartesian coordinate system has been intro- 
duced in space and suppose the function f (М) = 
= f (z, y, 2) is differentiable at the point М, (Zo, yo, Zo). 
Then 


ди = 


ðu У ди ди 
© |, = 52-1, cosa + эг , cos B+ 5; бозу» (1) 


M M, 


where cosa, cos B, cos y, the direction cosines of the 
vector 


1 = ай + a,j + ask, 
are found from the formulas 


UH. n NL. З 
cosa = 71} соѕ В ПГ, sys’ 
Ш = V a? - а? + a3. 
ди ди ди uses 
те symbols Be luc ay ih Oe |ы, signify that the par- 
tial derivatives are taken at the point Mo. 
For a plane field u = f (z, y), the directional derivative 
l at the point M, (zo, yo) is 


ди ди 


ЭГ |м, дт 


ди x 
cosa + ——| sina 2 
Mo + ду lu. , (2) 


where œ is the angle formed by the vector 1 and the 
z-axis. 

Remark. The partial derivatives ди/дт, ди/ду, du/dz 
themselves are derivatives of the function и in the direc- 
tion of the coordinate axes Ох, Oy, Oz respectively. 

Formula (1)—used to compute the directional derivative 
at a given point—holds true even when the point М 
tends to M, along a curve for which the vector l is the 
tangent line at the point Mj. 

Example 1. Find the derivative of the scalar field 


и = хуг 
at the point M, (1, —1, 1) in the direction from М, to 
M, (2, 3, 1). 

Solution. We find the direction cosines of the vector 
MM, = (1, 4, 0), the length of which is |MoM,| = 17, 
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and we have 


cosa = cos ут, cos y = 0. 


—. 
И17 ' 
The values of the partial derivatives of the function 
и = zyz at the point M, (1, —1, 1) are 
8| __{ ди) ua дш 
“Oz |м, — ' ду]м, ^ дт 
Using formula (1), we get 


ди 1 


м, у 


4 3 

a yu уп 5 уш? 
ди 

The fact that aE |m 


M, increases in the given direction. 
Example 2. Compute the derivative of the scalar field 


= — 1. 


Mo 


> 0 means that the scalar field at 


u = arctan ry 


at the point M, (1,1), which belongs to the parabola 
у = 2’, in the direction of the curve (in the direction of 
increasing abscissas). 

Solution. The direction of 1 of the parabola у = z? at 
the point M, (1, 1) is the direction of the tangent to the 
parabola at that point (Fig. 12). 

Suppose the tangent l to the curve at M, forms with 
the z-axis an angle a. We then have 


y = 2z, tana = у |2, = 2, 


whence the direction cosines of the tangent line are 
1 1 


cos 2 = — = — 
y 1-Ftan* с 5' 


cos B = sina 


_ tana _ 2 
У1- а y5' 


The values of the partial derivatives of the given function 
и (z, y) at the point М, (1, 1) are 


1 1 ди 1 1 


ди 
M, 1+z%yilm, 2' бум, 12мм, 2’ 


“Or 
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Substituting them into (2), we obtain 
ди 1 1 1 2 


WU Е EP 235. 
al^ 2 V5 2 y5 2y5' 
Example 3. Find the derivative of the scalar field 


Fig. 12 


и = zz? + 2yz at the point М, (1, 0, 2) along the circle 
х=1+ = 


y=sint—1, 
2=2. 


Solution. The vector equation of the circle is of the 
torm 


r (t) = (1 + cos t) i + (sin t — 1) j + 2k. 
We find the vector t tangent to it at any point M to be 
5 = —sint-.i+cost-j. 
The given point M, (1, 0, 2) is found in the zz-plane in 


the first octant and is associated with the value of the 
parameter { = n/2. At this point we have 


DINE c. л. : 
тм, = — sin > -i+ cos 3 j--i 
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From this we obtain that the direction cosines of the 
tangent to the circle are equal to cos а = —1, cos В = 0, 
cos ү = 0. The values of the partial derivatives of the 
given scalar field at the point M, (1, 0, 2) are: 


ди -2| EN NL =2z| =4, 
дг |м, м, ду IM, М. 
ди 
E м, = (222 + 2y) My 


Hence the desired derivative is 


ди 


ди_ _ ди 
д1 


M, ov |M 


=4.(—1)+4.044.0=—4. 


In the following problems, Ц is required to find, for 
the given functions, the derivative at the point 
М, (хо, Yo: 20) in the direction of the point M, (21, Yi» 2). 


Ж. i u-yz3 yz, Mq. 1,1), М, (3,2, 1). 
597 um 12у + 122—9, М. (1,1, —1), М, (2, —1, 3). 
„б. и = хе” + ye* — 22, М, (3, 0, 2), М, (4, 1, 3). 


61. и= 2—6, М о (1, 1), М, (4, 5). 
62. Find the derivative of the scalar field 
u = In (z? + y?) 


at the point М, (1, 2) of the parabola y? = 4z in the 
direction of the curve. 

63. Find the derivative of the scalar field u = arctan y/z 
at the point М, (2, —2) of the circle z? + у? — 45 = 0 
along an arc of the circle. 

64. Find the derivative of the scalar field и = z? + у? 
at the point М, (Zo, yo) of the circle z? + y? = R? in the 
direction of the circle. 

65. Find the derivative of the scalar field и = 2zy + y? 
at the point (V 2, 1) of the ellipse 22/4 + y?/2 = 1 in the 
direction of the outer normal to the ellipse at that point. 
66. Find the derivative of the scalar field и = 2? — у? 
at the point (5, 4) of the hyperbola z? — y? = 9 in the 
direction of thé curve. 


2^ 
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67. Find the derivative of the scalar field и = 
= In (xy + yz + zz) at the point М, (0, 1, 1) in the 
direction of the circle z — cos t, y — sin t, z — 1. 

68. Find the derivative of the scalar field u = z? + y? + 
+ 22 at the point M, that corresponds to the value of 
the parameter t = л/2 in the direction of the helical 
curve т = R cost, у = R sin t, 2 = at. 


Sec. 9. The gradient of a scalar field 


Suppose we have a scalar field defined by a scalar 
function 
u=f (z, Y, 2) 
where the function f is assumed to be differentiable. 
Definition. The gradient of a scalar field u at a given 
point М is a vector denoted by the symbol grad и and 
defined by the een 


grad u- Si = +k. (р) 


Using formula (1) of Sec. 8 for ‘he directional derivative, 
we have 


= (grad и, 1°), (2) 
where 1° is а unit vector in the direction of l, that is, 


M = = 10059 +1005 В+ К соз y. 


Properties of a gradient 


- 1. The gradient is in the direction of the normal to the 
level surface (or to the level line if the field is a plane 
field). 

2. The gradient is in the direction of increasing values 
of the function of the field. 

3. The modulus of the gradient is equal to the largest 
directional derivative at a given point of the field: 


max 57-:grad u| = и (2 Ay tm ) + (52 A). 


* The maximum is taken over all directions of l at the given 
point of the field. 


* 
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These properties yield an invariant characteristic of 
the gradient. They state that the vector grad u indicates 
the direction and magnitude of maximum change of 
a scalar field at a given point. 

Example 1. Find the gradient of the scalar field 


u = т — 2y + 3z. 
Solution. By (1) we have 
grad u = 1i — 2j + 3k. 


The level surfaces of the given scalar field are the planes 
z — 2y + 3z = C; the vector grad u = (1, —2, 3) is the 
normal vector of the planes of this family. 

Example 2. Find the greatest steepness (rate) of rise 
of the surface и = x” at the point ‚2,4 

Solution. We have 
grad и = yzV li + z” ln zj, gradu |м = 4i + 4 In 2j, 

GE). = grad u| = 4 V 1 4- (In 2)?. 

Example 3. Find the unit vector of the normal to the 

level surface of the scalar field 


и = 22 + yt 2. 
Solution. The level surfaces of the given scalar field 
are the spheres 


Y+y+2=C (С> 0). 


The gradient is directed along the normal to the level 
surface so that grad и = 2z-i + 2y-j + 2z-k defines the 
vector of the normal to the level surface at the point 
М (т, у, 2). For the unit vector of the normal, we obtain 
the expression 

gradu _ ук — г 
[ered ui Узи — 


ir} ^ 
Example 4. Find the gradient of the field и = (a, b, г), 
where a and b are constant vectors and r is the radius 
vector of the point. 
Solution. Let 


по = 


а = (a4, аз, аз}, b = (b, b, 03), т = {z, у, 2}. 
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Then 
а az аз 
и = b, bz b; А 
гу 2 
By the rule for differentiating а determinant" we have 
04 а аз 04 а, аз 0, а, аз 
a= b, b, 6з |, a by b bs |, M = by bz bs |. 
100 01410 001 
Hence 
аз 03|, |@1 аз а; 0 
ыы quis bs 
ij k 
= |а, а, аз | == [a, Б). 
b, by b; 


* Given a determinant D (t) whose elements a,; are differenti- 
able functions of t: 


arı (t) @1з(1) ... an (1) 
р (= азі (t) а, (f). ... Gen (t) 


ami (i ana (D -.- ann (0 


Then the derivative of the determinant, D' (t), is found from the 
formula 


aj (t) aig (0. ... Gin (1) 
D'(t)=| “1 (0. Gan 0. «++ Gan (t) 


апі (t) апз (1) -.. ânn (t) 
G1 (t) азаб)... ain (1) 
T Gs (t) ass (t)... а (0) 


ох 


апі (f) ans(t) ..- Gnn (0) 
Gy; (t) aig (t) ... ain (0) 
... +| 9n (@) аз)... asn (0) 


anı (t) analt) ... ann (t) 
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Example 5. Find the gradient of the distance 


where Р (z, у, 2) is the point of the field being studied 
and Р, (х0, yo, Zo) is some fixed point. 
Solution. We have 


Or . , Or or 
gradr = z ic 3-]t о Е 


—&—0)1+у—ю)1+@— К _ 
V (z—z9* + (у— Yo)? + 6— 20)? 
which is the unit vector of the direction P,P. 


1 


Fig. 13 


Example 6. Let us consider the scalar function 
= n + Г, 


where гу, г, are the distances of some point Р (т, y) of the 
plane from two fixed points, F, and F,, of the plane. 
Solution. The level lines of this function are ellipses. 
We have (see example 5) 
grad (г, + г.) = r? + го. 
This shows that the gradient is equal to the diagonal of a 
rhombus constructed on the unit vectors of the radius 


vectors drawn to point P from the foci F, and F, (Fig. 13). 
Consequently, the normal to the ellipse at some point 
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bisects the angle between the radius vectors drawn to that 
point. 

Physical interpretation: a light ray coming from one 
focus enters the other focus. 

Example 7. Find the angle Ө between the gradients of 
the functions 


и=У 22 у? and v=r+y+2V zy 


at the point M, (1, 1). 
Solution. We find the gradients of the given functions 
at the point M, (1, 1): 


_ m | _ iti 
втейи |м, = Vr | Yi 


grad v |n, — [ (1-- V/ 2) 1+ (14- =) j]],, - 2i-2i. 


The angle 0 between grad u and grad v at the point M, 
is found from 


2 2 
ЕТА (gradu, grad v) |y, — yat V2 24 
I gradu 1м, -|gradviy,  — 4.2V2 — — 
From this we have 
0 = 0. 
Example 8. Find the directional derivative of the 
radius vector г for the function и = sin г, where г = |r |. 


Solution. By (2), the directional derivative of the 
given function of the radius vector r is 


ди Я 
z = (grad зїп г, го). (3) 
We find the gradient of the function: 

]t——— 


hn Or. d(sinr) дг 
а Tt wr 


д en 9 (зір г). k 


grad sinr = 2682 ay Ce I 


__ d(sinr) дг. 
— dr дт | 


=(жї+ү itg Z k) cosr = r^ cosr. (4) 
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Substituting (4) into (3), we get 
E = (r° cos г, r°) = (го, r°) cos Г =cosr. 

Example 9. Find the derivative of the scalar field 
и = f (z, у, 2) at the point M, (zo, yo, zo) of the curve І 
specified by the system of equations 

Í (z, y, 2) =а, 
Ф (2, y, 2) =0 
in the direction of the curve. 

Solution. The direction of the curve l is given by the 
direction of its tangent vector t, which, by definition, 
із a vector that is tangent to the surface f (z, y, 2) = a. 
The surface f (z, y, z) — a is a level surface of the given 
scalar field u — f (z, y, z). Since 


| ; @=constant, 


дш ы (grad и, 1°) — (grad u, 1? 
д1 ) 


and the vector grad u is perpendicular to the level surface 
f (2, у, 2) = a, it follows that grad и is perpendicular to 
the unit vector то, and therefore 

ди 

ol |м, 

Example 10. At the point M, (1, 1, 1), find the direc- 
tion of the greatest change in the scalar field u = zy + 
+ yz + zz and the magnitude of that change. 

Solution. The direction of the greatest change of the 
field is indicated by the vector grad и (М). We find it 
thus: 

gradu (М) = (y+z)it@+2ji+yt+a2)k 
and hence grad u(M,) = 2 (i+ j+ k). This vector 
determines the direction of the greatest increase in the 


given field at the point M, (1, 1, 1). The magnitude of 
maximum change of the field at this point is 


шах 27 = | гад и (Мо) | —2y 3. 


69. Find the gradient of the scalar field и = ln (z? + 
+y? + z?) at the point M, (1, 1, —1). 

70. Find the gradient of the scalar field и = zex**v**:? at 
the point O (0, 0, 0). 


4-910 


= (grad и, T°) |m, — 0. 
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71. Find the angle ф between the gradients of the func- 
tion u — arctan z/y at the points М, (1, 1) and 
М. (—1, —1). 
72. Find the angle @ between the gradients of the func- 
tion и = (z + y) e*** at the points М, (0, 0) and 
M, (1, 1). 
73. Find the angle q between the gradients of the func- 
tions и = V2? + y? + 2 and v = In (22 + y? + 2) at 
the point M, (0, 0, 1). 
74. Find the points at which the gradient of the scalar 
field и = sin (х + y) is equal to i + j. 
75. Find the points at which the modulus of the gradient 
of the scalar fieldu = In V z? + у? + 2? is equal to unity. 
76. Letu = и (z, y, 2) and v = v (т, у, 2) be functions 
differentiable at the point M (z, y, z). Show that 

(a) grad (Au) = A grad u, A = constant; 

(b) grad (u + v) = grad u + grad v; 

(c) grad (uv) = v grad u + u grad v; 


(d) grad (=) =ou eee, v 3-0. 
77. Show that 
grad u(q) = grad Фф, 


where Фф = Фф (z, y, 2) is a differentiable function and 
и = и (Ф) has a derivative with respect to Ф. 
Find the gradients of the following scalar fields if 


r=zit+yjtzk, г=|г| = У 22 + у 22, 


and a and Ь are constant vectors. 


78. и = lnr. 

79. и = (a, г). 

80. u = (a, r)-(b, г). 
81. и = | [a, rl f. 


82. Show that 
(grad и (г), г) = u’ (r)-r. 
83. Show that 
[grad и (г), г] = 0. 


84. Let w = f (и, и), where и = и (z, у, 2), v = v (т, у, 2). 
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Prove that 

grad w = a gradu + СА grad v 
if f, u, v are differentiable functions. 
85. Suppose G is a convex region in space (that is, a 
region such that if two points M and N belong to G, 
then the whole line segment MN lies in G). Let there be 
given in G a scalar field u (M) which at all points has a 
gradient that is continuous and bounded in G: 


grad u (М) | < A, MEG, A = constant. 


Prove that for any points M and N of G we have the 
inequality 
[и (№) —u(M)|<A|MN|. 


86. Find the derivative of the function и = z?/a? + 
+ y?/b? + z?/ at an arbitrary point M (т, у, z) in the 
direction of the radius vector r of that point. 
87. Find the derivative of the function u = 1/r, where 
г = |r| in the direction of the vector l = cos a-i + 
+ cos B.j + cos y-k. Under what condition is the 
derivative equal to zero? 
88. Find the derivative of the function и = 1/r, where 
r — |r |, in the direction of its gradient. 
89. Find the derivative of the function u — yze* at the 
point M, (0, 0, 1) in the direction of its gradient. 
90. Find the derivative of the scalar field 
u = u (z, у, 2) 
in the direction of the gradient of the scalar field 
v = v (z, у, 2). 

Under what condition is it equal to zero? 
91. For the following scalar fields, find the direction and 
magnitude of greatest change at the given points Мо: 

(a) u (M) = ЖА ыды Ka (1, 0, 0). 

(b) u (M) = zyz; M, (2, 1 


4* 


CHAPTER III 
VECTOR FIELDS 


Sec. 10. Vector lines. Differential equations 
of vector lines 


Definition 1. We say that a vector field is given if a 
vector quantity a = a (M) is specified at each point M 
of space or of a portion of space. 

If a Cartesian coordinate system is introduced in the 
space, then specifying the vector field a= a (M) is equiv- 
alent to specifying three scalar functions of the point 
P (M), Q (M), R (M) so that 


а (M) = Р (z, y, 2) i + Q (z, у, 2) j + R (z, y, 2) К. 


Definition 2. A vector line of a vector field a is a 
curve at each point M of which the vector a is directed 
along the tangent to the curve. - 

Let a vector field be defined by the vector 

a = Pi + Qj + Rk, 
where 
P =P (z, У, 2), Q = 0 (z, У, 2), R =R (=, y, 2) 
are continuous functions of z, y, z that have bounded par- 
tial derivatives of the first order. 

Then the differential equations of the vector lines are 
of the form 

de dy ds (4) 

Integrating this system of two differential equations 

(1) yields a system of two finite equations: 


Фа (т, у, 2) = Cy, Ф, (2, у, 2) = С, 
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which, taken together, define a two-parameter family of 


vector lines: 
, ?, =C 1 
Ф, (2, у xm 1 ! (2) 
Ф (T, y, 2) = Cs. 
If the conditions of the theorem of existence and 


uniqueness of a solution are fulfilled in a certain region 
G for system (1), then a unique vector line 


Ф; (т, у, 2) = Фи (20, уо, 20), | 
Ф (2, у, 2) = 9» (20, Vo: 20) 
passes through every point М, (Zo, Yo: Zo) € G. 
Example 1. Find the vector lines of the vector field 
а = le, г], 
where c is a constant vector. 
Solution. We have 
e= ci + cj t+ cesk, r = zi + yj + zk 
so that 
ijk 
а = [с, г] = |с; с, сз 
ry z| 
= (C22 — сзу) i + (esz — 0,2) } + (ey — Сох) К. 
The differential equations of the vector lines are 
aa = ais = SOR £ (3) 


Multiply the numerator and denominator of the first 
fraction by т, the second by y, the third by z and add 
termwise. Using a property of proportions, we have 

dz dy dz zdz+ydy+zdz 

севу вая сузу 0 7 
whence 
д ах + ydy + 242 =0 

and this means that 

r+y+2=A,, A, = constant > 0. 
Now, multiplying the numerator and denominator of the 
first fraction of (3) by c, the second by с», the third by c; 
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and adding termwise, we obtain 
dz dy dz € dz+-¢, dy +e, dz 
0 


= = 
632 — сзу Суг — 611 Ciy — Сэ; 
whence 


су dz + с. dy + сз dz = 0 
and, consequently, 
сух + cay + 032 = Az, А, = constant. 
The desired equations of the vector lines are: 
P+y+2=A;, | 
E CiT + суу + C32 = As. 

These equations show that the vector lines are obtained 
via the intersection of spheres (having a common centre ut 
the origin of coordinates) 
with planes perpen- 
dicular to the vector е = 
=¢i + cj + ck. From 
this it follows that the 
vector lines are circles 
whose centres lie on a 
straight line passing 
through the coordinate 
origin in the direction of 
the vector c. The planes 
of the circles are per- 
pendicular to the indicat- 
ed straight line (Fig. 
14). 
Fig. 14 Example 2. Find the 

vector line of the field 


a = —yi + zj + bk, 


which line passes through the point (1, 0, 0). 
Solution. The differential equations of the vector 
lines are 


whence we find 
% + у® = cy 6-0 
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or, introducing the parameter ż, 
z=VC,cost, у= V C,sint. 


In this case, the equation 


takes the form 
VC,costdt — dz E 
Исен e OREST 


whence we obtain 


Thus, the parametric equations of the vector lines are 
r= УС, соз t, 
у= V C,sint, (4) 
$ = bt + Ca. 
If we require that the vector line pass through the point 
(1, 0, 0), we will have 
1=VC,cost, 0—VC,sint, 0-bt4-C,. 
The first two equations of this system are satisfied for 
t = 2kn, k = 0, +1, . . .and for C, = 1. Taking k = 0, 
we get t = 0 and the last equation of the system yields 
C, = 0. The desired vector line passing through the 


point (1, 0, 0) is 
zr=cost, 
y=sint, | 


2= bt. 

This is a helical curve. 

Find the vector lines of the following vector fields: 
92. r = zi + yj + zk. 
93. a = ai + a,j + ask, where a, а», аз are constants. 
94. а = (z — y) i + (z — z) j + (y — 2) k. 
95. Find the vector line of the field 

а = ri — yj + 2k, 


which line passes through the point (?/,, —!/,, 1). 
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A vector field is said to be plane if all the vectors of a 
are located in parallel planes and the field is the same 
in each of the planes. 

If a Cartesian coordinate system zOy is introduced in 
one of the planes, the vectors of the field will not contain 
any components along the z-axis and the coordinates of 
a vector will be independent of z, that is, 


a = Р (z, y) i + Q (z, у)}. 
The differential equations of the vector lines of a 


plane field are of the form 


dz dy _ 4 
P(z,y) Q(z, y) 0 


or 
dy _ Q(z, y) 
ат  P(zy' 
z= constant. 


From this it is evident that the vector lines of a plane 
field are plane curves lying in planes parallel to the 
zy-plane. 

Example 3. Find the vector lines of a magnetic field 
of an infinite current conductor. 

Solution. We will assume the conductor is in the direc- 
tion of the z-axis and that the current J flows in that 
direction. The intensity vector H of the magnetic field 
set up by the current is 


Н-П, г], (5) 


where I = Г.К is the current vector, г is the radius vector 
of the point M (z, у, т), and р is the distance from the 
axis of the wire to the point M. Expanding the vector 
product (5), we obtain 


H- — 2 21: 


or ita). 


The differential equations of the vector lines are 
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whence 
22 +02 = R, | 
2= С. 
That is, the vector lines are circles with centres on the 
z-axis (Fig. 15). 
Find the vector lines of the following plane vector 
fields: 


96. a = zi + 2yj. 
97. a = zi + zk. 
98. a = zi — yj. 
99. a = 22] + 4yk. 
100. а = z?i + yj. 
101. a = 2j — yk. 


The differential equations of 
the vector lines 


may be written as 


dz _ dy _ dz _ 
d ^P и =9, GaR 


or, in vector form, as 
47 a (M). (6) Fig. 15 


This form of the equations of vector lines turns out to be 
convenient in the solution of a number of problems. 
Example 4. Find the vector lines of the field a = [e, r), 
where с is a constant vector. 
Solution. Applying (6), we get 
dr 2 
dt = [е, rJ. (7) 
Forming the scalar product of both sides of (7) by € and 
using the properties of a mixed product, we find 


(s. x) = (e, r)=0. (8) 
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Similarly, forming the scalar product of both sides of (7) 
by r, we obtain 
(г, I) 9 (к, 9-0. (9) 
From equation (8) it follows that 
(с, г) = constant 
and from equation (9) it follows that 
(r, r) — constant. 


The vector lines are lines of intersection of the planes 
(c, r) = constant with the spheres r? = constant. 
Find the vector lines of the following vector fields: 
102. а = f (r)-r. 
103. а = (a, г) bo, where ao, b, are constant vectors. 


Sec. 11. The flux of a vector field. 
Methods of calculating flux 


J. The flux of a vector field. Suppose we have a vector 

field 

а (М) = P (z, y, 1+ 0(т, y, 2) j + R (z, y, z) К, 
where the coordinates P (z, y, 2), О (z, y, 2), В (z, y, 2) 
of the vector a (M) are continuous [the field a (M) is 
continuous] in some region G. Let S be a smooth or piece- 
wise smooth two-sided surface in which a definite side 
has been chosen (an oriented surface). 

Definition. The flux II of a vector field a (M) through 
an oriented surface S is defined as the surface integral 
of the first kind, over the surface S, of the projection 
of the vector a (M) by the normal n (M) to that surface: 


П = JJ proa dS = j (a, n°) 45, 


where n° is the unit vector of the normal n to the chosen 
side of the surface S; dS is the area element of the sur- 
face. 5. 

In the case of a closed surface, we will always choose 
the outer normal n that is directed outwards from the 
region bounded by the surface S. 
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If a, B, y are the angles that the normal n forms with 
the coordinate axes Oz, Oy, Oz to the surface S, then 


the flux may be expressed in terms of a surface integral 
of the second kind: 


П = M (a, n?) 45 = M [P (z, y, z) cosa 


+ Q(z, у, 2) cos + R (z, у, 2) cos y] dS 


or 
П = (а, по) 45 = P (x, y, 2) ду д2 
уена] 
+Q (z, y, =) dz 42+ R (т, у, 2) dz dy, 
where 


cosa dS = dy dz, cos B dS = dz dz, cos y dS = dz dy. 


Basic properties of the flux of a vector field 


(a) The flux reverses sign when the orientation of the 
surface is changed (that is, when the orientation of the 
normal n to the surface S is changed): 


| | а, по) 48 = — | | (а, nas, 


5+ S^ 


where 5+ is the side of.the surface 5 on which the normal 
n is chosen, and S- is the side of S on which the normal 
—n is taken (see [7]). 

(b) Linearity: 


1 (àa + pb, n?) dS = А Jj (a, n°) 45 + р jj (b, n?) 25, 


where À and p are constant numbers. 

(c) Additivity: if the surface S consists of several 
smooth parts S,, Sz, ..., Sm, then the flux of the vector 
field a (M) through S is equal to the sum of the fluxes 
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of the vector a (M) through the surfaces $,, $,,..., Sm: 
m 
n= > | | (а, n9as. 
kæmi "Sy 


This property permits extending the notion of flux 
to piecewise smooth surfaces. 


Fig. 16 


Example 1. Find the flux of the vector a — i through 
an area perpendicular to the z-axis and having the shape 
of a rectangle with sides 7 and 2 (Fig. 16) in the positive 
direction of the z-axis. 

Solution. According to the definition of the flux of 
a vector through a surface S, we have 


n- | (a, n°) dS. 


In our case, а = i, n° = i so that (a, n?) = (i, i) = 1. 
Taking into account that the area of the rectangle is 
equal to 2, we obtain 


A a 


Remark. If we had chosen the unit vector of the normal 
to the area S so that n? = —i, we would have got II = —2. 
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Example 2. Compute the flux of the vector field 
a — r, where r is the radius vector, through a right cir- 
cular cylinder of altitude k, base 
radius R and the z-axis. 

Solution. The surface S consists 
of a lateral surface o,, an upper 
base с, and a lower base о; of the 
cylinder. By the additivity proper- 
ty, the desired flux П is equal to 
П = П, + II, + Hs, where П,, П,, 
П, are the fluxes of the given field 
through O}, Oz, Os respectively. 

On the lateral surface о, of the cyl- 
inder the outer normal n° is paral- 
lel to the zy-plane and therefore 


(a, n°) = (г, 00) = praor = R 


Fig. 17). Hence 
П, = j | (ann) 45 = Н | f dS = R- 2n Rh = 2n R?h. 
g 


On the upper base o, the normal n? is parallel to the 
z-axis and therefore we can put n? — k (see Fig. 17). 
Then 

(a, n?) = (г, К) = pro, r = k 
and so 
П, = È Í (а, по) 25 =h ( È dS = л.да: лд. 
jene] 

On the lower base оу the vector a = r is perpendicular 
to the normal n? = —k. Therefore (a, n?) = (г, —k) = 0 
and 


n, = f f (а, n?) dS — 0. 


The desired flux is then 


n- $ (a, n°) dS = 3n R?h. 
S 
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Example 3. Find the flux of the vector field 
т 
a= те 

through a sphere of radius R with centre at the origin 
of coordinates. 

Solution. Since the normal n to the sphere is collinear 
with the radius vector r, we can take n? = г? = r/|r |. 
Therefore 


(et ee wo 
(а, по) = (тр rm) 7n P Tee = те. 
But on the sphere 5 we have |r | = R, and so (a, n?) = 


= 1/R?. 
The desired flux II is 


n= (a, п) dS = ar (b dS = 4n 
5 5 


since the area of the whole sphere 5 is equal to (as = 


S 
= 4лА:. 
104. Compute the flux of the vector a — 3j through an 
area having the shape of a triangle with vertices at the 
points M, (1, 2, 0), M, (0, 2, 0), M; (0, 2, 2) in the direc- 
tion of the coordinate origin. 
105. Find the flux of the vector 


a = ai + fj + yk, 


where a, В, ү are constants, through an area perpendicular 
to the z-axis and having the shape of a circle of radius R, 
in the positive direction of the z-axis. 
106. Find the flux of the vector a = r through the outer 
side of a circular cone whose vertex lies at the origin of 
coordinates; the base radius is equal to R and the altitude 
is h (the axis of the cone is along the z-axis). 
107. Find the flux of the vector a = f (|r |) г through 
a sphere of radius R with centre at the coordinate origin. 
II. Methods of computing the flux of a vector. 
1°. The method of projection onto one of the coordinate 
planes. Let an open surface $ be projected one-to-one 
onto the zy-plane into a region D,,. In this case, the 
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surface $ may be given by the equation z = f (х, y) and 
since the area element dS of the surface is 


dz dy 
[ сөзү], 


dS = 


it follows that computing the flux II through the chosen 
side of the surface $ reduces to computing a double integral 
via the formula 


n- Jie stas || (a, n0) аху. (1) 
S ry 


1 cos y | |2=Их, у) 


Here the unit vector n? of the normal to the chosen side 
of the surface 5 is found from the formula 


d [z— / (=, y)) dx ду 
n? = + gra = + ‹ 
| grad [z— f (=, y)] | 91 \2 тәу," 
И (ZHE +! 
(2) 
and cos y is equal to the coefficient of the unit vector К 
in formula (2): 


1 
cos y = + -————. (3) 
Of \2 of \2 
У (i) +(4,) +4 
If the angle y between the z-axis and the normal n° is 


acute, then in (2) and (3) the plus sign is taken, and if 
the angle y is obtuse, the minus sign is taken. The symbol 


(a, n9) 
| cos y | |z=f(x, у) 


signifies that f (z, y) must be substituted for z in the 
integrand. 

If it appears to be convenient to project the surface S 
onto the yz- and zz-planes, use is made of the following 
formulas to compute the flux II: 


(а, n9) 


= | cosa | 
Wt 


dy dz (4) 


х=Фф(у, 2) 
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or 


T. (a, n9) 
II == M ses T ene, o 2747 (5) 
x2 

Formula (4) is used when the surface S is projected 
one-to-one into the region D у, of the yz-plane, which means 
that it may be given by the equation z = 9 (y, 2); cosa 
is found as the coefficient of the unit vector i in the for- 
mula 


d [z— 9 (y, 2)] EE az 
ое 120 саан B 787: у шу RE 
| grad [z— Ф (y, 2)) | и дф \2 дф \2 
1+ (2%) * Ux) 
y д2 
That is, 


1 
cosa=+ 7/7 op \2, | Opie é 

V 1+ (32) +(2) 

The plus sign is taken if the angle a, between the z-axis 
and the normal n? is acute, and the minus sign if the 
angle œ is obtuse. 

Formula (5) is used in the case of a one-to-one projec- 
tion of the surface S onto the zz-plane; in this case, 
5 may be specified by the equation y = № (z, 2) and 
then 
ap дү 


grad [у (2, 2] _ Page 5er 


— А; 
| grad [y — 6 (=, 2)] | И i (yy 2 ' 


соз В is equal to the coefficient of the unit vector j in 
this formula, that is, 


cos Bp = + 


n? = + 


1 
„(н УР 
V «(Xy GE) 
If the angle В between the y-axis and the normal n° is 
acute, we take the plus sign, and if the angle В is obtuse, 
we take the minus sign. 


Remark. When the surface 5 is specified implicitly by 
the equation Ф (z, y, 2) = 0, the unit vector of the 
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normal 
n? = i соза + j соз В + К cos y 


is found from the formula 


0 grad Ф (z, у, 2) 
п == | grad Ф (z, у, 2) | 


9Ф aD 
de ita 


aD \2 | @ \2, ү 2’ 
и (2) tU) (ae) 
where the sign on the right is determined by the choice 
of the normal to the surface S. 

To compute the flux П of a vector field a through a 
surface S, it is necessary to project the surface one-to-one 
onto one of the zy-, zz-, yz-planes; this is possible if the 
equation Ф (z, y, 2) = 0 is uniquely solvable with respect 
to z(z = f (z, y). y (y =% (z, z) or z(z— e (у, 2) 
respectively. Then take advantage of one of the formulas 


(1), (4), (5). 
Example 4. Find the flux of the vector field 


а = (r—2z)i + (z + 3y + 2) j + (öz + y) Е 


through the upper side of the triangle ABC with vertices 
at the points A (1, 0, 0), B (0, 1, 0), C (0, 0, 1). 

Solution. The equation of the plane in which the triangle 
ABC lies is of the form z + y + z = 1, whence z = 
= 1 — z — у. The triangle ABC is projected one-to-one 
onto the zy-plane into the region D ,,, which is the triangle 
OAB (Fig. 18). 

It is given that the normal n? to the plane in which th 
triangle ABC lies forms an acute angle y with the z-axis 
and so we take the plus sign in (2) and obtain 

о Edu iD 9 

r= Tract 07 ys үз! ys (9) 

We find the scalar product 


(a, n°) =(2—22) — + (2+ By +4) Fe + (2 +4) т 


_ Tz+4y—2z 
Уз ` 


дФ 
jk 
EE 92 
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From formula (6) we find that cos y = 1/V 3 — 0 and, 
hence, 


dS = aa = V 8 dz dy. 


Using formula (1), we po the desired flux: 


- | | (а, n°) dS = M (7z -- 4y — 2) |z=1 -x-y dz dy 
S 


хи 
1-x 


= | | (д:+5у—1)4:ау= | dz | (8z+5y—1)dy =. 
Dxy 0 


Example 5. Find the flux of the vector a = у?) + zk 
through the portion of the surface z = 2? + y? cut off by 


Fig. 18 


the plane z — 2. The outer normal is taken with respect 
to the region bounded by the paraboloid. 

Solution. The given surface (a paraboloid of revolution) 
is projected one-to-one onto the zy-plane into the circle 
D,, (Fig. 19). We find the unit vector of the normal n? 
to the surface S: 


—2zi—2yj+k 


grad (z — z* — y?) 
V 429+4y?+1 ` 


[grad (223—8) ^ = 


по = + 
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It is given that the normal n? forms an obtuse angle y 
with the z-axis and therefore the minus sign is taken in 
front of the fraction. Thus, 

пак 
У 4-41 ' 


whence 
1 
wr yu < 
and so 
dzd IZ 17.231434 
ETT = y 422+ Ay? 4-1 dz dy. 
We find the scalar product 
(a, n°) dad 2y3—z 


Иа Hi ° 
The desired flux is, by (4), equal to 
П = | f (a, n°) dS = № (902—2) |з dz dy 


5 к 


= | | (203 — y?— 12) dz dy. 


Dxy 
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The domain of integration О, у is a circle of radius А = 
= ) 2 with centre at the coordinate origin. Introducing 
polar coordinates z = р соз Фф, у = p sin g, we have 


П = | | (2p? sin? ф— p?) p dp d 
Dxy 
2n V2 
=| de j (2p* sins ф — p?) dp = — 2n £t "n — 2n. 
0 0 


Example 6. Find the flux of the vector field a = 
= į — j + ту through the circle S obtained by cutting 


NN 
N 


N 


АА 
NN 


A\ 


Z 


the sphere z? + у? + z? < А? with a plane у = =. Take 
the side of the circle facing the positive portion of the 
z-axis. 

Solution. Since the plane y = z is perpendicular to the 
zy-plane, the circle S lying in that plane is projected onto 
the zy-plane into the line segment A,A, and so the one- 
to-oneness of the projection is disrupted. The circle S 
is projected one-to-one onto the other coordinate planes. 
Projecting the circle onto the zz-plane, we obtain a region 
D,, bounded by an ellipse (Fig. 20). The equation of the 
ellipse can be found by eliminating y from the system 


See. 11] The flux of a vector field 69 


of equations 
д у? #2 = Ri, | 


у=т, 
whence 
233 pg = R or 24 a1, 
I 


It is given that the normal to the circle S forms au 
obtuse angle B with the y-axis (see Fig. 20) and so we 
take 

n= — grad (у —z) =i—j, 


From the latter equation we have cos p — —1/V 2 — 0. 
The area element dS of the circle is equal to 
45=- V3dz dz. 


= TeosB | 
We find the scalar product: (a, n°) = y 2. 
The desired flux is, using formula (5), 
= = [ = лн? == 9 2 
П | 24:41=2 || drdie 2E = УЗ Ет 
since the area О of the region D,, bounded by an ellipse 
with semi-axes a = R/|/2 and b = R is equal to 


Q= M dz dz = nab = УЗ . 


Example 7. Compute the flux of the vector a — 
= zi + yj + zk through the outer side of the lateral 
surface of the circular cylinder т? + y? = А? bounded 
by the planes 2 = 0 and z = Н (Н > 0). 

Solution. The given cylinder is projected onto the 
rzy-plane into a line, namely, into the circle (Fig. 21) 


z у? = R | 
z=0. 


zI 
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We will therefore project the cylinder onto the other 
coordinate planes, for instance, the yz-plane. Since the 
cylinder does not project one-to-one onto the yz-plane, we 
take advantage of the additivity property of the flux of 
the vector and represent the de- 
sired flux П as a sum of 
fluxes: П=П,-ЕП»,, where П, is 
the flux ofthe field through the 
portion S, of the cylinder located 
in the region where y > 0, and 
П, is the flux of that field 
through the portion S, of the 
cylinder located in the region 
where y < 0. On 5, we have 


ne = AM ‚ (a, n?) — т zte =R 


and so 
П, = RdS=R dS = RS, 
[^-^] 
Fig. 21 where 5 is the area of the por- 


tion S, of the cylinder. Since 
5 = xRH, it follows that П, = x R?H. 
On S, we again have 
_ 1+] т-у? _ 
ee dem eR. 


and so 


ll, = || В 45 = RS =n. 


5, 


The desired flux is П = 2x R?H. 

Remark. The solution is made simpler if we introduce 
curvilinear coordinates т = В соз Ф, у = Rsinq,z—z 
on the cylinder (see item 3° below). 

To find the flux of the vector field a = P (z, y, z)i + 
+ Q (z, y, z) 3 + В (z, y, 2) К through the surface 5 
specified by the equation 2 = } (т, y) by the method of 
projection onto a coordinate plane, it is not necessary to 


Sec. 11] The fluz of a vector field 71 


find the unit vector of the normal n°, and we сап take 
the vector 
д of . 
n= + атай [z— f (т, у)] = + (а 3x). 


Formula (1) for finding the flux II becomes 
П = | f (a, n°) dS = ff (a, n) dz dy. (7) 
5 Dry 


In similar fashion we obtain formulas for computing 
the fluxes through surfaces given by the equation z = 
= Фф (y, 2) or y= (=, 2). 

Formula (7) is written thus in coordinate form: 


n= + || {—Р{г, у, fie ora — 01, у, fe, У 


Dxy 


т=/(х, y) 


+ Riz, у, f(z, у} dz dy. 
Example 8. Compute the flux of the vector field 


а= у НУ -у?—1К 
through the outer side of the hyperboloid of one sheet 
z = Vz? + y! —1 bounded by the planes z = 0,2 = V3. 
Solution. The given surface is projected one-to-one onto 
the zy-plane into the region D,, bounded by the circles 
P+ y=, | and 77775 | 
z=0, z=0. 
We find the outer normal n: 
= —~Vriyp—i)= zzu 
n= + grad (2— V 2 F y2 —1)= + ( varpi +k). 
Since n forms an obtuse angle y with the z-axis (Fig. 22), 
we take the minus sign and, hence, 
äts _ 
Vapi 
We find the scalar product 


Eu yz {1...1 _. 
Mun eRpaeq oio т 
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Using formula (7), we obtain 


n= | [ (а, п) 25 | | drdy 


Z3 534—414 * 
: D. V zi y1—1 


Passing to the polar coordinates z = p cos ф and у = 
= p sing, we have 


л 
= pdpdp — pdp 
п- |} V e1—1 = | de | y p*—1 
Dry 0 1 
=2n V g—1]|; =2V 3n. 
Example 9. Compute the flux of the vector field 
a = yi + zj + zk 
through a closed surface bounded by the cylinder z? + у 
= R? and the planes z = z, z = 0 (z > 0). 


Fig. 22 


Solution. The surface S is piecewise smooth and so we 
take advantage of the additivity property of a flux re- 
presenting the desired flux П as a sum of fluxes П,, JI,, П, 
through the smooth portions S, (semicircle z? + y? < А?, 
0<z< R, 2 = 0), S, (part of the plane 2 = z), and 
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5; (part of the cylinder z? + y? = R*): П = II, + II, + 
+ II,. Since S is closed, we take the outer normal to it 
(Fig. 23). 
(1) On S,, where 2 = 0, we have n? = —k and so 
(a, n?) = —z. 
This means that the flux 


ese MET = — | 2424. 


S, S, 


Passing to the polar coordinates z = р cos ф and у = 


= р sin 9, we find 


Il, — — | | p cos pp dp dp 
+ л/2 R 
= = eee | (ар = -5 Re. 


(2) On S,, where z = z, we have 
n = +grad (z — х) = + (—i +k) 
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and since the normal n to S, forms an acute angle with 

the z-axis, we take the plus sign in the right-hand member. 

Thus, n = —i + k and, hence, (a, n) = z — y. 
Projecting S, onto the zy-plane, we get the semicircle 


Day: << R?— уг. 
Then by (6) we have 


T= | f (а, m] ded, 


Dry 
and again passing to polar coordinates, we find 


2/2 R 
П, = | (cos p— sin Фф) dp f p? dp = 2 R. 
-л/2 0 


(3) On $у, where 22 + y? = R?, that is, on the lateral 
surface of the cylinder, we have 


o_ втаб (212 -у2— R?) — ау _ и) 
п тшнк) = узи Ж ОН oC 
In this case, the normal п forms a right angle with the 
z-axis and therefore cos ү = 0 and so the choice of sign 
in the right-hand member is arbitrary. Take the plus 

sign and then 


n= 1+0, (a, n9) = ery 


and so 


It is impossible to project the surface S, (right cylinder) 
onto the £y-plane since it projects into a line, a semicircle 

“(the one-to-one nature of the projection will be upset). 
The same occurs when projecting onto the rz-plane. We 
therefore project the surface S, onto the yz-plane, onto 
which it is projected one-to-one into the region D,, 
bounded by the line 


a+ у= Re, | 
Z=T.- 


Eliminating z from this system, we obtain the equation 
for the projection of this line onto the zy-plane: z? + y? = 
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= А? is a circle. Since 


zz; (20. 


РРР i) | = | (n5, i) | = | R 


we will have 


п) Б. сват аш 
y: yz 
= {| 204-2) | y dy dz. 
Dy: yz 


Using the polar coordinates у = p cos ф and z =p sin ф, 
we get 


R 
1-2 2 d dp = 0. 
3 SJ pcos pp dp dp = f coso ofe p 


Thus, 
2 2 3 
П = — л+л +0=0. 


108. Compute the flux of the vector field a = yi + 2j + 
+ zk through the upper side of a triangle bounded by the 
planes 


ttytz=a, т=0, у= 0, 2 = 0. 


109. Compute the flux of the vector field а = zzi through 
the outer side of the paraboloid 2 = 1 — z? — y? bounded 
by the plane 2 = 0 (2 > 0). 

110. Compute the flux 7 the vector field а = zi + zk 
through the lateral surface of the circular cylinder y = 
= V R? — z? bounded by the planes z = 0, z = h (h > 0)» 
111. Compute the flux of the vector field a = zi + yj + 
+ zk through the upper side of a circle cut out of the 
plane z = k (k > 0) by the cone z = Из + y*. 

112. Compute the flux of the vector field a = 3zi — yj — 
— zk through the outer side of the paraboloid z? + y* = 
= 9 — z located in the first octant. 

113. Compute the flux of the vector field 


а = (1° + у)і + (2 + 28) ј + (2 22) К 
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through the portion of the plane z = 0 bounded by the 
circle z? + y? = 1 in the direction of the unit vector k. 
114. Compute the flux of the vector field a = yzi — zj — 
— yk through the total surface of the cone z? + y? = 2? 
bounded by the plane 2 = 1 (0 5 2 « 1). І 
115. Compute the flux of the vector field а = 2zi + 
+ (4 — 2y) ј + 2zk through the closed surface bounded 
by the paraboloid z? + z = 1 — 2y (y > 0) and the 
plane z = 0 (z > 0). 
116. Compute the flux of the vector field a = z?i + 
+ у + 2k through the total surface of a pyramid 
bounded by the planes z+ у + 2 = 1, х=0, у = 0, 
2 = U. 
117. Compute the flux of the vector field a = zi + yj + 
+ zk through the sphere z? + у? + 22 = Rè. 

2°. The method of projection onto three coordinate planes. 
Suppose a surface $ is projected one-to-one onto all three 
coordinate planes. Denote by Dzy, D,,, Dy, the projec- 
tions of S onto the zy-, zz-, yz-planes respectively. 

In that case the equation F (z, y, z) = 0 of the surface 
S is uniquely solvable for each of the arguments z, y, z 
so that 


z = z (y, 2), y=y(z,z) z = z (z, y). 
Then the flux of the vector 
а = P (т, Y, z)i +Q (z, Un 2) j + R (z, У, 2) К 


through the surface 5, the unit vector of the normal of 
which is 


n? = cos a.i + cos B-j + cos y-k, 
can be written thus: 


П = | (a, n°) dS 


= | | [P (x, у, 2) cosa +Q (z, у, 2) соз В 
5 


+ R (z, у, 2) созт] 45. (8) 
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We know that 
dS cos а = +dy dz, 


dS cos В = + ах dz, (9) 
dS cos y = +dr dy, 
the sign in each of the formulas of (9) being chosen to 


coincide with the sign of cosa, cos B, and cos y on the 
surface 5. Substituting (9) into (8), we get 


П = + | [Plz (у, 2), у, z] dy dz + | Jets. y (z, 2), 2] dz dz 


Dy: Dxz 


+ № R{z, у, z (z, y))dzdy. (10) 
Dry 
Example 10. Find the flux of the vector 
a = zyi + yzj + zzk 


through the portion of the external side of the sphere 
z? + y? + 22 = 1 located in the first octant. 
Solution. We have 


= [gud (рина) = урар АНИ, 
whence, taking into accqunt that the surface 5 lies т 
the first octant, we obtain 


соз а =z > 0, соз В = у >> 0, созү = 2 > 0. 


We therefore take the plus sign in (10) in front of the 
integrals, and putting 


P = ту, Q = yz, В = zz, 


we obtain 
П = | | zy dy a) | иби | | 12атау. (41) 
Dyz xz xy 


From the equation of the sphere z? + y* + z? = 1 we get 
z=z(z, y)-V1—z—y, y=y(z, 2) = И 1—22 — 22, 
2=2 (у, 2) = 1—22. 


Substituting these expressions for z, у, 2 respectively 
into the third, second, and first integrals on the right of 
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(11), we get 


П = 1—z2— y?dzd. 1—32— zi drd 
je" 12—у rat | psy 12 — 22 dz dz 


+ M yVi—y-—zdydz. (19) 
Dy: 


Let us compute the first integral on the right and pass to 
the polar coordinates z = р соз ф and у = р siny, 
where 0 < 9 < x/2, 0 x: p « 1. This yields 


І, = | | ху 1— 12— узагду = | f p? V 1— р? cos p de dp 
Dry Dry 
n/2 1 1 
- | cos e de | eV 1—pidp-— f p? Y 1— р? dp. 
0 0 0 
Setting р = зіп t, dp = cos t dt in the last integral, we 
have 
л/2 n/2 
1, = | sin? t cost dt = + | sin?2tdt = те. 
The second and third integrals in (12) аге computed in 
similar fashion and we obtain 


I= | | 2 И1—ж#—4х4:=ць, 


Dxz 


һ- {зуи 
yz 
The desired flux is 


П=1+1,+1,= 3. 


118. Use the method of projecting onto all three coordinate 
planes to compute the flux of the vector field through a 
surface S. 
(a) a = zi — zj + yk; 
S is the upper side of a bounded portion of the plane 
3х + бу — 22 = 6 cut out by the coordinate planes. 
(b a=(z +y +2)i + (@т+у+—1)}— 2; 
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5 is the upper side of part of the planez + у + 2 — 1 = 0 
lying in the first octant. 


(с) a=(1— Уу 22) i+ 5+ (Vy —22—2)k; 


S is the outer side of the paraboloid of revolution y = 
= 22 + 2? bounded by the plane у = 4 and lying in the 
first octant. 

3°. The method of introducing curvilinear coordinates 
on a surface. In certain cases, when calculating the flux 
of a vector field through a given surface S it is possible 
to choose a simple coordinate system on the surface 
itself to compute the flux instead of projecting onto 
coordinate planes. 

Let us consider some special cases. 

Case (1). Suppose a surface S is part of the circular 
cylinder z? + y? = R* bounded by the surfaces z = 
= fı (z, y) and 2 = f, (т, y), and we have f, (z, y) < 
< fa (z, y). 

Setting 

z = Rcosọ, y=Rsing, 2 = 2, 


we have for the given surface 
0< p< 2л, }, (Reosg, Rsing) S 2 Ф 
< f, (R cos ф, R sin 9), 
and for the element of area dS we obtain the following 
expression (Fig. 24): 
dS — R dq dz. 
Then the flux of the vector field a through the outer 
side of the surface S is computed from the formula 
2n /«(R cos 9, R sin 9) 
П= А \ dọ (a, по) dz, (13) 
0 f (R cos Ф, R sing) 
where 
o... grad (z3-Fy?— R?) _ zi-ryj 
— [grad (z?+y8—R%)| А c 
Example 11. Find the flux of the vector 
r = zi + yj + zk 
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through the outer side of the lateral surface of the cir- 
cular cylinder z? + y? = R? bounded by the planes 
z = 0 and z = H (Н > 0). 
Solution. Here we have 
0< p < 2л; f,(Reosg, R sin ф) = 0 
fa (R cos o, R sin 9) = H 


Fig. 24 


Introducing curvilinear coordinates on the cylinder, 
we get 
z= R соз ф, у = Rsin g, 2 = 2, 

Then the desired flux of the vector г is 

2л H 

П=В j dg | (г, n) dz. 

0 0 

But since 
r = zi + yj + zk = Reosqg-i+ А зір ф:ј + zk 

and the normal n° on the cylinder is 


p — ЖИ. — Ass ght Reet u cosg. -i+sing: j 


it follows that the scalar product on the cylinder will be 
(г, n°) = R cos? ф + R sin? ф = В. 
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Finally, we obtain 
2n H 


П= Аз | do | dz= 22H. 
0 0 
Example 12. Compute the flux of the radius vector 
r = zi + yj + zk 

through the lateral surface of the circular cylinder 
1? + y? = 1 bounded from below by the plane z + y + 
+z = 1 and from above by 
the plane = +y + 2 = 2. 


Solution. Here (Fig. 25) 
we have 


R = 1,f,(z, у)=1—— 
— у,ў. (у) = 2 — = y. 


Passing to coordinates оп 
the cylinder 
т = созф, y=sing, 


2 = 2, 
we get 
fı (т, y) = 1 — созф — зіп Ф, 
В (z, у) = 2 — соз o — 
— Sin Ф. 


According to (13), the flux 
of the vector r is 


2n 2-cos ф-віл Ф 
П = | de | (г, n°) dz. 


ù 1 -cos 9-3in Ф 
But since on the cylinder z? + y? = 1 we have 
n? = zi + yj = cos g-i + sin Ф.ј, 
it follows that 
(г, n) = 22 + у? = cos? pt sin? ф = 1 
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and, hence, 


2л 2—cos ф—51п 9 2л 
T= f dọ dz = | йр=2л. 
0 i-cos ф—81п 9 0 


119. Find the flux of the vector 
a = yi + zj — ek 


through the outer side of the lateral surface of the cylinder 
a? + y? = 4 bounded by the planes z = 0 and z -+ y + 
+ 2 = 4. 
120. Find the flux of the vector 

a = zi — zyj + zk 
through the outer side of the cylindrical surface z? + 2? = 


= А? bounded by the planes y = 1 and z + y = 4. 
121. Find the flux of the vector 


а = 231 — уЗ} + zzik 


through the outer side of the cylindrical surface z? + y? = 
= 9 bounded by the sphere x? + у? + 22 = 25 
122. Find the flux of the vector field 


a = zi — yj — хук 


through the outer side of the lateral surface of the cylin- 
der 22 + y? = 1 bounded by the plane 2 = 0 and the 
hyperbolic paraboloid 2 = z? — y’. 

123. Find the flux of the vector field 


а = (zy — у) i + (22 — z? + zy) j + zk 


through the outer side of the lateral surface of the cylin- 
der z? + у? = 1 bounded by the elliptic сопе 22 = 
= 2/2 + у. 

Case (2). Suppose the surface S is a part of the sphere 
z? + у’ + 22 = А? bounded by conical surfaces whose 
equations in spherical coordinates have the form Ө = 
= 7, ($), 0 = fa (Ф) and by the half-planes Фф = Фф, 
ф = Po. 

For the points of the given sphere, set 


z= R cos фзіп Ө, у = Д ѕіп фѕіп Ө, 2 = R соз Ө, 
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where 9, < 9 < Фф, and 0, < 0 < Ө,. Then for the 
element of area dS we obtain (Fig. 26) 

dS = Ң? sin Ө d0 dq. 


In this case, the flux of the vector field a through the 
outer part S of the sphere is found from the formula 


Ф, е, 
П = R? | de | (a, n?) sin 040, (44) 
Ф, 0, 


where 
o... grad(z?--y3-Fz?— R?) _ zi-ryjd-:k 
“| grad (z3++ y?+ 2?— А?) | — R d 


Example 13. Find the flux of the vector 


Fig. 26 


а = (z — 2y + 1) i + (21 + y — 32) i + By + 2) К 


through a part of the surface of the sphere z? + y* + 2? = 
— 1 located in the first octant into a region where 
а oy* £z 1. 


6* 
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Solution. Here we have 
R=1, ф = 0, Ф, = 7/2, 0, = 0, 0, = n/2, 

n? = zi + yj + zk, (a, n?) = 22 + y+ 24+ =. 
On the sphere z? + у? + 22 = 1 we introduce coordinates 
ф and Ө so that 

т = cos ọ sin Ө, у = зіп фзшб, 2 = cos 0. 
We then have 
(а, n?) = 1 + cos y sin 0 
and, using (14), we obtain 


л/2 л/2 
П= ( dq | (1+ cos ọsin Ө) sin 020 
0 
п/2 n2 n/2 n/2 
- Top sas Posed] aoai 
= ф | sin + | созфаф | sin 040 — 7 n. 
0 0 0 0 
124. Find the flux of the vector field 
а = Fi — у] + zk 
through the outer side of that part of the sphere z? + y? + 
+ 2 = 1 cut out by the conical surface z? = 22 + у, 
z> Va + И. 
125. Find the flux of the vector field 
а = ya + zzj + zyk 


through the outer side of that part of the sphere z* + y* + 
+ 22 = А? located in the first octant. 
126. Find the flux of the radius vector 
= zi + yj + zk 

through the outer side of that part of the sphere z? + y? + 
+ 2? = 2 that is bounded by the planes z = 0 and z = y. 
127. Find the flux of the vector 

а = zzi + yzj + 2k 
through the outer part of the sphere z? + y? + 22 = 9 
cut off by the plane z = 2 (z > 2). 
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Sec. 12. The flux of a vector 
through a closed surface. 
The Gauss-Ostrogradsky theorem 


Theorem. If in some region С of space the coordinates 
of a vector 


а = P (z, у, 2) i + Q (z, у, 2) } + R (z, у, 2) К 


are continuous and have continuous partial derivatives 
aðPlðz, д0!ду, AR/dz, then the flux of the vector a through 
any closed piecewise smooth surface = located in G is equal 
to the triple integral of OP/dx + 3Qlðy + AR/dz over 
the region V bounded by the surface 2: 


п= (а, n?) do = 11 (52450-32) ао а) 
r 


(the Gauss-Ostrogradsky formula). 

The normal n to the surface X is taken to be the outer 
normal. 

Example 1. Compute the flux of the vector 


а = ri + yj + 2k 
through the closed surface 
і? + у + 22 = А, 2 = 0 (2 Р 0). 
Solution. By formula (1), 
П = (2x + 2y + 22) dv. (2) 
И 


The integral (2) is conveniently computed т the spherical 
coordinates г, Ө, ф. We have 


z=rsinOcosg, у = гзіп Өзіп ф, 2 = гсоз 0 
and the element of volume is 


dv = т? sin Ө dr dð dq 
so that 


П=2 Na (г sin cos q + rsin8sin 9 
Y 


+r cos 0) г2 sin 0 dr dô dọ 
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2л л/2 R 
= 2| dọ | sin Ө (sin Ө cos ọ + sin Ө sin ọ + cos 0) аө | r3 ar 
0 0 0 


2л 
= 28 {а of соз Ө sin @ 40 = 7" | 
0 0 
Example 2. Compute the flux of the vector 
а = Áázi — yj + zk 


through the surface of a torus. 
Solution. Taking advantage of the Gauss-Ostrogradsky 
theorem we find that the desired flux is 


n= Poo по) do = e (92.50 2E) av — àv, 


where V is the volume of the torus. To compute the 
volume V, let us take advantage of the Guldin theorem 


Fig. 27 


on the volume of a solid of revolution, by virtue of which 
the volume is equal to the product of the area of the 
rotating figure into the path covered by the centre of 
mass of the figure during the rotation. 

Let A, and А, be the inner and outer radii of the torus 
(Fig. 27). The area S of the circle, which during rotation 
forms the torus, is equal to 


$— 5 (Sz Pe y. 
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The path length described by the centre of mass (the 
centre of that circle) is the circumference l of a circle 
of radius (R, + R,)/2, that is, 


l= 2a CET s (n, ER). 


Thus, the volume V of the torus is equal to 


= 2 
Vos (He) s e Ry) = 97 (Н,— Ry: (Rot RJ. 
The desired flux is 
=? (Ry — Ry)? (Bs Y Ry). 


Example 3. Using the Gauss-Ostrogradsky theorem, 
compute the flux of the vector field 


. 2r 1T y? 
a= (455 Ji i+ 2z arctan y-j— CFP IES k 
through the outer side of that part of the surface z = 
= 1 — z? — y? located above the zy-plane. 
Solution. In order to be able to apply the Gauss- 
Ostrogradsky theorem, close the given surface from 


below with a portion of the zy-plane that is bounded by 
the circle 


gi у? =1, | 
z=0. 


Let V be the volume of the resulting solid bounded by 
a closed piecewise smooth surface o consisting of a part 
01 of the paraboloid of revolution z = 1 — z? — у? and 
a part о, of the plane z = 0 (Fig. 28). 

The flux of the given vector through the surface o is, 
by the Gauss-Ostrogradsky theorem, equal to 


= Qe "ше (=> T. +2) av. 


We find a sum 


Р, 20 , ðR 2z (1+) 
++ р Si: 


И Чу три 
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Consequently, the flux is 
n-$ (a, n?) do — 0. 
o 


By virtue of the additivity of the flux we have 
П = M (a, n?) do+ f | (a, n?) do = 0. 


Fig. 28 


From this the desired flux is 


П, = (a, по) do = — (a, n°) do. 

J) j! 
The flux П, of the vector a through the circle z? + y? < 
<1, 2 = 0 is equal to 


П, = (a, n°) do. 
j 
Since on the plane z = 0, we have 
даз 


ат Ра 


i+ 22 агсіап у-ј — К, по = — К, 
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and hence (a, n°) = 1, it follows that the flux IT, through 
the circle с, is equal to the area of the circle до: 


I, = | }4о=л. 


о, 


The desired flux П, = —П, = —л. 

By appropriately closing the given unclosed surfaces 
and making use of the Gauss-Ostrogradsky theorem, com- 
pute the fluxes of the vector fields through the indicated 
surfaces (the normal to the closed surface is the outer 
norma)). 

128. а = (1 — 2z) i + yj + zk; 5: z*- yp = 22 

(0 x z « 4). 

129. a = zi 24 + yk; 5: 22 + у? = 4 — 2 (229). 
130. а = (у? + 22) і — yj + 2yzk; S: -+ 22 = у? 
(0 x y « 1). 


Sec. 13. The divergence of a vector field. 
Solenoidal fields 


The notion of the flux of a vector through a closed 
surface leads to the concept of the divergence of a field. 
This concept yields a certain quantitative characteristic 
of a field at each point in the field. 

Let M be a point of the field under study. Surround it 
by a surface 2 of arbitrary shape, for instance, a sphere 
of sufficiently small radius. Let the region bounded by 
the surface = be (V) and its volume V. We consider 


the ratio 
$ (8, n?) da 
I 


——. (1) 


Definition. If the ratio (1) has a finite limit when the 
region (V) contracts to the point M, then this limit is 
termed the divergence of the vector field (the divergence 
of the vector a) at the point M and is designated as 
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M (a, n9) do 


diva(M) = lim —,——. (2) 
(V)- M 


div a (M). We thus have 


Formula (2) yields an invariant definition of diver- 
gence. This definition means that the divergence of the 
field a at the point M constitutes the volume density of 
the flux of the vector a at that point. 

The points M of the vector field a (M) at which div 
а > 0 are termed sources, while the points at which 
div a « 0 are termed sinks of the vector field. 

The divergence of a vector field is a scalar function of 
the points of the field. 

If the coordinates of the vector 


а (M) = P (=, у,2)і + Q (т, у, 2) і + R (z,y,2)k 


have continuous partial derivatives дР/дх, dQ/dy, AR/dz 
in the neighbourhood of the point M (z, y, z), then, using 
the invariant definition of divergence, we find from the 
Gauss-Ostrogradsky theorem that 


; P 
diva- AXE. (3) 


All the quantities in (3) are considered at the same point 
M (z, y, 2). 

Using (3) for divergence, we can write the Gauss- 
Ostrogradsky theorem (see Sec. 12) in vector form: 


(ў (a, по) do = У div adv (4) 


> 


Example 1. Using the invariant definition, compute 
the divergence of the vector a = zi at the point O (0, 0, 0) 
taking for the surfaces o surrounding the point O the 
spheres o, of radius e centred at that point. 
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Solution. By the definition of divergence, we have at 


the given point 
e (a, 00) do 


diva(0)= lim = 


(0,)-+0 ve ' 


where v, is the volume of a ball bounded by a sphere 


Oe, ОГ 
q (a, n9) do 


div a (0) = lim —* 
e-0 


But since the volume of the ball is v, = 4ne?/3, it fol- 


lows that 
А) (a, n^) do 


div a (0) = lim —* 
£0 nes 


w| 


Let us compute the flux faa, n°) do of the given vector 


О, . 
through the sphere де. The unit vector of the normal n? 
to the sphere о, is directed along a radius of the sphere, 
and so we can put 
A r 
EET 


where r^ is the unit vector of the radius vector r — 
= zi + yj + zk, or 


по = го 


— rI 
ыг. 


zi+yj+zk 


e 


n’ = 
The desired flux is 


М) (а, по) do = А) ra do. 


о, 
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Passing to coordinates on the sphere o,, 
xz = e cos ọ sin Ө, у = езіп фзш Ө, 2 = е соз Ө, 
we get 


3 1 ina 3o 
@ в, 19) dom ( [ steatgsinteeruneceas 
о, с 


е 
2л л 4 
= 23 f cos? ф dọ f sin’ ọ dp-—4 лез. 
0 0 
Consequently 


A cs 
з ЛЕ 
div a (0) = lim – 


£40 4 дез 
3 


=1. 


Example 2. Compute div г. 
Solution. We have г = zi + yj + zk, so that P = т, 
Q = y, R = 2 and, hence, by formula (3), 


берс ФЕ у дуу И 
divr= a t Фу + = = 3. 
Example 3. Compute div (и. a), where и (М) is a scalar 
function and a (M) = Р (z, y, z)it+Q(z, у, 2j + 
+ В (z, y, 2) К is a vector function. 
Solution. Using formula (3), we get 
div (ua) 
..9(uP) , 8(uQ) , 9(uR) , ӘР ди д0 ди 
ee im Е N 
ôP 6Q ôR ди 
- utr vb 
ди ди 


+09 + RA = ші а + (а, grad и). 


OR ди 
Tuy Ta a 


Thus 
div (ua) = u div a + (a, grad и). (5) 


Example 4. Find the divergence of the vector 


а-=Ф (г) r° = 70) г, 
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where r = |r | is the distance from the coordinate origin 
to the variable point M (z, y, 2). 
Solution. Using formula (5), we obtain 


diva = divr+ (r, grad 20 ) | 
Furthermore, 
divr=3, grad at) — ( 90 ) eee (900 " 
And so 
div а= 20.3, (X =e) го, г) 


=3 20) je 0009 e) +9’ (r). 


131. For what function %р (г) will we have div ẹ (г) г = 
= 24 (г)? 
132. Find div (rr). 
133. Find the divergence of the vector field 
а = (е, г], 

where € is a constant vector. 
134. Find 

div (г [w, rl), 


where w is a constant vector. 
135. For what function :p (z) will the divergence of the 
field 
а = zzi + yj + v (2) К 
be equal to z? 
136. Find the flux of the radius vector r through the 
surface of a sphere. 
137. The electrostatic field of a point charge q is 


= 2 
Compute div E. 
138. Show that 


+6 (г, по) do = V, 
E 
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where V is the volume bounded by a closed surface X. 
139. Prove that if У is a closed piecewise smooth surface 
and e is a nonzero constant vector, then 


§ cos (u, е) da =0, 

z 
where n is a vector normal to the surface 2. 
140. Prove the formula 


9 (pa, n?) da = | | | (p div a+ (а, grad ф)) dv, 


where ф = ф (z, у, 2) and È is a surface bounding the 
volume V. Establish the conditions under which this 
formula is applicable. 
141. Prove that if the function и (т, у, 2) satisfies the 
Laplace equation 

д%ц Ou ðu 

эз 3395 


e 25 do — 0, 
> 


where ди/дп is the derivative with respect to the direction 
of the outer normal to the piecewise smooth closed sur- 
face 2. 

142. Prove that if the function м (z, у, 2) is a second- 
degree polynomial and 2 is a piecewise smooth closed 
surface, then the integral 


then 


is proportional to the volume bounded by the surface 2. 

Find the flux of a vector field through the indicated 
closed surfaces: (1) directly, (2) via the Gauss-Ostrogradsky 
theorem. 


= х®+ y?, 


143. а= ті + zk; 5: | ud 
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144. 


145. 


146. 


147. 


148. 


149. 


150. 
151. 


152. 


. a=zi—2yj—zk; 5: { 


29 = 22+ yt, 
z—H, Н > 0. 
2=6— 12 — y?, 
22 = 12 02, 2220. 
z+ 22 = у, 
22 = z+ у, 
2 =1 +02. 

gi yt 4, 
а= 2гі— (2—1) К; 5: ee 2= 4. 
+ y+ 22 = 4, 
32 = 12+ у? (i2 T1). 


a—2zi4-2yj—zk; S: { 
a= zi— zj; s: | 
a—yzi—zj—yk; 5: { 


a=zi+2yj—zk; 5: { 


a=2zri—yj+ zk; 5: { 


a-—yzi--2yj—zk; S: z?4-y?--22— 4. 

а= (2+ 2) i-F (y - 2) }+ (z+ 2) Е; 

8: didi: 

* [z=y, 2;>0. 
9—2= 12 y, 

а = 3:1 — уј— zk; S: 4 == 0, у= 0, z=0 (first 
octant). 


r+y+z=1, 


. a=(y—z)i+(z—y)j+(z—z)k; 5: [set 


х=0, z=0. 
1—2= 22402, 
2= 0. 


Solenoidal fields 


Definition. If at all points M of a certain region G 
the divergence of a vector field (specified in G) is zero, 


div a (M) — 0, 


then we say that the field is solenoidal in that region. 
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Thus, a solenoidal field is, by definition, without 
sources and sinks. 
From the Gauss-Ostrogradsky theorem it follows that 
іп a solenoidal field the flux of a vector а = a (M) 
through any closed surface c 
n lying in the field is zero: 


(a, по) do = 0. 
3 


In a solenoidal field G, the 
vector lines cannot begin or 
end. They can either be closed 


L curves or have ends on the 
boundary of the field. 
Fig. 29 The equation 
div a (M) 20 


is encountered in hydrodynamics and is termed the con- 
tinuity equation of an incompressible fluid. 

In this case, the amount of fluid exiting through some 
closed surface o is always 'equal to the amount entering, 
and the total flux is zero. 

Which of the following vector fields is a solenoidal 
field? 

155. а = z (z? — 0?) і + y (22 — 2°) j + 2 (у? — z?) k. 
156. a = y — (zt + y) j + 2 (3у° + 1) К. 

157. а = (1 + 22у) i — ия + (zy — 22у + 1) К. 
158. Show that the field of the vector 


Е= 21° (r-Vzcryrz) 


is solenoidal throughout the region that does not contain 
the coordinate origin О (0, 0, 

159. Under what condition will the vector field a = 
= Ф (ғ) г be solenoidal? 

Suppose we have a field a (M), which is not necessarily 
solenoidal. In the field, consider a closed oriented con- 
tour L. The surface È containing the curve L as its edge 
will be called a surface spanned by the contour L. Let us 
agree to orient the normal n to the surface Х so that the 
chosen circuit about the contour L will be seen from the 
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end of the normal as being counterclockwise (Fig. 29). 
160. Show that in a solenoidal field the flux of a vector 
a (M) does not depend on the type of surface 2 spanned 
by the given contour L and depends solely on the contour 
itself. 


Sec. 14. A line integral in a vector field. 
The circulation of a vector field 


Suppose we have a continuous vector field a — a (M) 
and a piecewise smooth curve L on which a positive 
direction has been chosen (in other words, L is an oriented 
curve). 

Definition 1. The line integral from the vector a = 
= а (М) along the oriented curve L is termed a line 
integral of the first kind (an integral over the arc length 
of a curve) of a scalar product (а, т): 


j (a, 1^) ds, 

L 
where 1? = x? (M) is the unit vector of the vector tangent 
to the curve L whose orientation coincides with that 
of L; ds is the differential of arc length s of L. 

If r — r (M) is the radius vector of an arbitrary point M 

of the curve L, then the line integral in the field a (M) 
may be written thus: 


MOL (a, dr), (1) 
L L 


If a rectangular- coordinate zyz-system is introduced 
in the vector field, then г = zi + yj + zk, 
а = P (z, у, 2) 1 + Q (т, у, 2) і + В (т, у, 2) К 


and the line integral (1) can be expressed in terms of а 
line integral of the second kind: 


| (a, dr) = f Р (т, y, 2) dz + О (т, y, z) dy + R (т, у, 2) dz. 
L L 
When a — a (M) is a force field, the line integral (1) 


yields the magnitude of the work of that field over the 
curve L. 


7-910 
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Properties of a line integral 
(a) Linearity: 


[Oe pas dr) 55. | (ay, dr) + | (аз, dr), 
L L 


where À and p are constants. 
(b) Additivity: 
(a, dr) = Í (а, dr) + | (a, dr), 
[+13 Li Lo 
(c) The integral reverses sign with a change in the 
orientation of L: 


) (a, дг) = — f (a, dr), 


where A is the initial and B is the terminal point of the 
curve L. 


Calculating a line integral in a vector field 
Let the curve L be specified by the parametric equations 
z = 9 (t), у = № (t), z = x (t) (to tx t5). 
Here, at the initial point A of L the parameter ¢ assumes 
the value ¢ = і, and at the terminal point B of L it 
assumes the value ? = 1, (the direction on L corresponds 
to increasing values of the parameter t from tọ to t); 


the functions ф (t), sp (t), x (t) have continuous deriva- 
tives on the interval [fy, ¢,). Then 


ty 
| а.а) = | (а,й) = | (Pie, 0, х(@1Ф' (0 
L AB to 
+Q [o (t), p (£), x (0) $ (+ Rig (2), p (£), хх (0) dt. 
If the curve L is given by a system of equations 
y = (2), z = x (2), а < z < b, then 
b 
(a, dr) = | (Рі, (2), (21+ Q Ens $ (2), х (2) 9 (2) 
AB a 


+ Riz, ф (2), x (2)) x (z)} dz. 


Sec. 14] A line integral in a vector field 99 


Similar formulas may also be written for cases where 
the line is specified by one of the following systems 
of equations: 


х=Ф(у), z—-x(y (о <у< И) 


z= Фф (2), у= № (2) (2 <2< 4). 


Example 1. Find the line integral of the vector а = 
— r/ | r |, where г is a radius vector on the line segment 
from point A (r4) to point B (гв). 

Solution. The desired line integral is 


ја, а) = D (1) 
AB AB 


or 


From 
d (r, т) = (dr, т) + (т, dr) = 2 (r, dr) 
we find 
(г, dr) - 4(т,т) = d (Ir) = 3-2 |r] d |e] = Ir] d [e], 
whence 


(г, дг) — 
A — d |r]. (2) 


Substituting (2) into the integral (1), we get 


Tg 


| ка, ао) = fairl | an = еа — rna. 
АВ АВ ГА 
Note that | 
| dr | 52а jr |. 


Find the line integral over the line segment bounded 
by the points A (rj) and B (r,) for the following vector 
fields: 

161. a — r. 


162. a = = 
r 


163. a= TIL ‚ r is the unit vector, 


164. Compute the line integral over the straight line 
passing through the points O (0, 0, 0) and M, (1, 1, 1) 


7* 
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from point О to point M, if a = [b, r) , where b is a con- 
stant vector. 
165. Prove that 


| (grad и, dr) =u (B) —u(A). 
АВ 


Example 2. Find the line integral of the vector 
a = zi + zj + yk 
over an arc L of the heli- 
cal curve 
r-—Hcost, y=Rsint, 
t 
=з 
from point A, the point of in- 
tersection of the curve with 
the plane z = 0, to point B, 
the point of intersection with 
the plane z = 1 (Fig. 30). 
Solution. Here, the line in- 
tegral is of the form 


| (a, dr) 
L 
= f zdz + zdy + y dz. 
L 


The helical curve is located 
Fig. 30 on a circular cylinder z? + 
+ у? = R. At point А we 

have t = 0, at point В we have t, = 2л. Since 


4х= — Rsintdt, dy = Rcostdt, di- 7, 


it follows that the integral is equal to 
2л 
f (a, dr) = | ( — R sin t + А? cos? + + sin t) dt 
1 0 
2n 2n 
= А? | cos? t dt — .É- | tsintdt=nR2+R 
0 


2л 
0 


Sec. 14] A line integral in a vector field 104 


because 
2л 
| cos? t dt = л; 
0 
2n 
| tsintdt= —2n. 
0 


Example 3. Find theline integral of the vector (see 


example 2) Е 
а = 2 у 


over the straight line AB (see Fig. 30) in the direction 
from point A to point B. 

Solution. Since the straight line AB (the generatrix 
of the cylinder z? + y? = R?) is located in the zz-plane 
and passes through the point A (R, 0, 0), it follows that 
у= 0, z = R, dz = 0, and for the radius vector г of 
the points of AB we will haver = Ri + zk, dr = dz-k. 
Therefore the scalar product 


(a, dr) = z йт + z dy + y dz 


on AB will be zero. 
Hence, the desired line integral 


| (a, dr) = } (a, dr) 


on AB will be zero. 

From examples 2 and 3 it follows that in the general 
case a line integral in a vector field depends not only on 
the initial and terminal points of the path of integration 
but also on the shape of the path. 

Example 4. Compute the work done by the force field 

F = yi + rj + (z+y+z2)k 
along the segment AB of a straight line passing through 
the points М, (2, 3, 4) and M,(3, 4, 5). 

Solution. The work done by the force field will be 

equal to the line integral along М,М,: 


A= йе f ydzJ4-zdy--(rJ-y + 2) dz, 
М1 Мз 


MiM: 
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Let us find the canonical equations of the straight line 
M,M;,. We have 


4 o4 op 
whence 
y=2+1, 
m 


Here, т varies from 2 to 3 (since the abscissa of М, is 
equal to 2 and the abscissa of M, is 3). The desired work is 
3 


A- Горен toe 


3 
= | (5z4-4) dr 5. 

2 
166. In the plane vector field 
у — 23] 
Vite 
compute the line integral along the semicircle 

z= Rcost, y= Rsint (0x t x n). 
. In the plane vector field 
а = (1% + y)i+@—y)j 
compute the line integral over the curve у = | = | from 
the point (—1, 1) to the point (2, 2). 
168. In the plane vector field 
а = (2? — 2zy) i + (у? — 2zy) j 

compute the line integral: 

(a) along the parabola у = т? from the point (—1, 1) 
to the point (1, 1); 

(b) along a segment of the straight line joining the points 
(—1, 1) and (1, 1). 
169. Compute the work of the force field F = 2zyi + z*j 


along an arc of the circle z? + у? = 1 from the point 
(1, 0) to the point (0, 1) counterclockwise. 
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170. Compute the work of the force field 

F = (2 + 2zy) i + +} 
along the parabola у = 2? from the point (0, 0) to the 
point (4, 1). 
171. In the vector field 

"M zi-Fyj-4- zk 

compute the line integral along a segment of the straight 
line from the point (1, 1, 1) to the point (4, 4, 4). 
472. In the vector field 

а = (у? — 22) i + 2yzj — z'k 
compute the line integral over the line 


r=t 
inet (0<t<1) 
2 = 13 


in the direction of increasing values of the parameter t. 
173. In the vector field 

a = yl + zj + zk 
compute the line integral along a turn of the helical curve 


z=acost, 
y —a sint, (0<t<2n) 
z—bt 
in the direction of increasing values of the parameter /. 
174. In the vector field 

а = ri + yj + zik 
compute the line integral in the direction from the point 
(0, 0, 0) to the point (1, 1, 1) along the line segment 
between the two points. 


Computing the circulation of a vector field 


Definition 2. The circulation C of a vector field a = 
= а (M) is a line integral taken around a closed oriented 
curve L Thus, by definition, 


C=§ (а, dr), 
L 
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where the symbol $ denotes the integral around the 


> 


closed curve L. 
If the vector field a = а (M) is given in coordinate 
form, 


a = P (z, Un 2)1 + Q(z, y z) j + R (z, У, 2) К, 
then the circulation of the vector field is 


C=§ Pdz+Qdy+ Raz. 
L 


For the positive direction of traverse around the closed 
curve L we take the direction in which the region bounded 
by the curve is on the left. 

Example 5. Compute the circulation of the vector 


field 
а = —yli + 23] 


around the ellipse L: 22/а? + у?/ = 1. 
Solution. By the definition of circulation we have 


C =$ (a, dr) = 6 — dz + z dy. (3) 
L L 
The parametric equations of this ellipse are of the form 
z=acost, 
у= | (0<t<2a), (4) 
whence 
dz = —a sin t dt, dy = b cos t dt. (5) 


Substituting (4) and (5) into (3), we get 
2л 
С =ab | (62 sin‘ / а? cos* t)dt= X лар (a? + 5?) 
0 


since 
2n 2л 


j entrare j (1— cos 2t)? dt 
D 


-+f (1—2 cos 21 + 159284) at 
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T [3 1 063 3 
f а: тше л. 
Similarly, we find that 

2л 

| cost tdt = > л. 

0 

Example 6. Calculate the circulation of the vector 
field 
a = уе" + хе?) + хуг 


around the curve L obtained by cutting the cone z? + y* = 
= (z — 1 with the coordinate planes (Fig. 31). 


Fig. 31. 


Solution. The curve L consists of two parts BC and CA 

located in the coordinate yz- and zz-planes respectively, 
~ 2? + у = 1 

and the arc AB of the circle 2 -0f[^ 


the circulation of the given vector field will be 


Therefore 


C= $ (a, dr) = | а.а) + | (a, dr) + | (a, dr). 
L BC CA АВ 
(1) Оп the line segment BC we have 


z = 0, dtr = 0; z = 1 — y, dz = —4у; 1 р у р 0. 
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Consequently, 
| (а, 4г) = | y dz = 0). 
BC BC 
(2) On the line segment CA we have 
у = 0, dy = 0; z =1 — z, dz = —dz; 05251 
and so 


f (a, dr) = | zdy =0. 


CA CA 
~ 12 + у? =1, 
(3) On the arc AB of the circle _0 | we have 
z=0, dzz 0, which means that 
| (а, 4г) = j e™ (y dz + z dy) = | era (zy) 
АЙ AB АВ 
E "ETT 
= {ае ue „=1—1=0. 


AB 
The desired circulation of the vector field is zero. 
Example 7. Compute the circulation of the vector 
field 
а = ryi + уг) + zzk 
if 
L: | т? + у? = 1, 


t+y+z=1. 
Solution. We have 


C=O (a, dr) = f zy dz + yz dy + zz dz. 
L L 


The curve L is an ellipse obtained by cutting the cylinder 
z*--Fy*-—1 by the plane z + y + z = 1. Let us find 
the parametric equations of the curve. The projection 
of any point of the curve on the zy-plane lies on the 
circle z? + y? — 1. From this we obtain z = cost, 
у = sin t. But the ellipse lies in the planez + y + 2 = 1, 
whence 2 = 1 — z — y or z = 1 — cost — sin t. Thus 


Sec. 14) A line integral in a vector field 107 


the parametric equations of the curve L are: 


r-cost, 
y — sint, (0x t «— 2л). 
2—1-—cost— sint 


From this we find 
dr = —sin t dt, dy = cost dt, dz = (sin t — cost)dt 
and so the circulation is 
2л 
С= | [— cos і sin? £ + sin # (1 — cost — зіп /) cost 
0 
+ cos t (1 — cos t — sin t) (sint — cos ¢)] dt 


2л 
= f (—3 sin? z cost + sin 21— cos? ¢ sin t—cos* ( + cos? t) dt 
0 


2л 
= - | cos? ( dt = — д. 
0 


Compute the circulation C of a vector a around the 
given curve L: 
175. a = (zz + y) i + (yz — х) j — (z* + y?) k; 
т?-+ y? =1, 
L: { ty 


2=3. 
176. а = yi + 22} + 27k; 


L z+ у2-{ 2? = R, 
4 22 +02 = Нг 


177. а = (2z + z) і + (2y — 2) } + zyzk. L із the line of 
intersection of the paraboloid of revolution z? + у? = 
= 1 — z with the coordinate planes. 

178. Show that if in a vector field the circulation of a vec- 
tor around any closed circuit is equal to zero, then there 
can be no closed vector lines in such a field. 


(z>0). 
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Sec. 15. The curl (rotation) 
of a vector field 


Suppose we have the vector field 

a (M) = P (z, y, 2) 1+0 (5, y, 2) j + В (z, y, 2) К. 
We will assume that the coordinates P, Q, R of the vector 
a (M) are continuous and have continuous partial deriv- 
atives of the first order with respect to all its arguments. 

Definition 1. The curl (or rotation) of a vector a (M) 
is a vector (symbolized: curl a (M) or rot a (M)) defined 
by the equation 


ƏR \. 

curl а= (5 —) i+ (4-4) 

90 _ 

+ (52 or — r2 k 

or, in easy-to-remember symbolic form, 
ij k 
д 0 д 

curl a= a ду 95 |: (2) 

P Q R 


This determinant is ordinarily expanded in terms of 
elements of the first row, in which case the operations of 
multiplication of elements of the second row by elements 
of the third row are regarded as operations of differen- 
tiation; for example, 


Definition 2. If in some region G we have curl a = 0, 
then the field of the vector a in G is said to be irrotational. 
Example 1. Find the curl of the vector 


а = (z + 2)i + (y + 2)j + (22 + z) К. 
Solution. Using formula (2), we obtain 
i j k 
д 9 9 
дт әу д: 
r+zy+2 22+2 


curl a= 
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Expanding the determinant in terms of elements of the 
first row and .regarding the operation of multiplication, 
say д/ду by x? + 2, as an operation of partial differen- 
tiation, we obtain 

curl a — —i— (2z — 1) j. 


Example 2. Find the curl of H, the intensity vector of 
a magnetic field under the conditions of example 3 of 
Sec. 10. 

Solution. The magnetic-field intensity vector H is 


2 
H- ilr 
or 
i jk 
2 2I... 2b. 
H=- 001 = И oF 9, 
ху 
where р? = z? + у?, whence, by (2), 
i j k 
д д д д 2Iz 
car Не) 3 x = (age) 
21у 21т 0 


eye GA) 

д 21у T34 у2— 21? x24 уз— 2? B 

+a (apa) k= 2 LE e e [К—0 
(22 + у? # 0). 


Thus, curl Н = 0 everywhere except the z-axis, at the 
points of which the last formulas are meaningless (the 
denominator vanishes), that is, the field of the vector H 
is irrotational everywhere outside the points of the z-axis. 

Find the curl of the following vectors: 

179. a = ane лы и + 22) ] + (22 + 27) К. 
180. а = a + + 
181. a = E (—y'i + Mm 
182. Show that if the coordinates of the vector a (M) 
have continuous partial derivatives of the second order, 
then 

div cur] a = 0, 
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that is, thefield of the vector curl a (M) is a solenoidal 
field. 
183. Show that 

(a) curl (a + b) = curl a + curl b, 

(b) curl (Aa) = А curl a. 
where À is a numerical constant. 
184. Show that if и = u (M) is a scalar function and 
а = a (М) is a vector, then 


curl (ua) = u curl a + [grad u, al. 


185. Show that if a and b are constant vectors, and r 
is the radius vector of the point M (z, у, z), then 


curl (т, a) b = (a, b]. 
186. Show that 


curl (га) = + (г, al, 


where a із а constant vector and r= |r| = 
=V2+/+ 2. 
187. Show that curl (f (r) a) = Lo [r, al, where f (r) 


is an arbitrary differentiable function of its argument and 
a is a constant vector. 

188. Show that the vector field а = f (r) г is irrotational 
that is, curl a = 0. 

189. Show that 


div [a, b] = (b, curl a) — (a, curl b). 


190. Show that the curl of a field of linear velocities v of 
a rotating solid is a constant vector in the direction paral- 
lel to the axis of rotation, the modulus of which is equal 
to twice the angular velocity of rotation: curl у = 20. 
191. Determine the angular velocity о of rotation of 
a solid about a fixed axis passing through some point of 
the solid if its linear velocity is 


у = 2zi + y*j + zzk. 


192. Show that the field of the curl of the vector a (М) 
is free of sources and sinks. 
193. What must the function f (т, 2) be so that the curl 
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of the vector field 
а = yzi + f (z, 2) j + zyk 


is coincident with the vector k — i? 


Sec. 16. Stokes’ theorem 


Suppose the coordinates of a vector 
a (M) = Р (z, y, z3) i + Q(z, у, 2) j + R (z, y, z)k 


are continuous and have continuous partial derivatives. 

Theorem. The circulation of 
the vector a arounda closed con- 
tour L is equal to the fluz of the 
curl of the vector through any 
surface Ў spanning the con- 
tour L: 


$ (a, dr) = |] (curl a, n?) do. 
(1) 


It is assumed that the orien- 
tation of the normal n? to the 
surface X is matched with the 
orientation of the contour L Ф 
so that, when viewed from the Fig. 32 
end of the normal, thecontour 
is traversed in the chosen direction counterclockwise. 

Example 1. Compute the circulation of the vector 


a = yi + zj — zk 


2+ 
around the contour L: | r 
(2) via the Stokes theorem. 

Solution. (1) The contour L is a circle of radius R = 2 
lying in the plane z = 3 (Fig. 32). We choose the orienta- 
tion as shown in the drawing. The parametric equations of 


B (1) directly and 
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the curve L are 


z-2cost, 
y —2sint, (0« t « 2л), 
2= 3 


so that 
аг = —2sintdt, dy —2costdt, dz=0. 
For the circulation of the vector a we have 
2л 
C= | 12 5щтЕ(—2 sin t) + 4 cos? 12 cost — 3.0] dt = — 4л. 
0 
(2) To compute the circulation via the Stokes theorem, 
choose some surface У spanning the contour L. For È 
it is natural to take a circle having. L as its boundary. 


According to the chosen orientation of the contour the 
normal n? to the circle has to be taken equal to k: n? — k. 


Then 
ij k 
д д 
cur] а = p Эу E = (2z—1)k 
y 22 —z 


Therefore, by the Stokes theorem, 


C= 1] (curl a, n°) do = |] (2=—1) do 
2л 2 
2 


2 
= | ae | (2p cos p— 1) p dp = — 2n E = — 4n. 
0 

194. Show that the flux of the rotor through an open 
surface spanning a given contour does not depend on the 

shape of the surface. 
Find the circulations of the vectors around the indi- 
cated contours (1) directly and (2) via the Stokes theorem. 


24 y? — 4, 
195. a=zitajtyk; L: |" ш. 
2= 0. 
P+ y+ = 4, 


196. a—yi—zjd-zk; Г: { ура (2>0). 
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197. а = 2xzi — yj + zk around a contour formed by 
intersection of the plane z + y + 2z = 2 with the coor- 
dinate planes. 


198. а=ш— 2} - (24у) К; Г: кш 

Pry? +22 = 16, 

z=0,y=0,z=0. 

2= 0 + 2, 

2=9. 

224-02 = 9, 

3y + 42 5. 

у z?—1, 
І = 2. 


199. а = 221; L: | 
200. а = zy?i + 222) + 22yk; L: | 
201. а = y*i--z?j; L: { 


202. a=yi—azj+z2k; L: | 


203. Given a vector field of velocities у of the points 
of a solid rotating with a constant angular velocity 
about the z-axis. Calculate the circulation of the field 
around the circle 


r-—acost, 
L: $y —asint, (0<1< 2л), 
z=0 


directly and via the Stokes theorem. 

From the Stokes theorem we find that the projection 
of the vector curl a on any direction n is independent of 
the choice of the system of coordinates and is equal to 
the surface density of the circulation of the vector a 
around the contour of the area perpendicular to that 
direction: 


$ (a, dr) 


рга curl ala = (curla, n°)| = lim + 5 
м (Х)-+М 


(2) 


Неге, (2) is a plane area perpendicular to the vector n; 
S is the area of (У); L is the contour of the area and is 
oriented so that the traverse of the contour is counter- 
clockwise as seen from the end of.the vector п; the nota- 


8-910 
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tion (Z) — M means that the area (2) contracts to the 
point M at which we consider the vector curl a, and the 
direction of the normal n to that area is always the same. 

Example 2. Compute the density of the circulation of 
the vector a — yi around the circle 


xr=acost, 
L: $y-asint, (O<t<2n), 
2-0 
at the centre of the circle, M (0, 0, 0), in the positive 
direction of the z-axis. 
Solution. Here, (Z) is a circle of radius a with centre 
at M, so that S = na*. 
The desired density of the circulation is 


- 1 А 1 
nte ti oe 


2л 


= lint ә | (— a?) sin? £ dt = — 1. 
On the other hand, 
i j k 
д д 
сиг] а= Jz ду ds = —k 
у 0 0 


and, hence, 
(curl a, n°) |м = (—k, К) = —1, 


which, by virtue of (2), confirms the correctness of the 
result. 

204. Compute the density of the circulation of the vector 
а = 21 + 2) + yk around the circle L: {у = a cost, 
z=asint, z = 0, (0 < t< 2n)} at the centre of the 
circle, C (0, 0, 0), in the positive direction of the z-axis. 
205. Compute the density of the circulation of the vector 
а = 2yi + 5zzj around the ellipse Г: {хт = a созі, y = 
= b sin t, z = 1(0 x t < 2л)} at the centre С (0, 0, 1) 
of the ellipse in the positive direction of the z-axis. 
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Sec. 17. The independence 
of a line integral of the path 
of integration. Green's formula 


Definition. A region G of three-dimensional space is 
said to be simply connected (more precisely, it isa simply 
connected plane region) if any closed contour lying in the 
region can be spanned by, a surface lying entirely in G. 
For example, the whole of.three-dimensional space and 
the interior of a sphere are simply connected regions; the 
interior of a torus and three-dimensional space with 
a straight line deleted are not simply connected regions. 

Theorem. In order that the line integral 


| (a, dr) = j P (z, у, 2) dz +Q (т, y, z) dy +R (т, y, 2) dz 
1. L 


should not depend on the form of the path of integration L, 
it is necessary and sufficient that the vector field 


а = Р (=, У, z) i + Q (z, y 2) j + R (z, y z)k 
be irrotational, that is, 
curl a (M) = 0. 
It is assumed here that the coordinates P (z, y, 2), 
О (z, у, 2), R (z, y, 2) ofthe vector a have continuous 


partial derivatives of the first order and the domain of 
definition of the vector a (M) is simply connected. - 


In that case, the line integral | (a, dr) will depend 
solely on the position of the initial and terminal points 
of the path of integration 1. 

If the theorem is complied with, the circulation oj 


the vector a (M) around any closed contour C located in 
the field of the vector a(M) is equal to zero: 


фка, dr) = 0. 


С 
Example 1. Show that in the field of the vector 


a = 27 + хул] + Toy 
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the line integral | (a, dr) is independent of the shape 
L 
of the path of integration L. 

Solution. The coordinates of the vector a are everywhere 
continuous functions so that the domain of definition G of 
the vector a is the entire space (a simply connected region). 
In this region we have 


i j k 
д д д 
curl a=| 9 ду д: |=0 


zy?z xyz iy 
Consequently, the line integral 


je. dr) — E dz + z?yz dy + i a?y? dz 

L 

is independent of the shape of the path of integration L. 
In particular, for the plane vector field 


a (М) = P (z, y) i О (z, y)j (1) 
we have 
i j k 
curl a(M)= Xxx = (3—32). 
Р О 0 


Therefore, for the plane vector field (1) that is defined 

in asimply connected region G, the condition curl a ( 

= 0 is written in coordinate form thus: дР/ду = dQ дг. 
This means that in order јог а line integral 


ki y)dz+Q (2, y) dy 


in a plane field defined in a simply connected region G to 
be independent of the shape of the path of integration, it 
is necessary and sufficient that the following relation hold: 


aP 20 
у  0:° 

Remark. The requirement that the region G, where the 
vector a = a (M) is defined, be simply connected is essen- 
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tial. If the region G is nonsimply connected, then, provid- 
ed that curl a (M) = 0, the line integral may depend 
on the shape of the path of integration. 


Example 2. Let us consider the line integral 
—ydz zdy 
eet si» 


Solution. The integrand is meaningless at the point 
O (0, 0); and so we eliminate that point. In the remaining 
part of the plane (which is then a nonsimply connected 
region), the coefficients of dz and dy are continuous and 
have continuous partial derivatives, and the following 
identity holds true: 


д т __ д y 
9: (=>) = ду ( — nia) , 
On the other hand, if we compute the integral around the 


circle L: z? + y? = R?, then by parametrizing the equa- 
tion of the circle we get - 


2n 2л 
—yd d Н? sin? t+ R? cos? t 
= Ер n Ami M у (ш 22. 


We find that the circulation is nonzero and, hence, the 

line integral depends on the path of integration. 
Determine in which of the vector fields indicated below 

the integral is independent of the shape of the path of 


integration: 
206. а = 24 + z?j + yk. 
207. а= itu te 


утрату, 
_ yi—zj+zk 
208. а= руа 


Green's formula 
Given in a region D with boundary L a plane vector 


field 
a = P (z, у)і+0 (5, у)і, 


where the coordinates P (z, y), Q (z, y) are continuous 
and have continuous partial derivatives ОР/ду, 0Q/Oz. 
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Then Green’s formula, 
фраг Qay= | TE) 9244. (2) 
D 


holds true. Here, the boundary L is traversed so that the 
region D is on the left. 


— 


— 


Li 
Fig. 33 


The region D may also be nonsimply connected so that 
the boundary may consist of several components (Fig. 33). 
In that case, the integral э 


JPas+ Ody 


is understood to mean the sum of the integrals over all 
components of the boundary of D. 

Green's formula (2) is a special case of the Stokes theo- 
rem (see Sec. 16). 

In certain cases, Green's formula permits simplifying 
the computation of the circulation of a vector field. 

Example 3. Compute the circulation of the vector 


a—y 1-4 2 - у? 1+ ylzy + № (c+ УТ 223 y?)]j 


around the circle z? + y? = R. 
Solution. The circulation of the given vector is 


C= фка, аг) = ф VIF +y dz 
L L 


у [зу № (z+ Ут 22 у )] dy. 
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Here 


P=Vitety, Q=aytyln(c+Vite+y*). 
We find the partial derivatives 
aE, ANO 
A  yipsyg 9: 57 
Using Green's formula, we obtain 


—— ПИСЕ 
VIAF) 


c- утрат утрати) 9 


= | jy dz dy. 
D 
Passing to polar coordinates, 


= р соз ф, y=psing, 
we have 


C= | \p2sin? фр ар аф= | | р? sin? pdpdg. 

| [| 

Since 0 < 9 < АЫ О<о = R, it follows that 
C= Jintao ot = 


Compute the circulation of the following vectors around 
the given contours using Green's formula: 
209. а = (у + z)i + у— 2) }; L: х+у=1 ғ = 0, 
y — V. 
210. а = (z — y?) i + 22уј; L: y =z, y = 22. 
23у 
211. а = = 1р (1 ар ний L: а? + у? = 22. 
212. а = уі — 22]; L: = —1, 2 = 0, у = 0. 
(Ву V 1 21-448) 1+ (1вув 2° Viet tay) | 
213. a= = ———— 
3y id :1+ 49 
L: 22+y?=1. 
214. Use Green’s formula to compute the difference be- 
tween the integrals 
I= | (a ye dz— (z— dy 


AmB 
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and 
= | (z+y)}}dz—(z— y)dy, 
AnB 


where AmB is a line segment joining the points A (0, 0) 
and B (1, 1), and AnB is an arc of the parabola y = 22. 
215. Prove that the integral 


(2z +y) dz -- 2x dy, 


L 


where L is a closed contour, yields the area of the region 
bounded by that contour. 
216. Using Green's formula, compute the line integral 


(a, dr) in the vector field 


а = (ех sin y — y) i + (ех cos y — 1) j, 


where the curve L is the upper semicircle z? + у? = 2x 
traversed in the direction from the point A (2, 0) to the 
point O (0, 0). 


CHAPTER IV 
POTENTIAL FIELDS 


Sec. 18. The criterion for the potentiality 
of a vector field 


Definition. A vector field 
а (М) = P (z, y, z)i + Q (z, у, 2)} + В (т, у, 2) К 
specified т a spatial region У is said to be potential И 
there exists a scalar function q (M) such that at all 
points of V the following equality holds: 
a (М) = grad o (M). a) 

The function q (M) = ọ (z, y, z) that satisfies (1) in 
V is termed the potential (or the potential function) of 
the vector field a. 

The relation (1) is equivalent to the following three 
scalar equations: 

OP: =# = 2Ф 
Pir y=, 0(2,у,2)= 50, Ну 1) = 5. 
(2) 

The potential of a field is not defined uniquely but 
only up to an additive constant. 

Remark. For force fields, the function ф (M) is ordi- 
narily called a force function, and the potential is the func- 
tion —9 (M). 

Example 1. (The electrostatic field of a point charge.) 
Show that the field of electric intensity E set up by a point 
charge 4 located at the coordinate origin, 


P Е=- г, г= У 22+ у2 + 22, 


is a potential field. 
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Solution. The problem is posed thus: show that there 
exists a function ф (т, у, 2) such that relations (2) hold. 
In our case, we have 


P(zy, 2-25, Oy 294. Rey, 1) = 4. 


Since 


à (1 1 er z 
x)- Dmoz Pr 
and, analogously, 
£)--Ri)-- 


it follows that the function 
— Dr cM, PER 
P (т, y, 2) — r y + ул 
is the potential of the given field: 
grad ( —+) =E. 


In this case, the coordinate origin (this is where the 
charge q is located) is a singular point of the field E. 
Theorem. Рог a vector field а (M) specified in a simply 
connected region V to be potential, it is necessary and suf- 
ficient that the following condition hold at every point of V: 


curl a = 0. (3) 


In other words, for a vector field specified in a simply 
connected region to be potential, it is necessary and suf- 
ficient that it be irrotational. 

The potential ф (z, у, 2) of the vector field 

а = P (z, у, 2) 1 + О (z, у, z2j+RG, у, 2) К 

is defined by the formula 

(x, y, 2) 
ey 2- | Pdr+Qdy+Rdz, (4) 

(хо, Vo. 20) 

where (zo, Yo, Zo) is some fixed point of the field and 

(r, y, z) is an arbitrary current point. 

Example 2. Show that the field of the vector 
а = 221 + у} + 2k 


is a potential field. 
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Solution. The coordinates P = 22, Q = y*, R = z* of 
the vector a are infinitely differentiable functions through- 
out the space so that a is an infinitely differentiable 
vector defined throughout three-dimensional space. We 
have 


i j Е 
д д ô д д 
сиг] а= Jz ду 5+ = (52—95)! 
22 у 22 
д д д 
+(¢2-Z2)1+(Zy—-Ze)k=0 


By virtue of the theorem given on p. 122 the field of the 
vector a is a potential field. It is readily seen that the 
function 


3 
9(z, y, = EET + с, 


where C is an arbitrary constant, is the potential of the 
given field. 

Check to see whether the following vector fields are 
potential fields: 
217. a = zzi + 2yj + хук. 
218. а = (22у + z?) i + (2yz + 23) j + (222 + y?) К. 
219. а = 4 (zi + yj + 27%. 


220. а = yz cos zy-i + zz cos zy-j + sin zy-k. 
221. a = In (1 + 22) i + In (4 + 2?) j + zzk. 
a 


a a= (3c 3) re ( +1) ie ee 


z 
223. H= ( yi o zi), Пур, гео. 
224. Prove that the field a = f (r)-r, where f(r) is a 
differentiable function, is a potential field. 
225. Show that the vector lines in the potential field 


а = grad @ are perpendicular to the level surface of 
the function q. 
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Sec. 19. Computing a line integral 
in a potential field 


Theorem. A line integral in a potential field a (M) is 
equal to the difference between the values of the potential 
Ф (M) of the field at the terminal and initial points of 
the path of integration: 

M, 

| (а, dr) =9 (M2) —9(M)). (1) 

M, 

Example 1. Compute the line integral in the field 
of the vector 
r = zi + yj + zk 

along the line segment bounded by the points M, (—1, 0, 
3) and M, (2, —1, 0). 

Solution. We will show that the field of the given vector 
is a potential field. Indeed, 


i j k 


It is easy to see that the potential of this field is 


9 (т, у, = Stet? үс, 

Using formula (1), we obtain 

м, 

f а, dr) =9(2, —1,0—9(—1,0,3) 2 2—5- —$. 
Mi 


Note that it is immaterial what line joins the points M, 
and М,; for fixed M, and M, the integral 


M, M, 
| (a, dr) = | тах + уду-+ zdz 
M, M, 


always has the same value. 
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Computing the potential of a field in Cartesian coordinates 


The formula 
ez, Y, 2) 
(z, V, 2) 


= j P (z, y, 2) dz +Q (т, y, 2) dy + (т, у, 2) dz 
(Xo, Ио, 20) 
(2) 


may be used to find the potential function q (M) = 
= Фф (х, y, 2) of a specified potential field 
а (z, y 2) = Р (х, у, 2)1+ 0 (2, у, 2) 
+ В (z, у, 2) К. 


To do this, fix theinitial point М, (zo, yo, Zo) and join 
it with the current point M (z, y, z) with a broken line 


M(2,Y,2) 


Fig. 34 


M,ABM whose segments are parallel to the coordinate 
axes; namely, M,A || Oz, AB || Oy, BM || Oz (Fig. 34). 
Then formula (2) takes the form 

x у 
g (z, у, 2) = | P (z, Yo, 20) dz Í Q (z, ys zo) dy 


Xe Vo 
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+ {в (z, у, 2)аг, (3) 
1. 


where z, y, z are the coordinates of the current point on 
the segments of the broken line along which the integra- 
tion is performed. 

Example 2. Prove that the vector field 


а = (y + 21+ (z + 2)ј + (= + ук 
is a potential field and find its potential. 
Solution. Ist method. A necessary and sufficient condi- 


tion for the potentiality of a field a (M) is that curl 
a (M) be zero. In our case, 


i j k 
curl a= ыйы Bs, 9; 
= Oz ду д2 
у+2 2+2 ху 


= (1—1) 14 (1—1) } + (1—1) К =0. 
That is, the field is a potential field. The potential of 
the field can be found with the aid of formula (3). For 
the initial fixed point we take the coordinate origin 
О (0, 0, 0). We thus have 


x Р у 2 
e. y, 2) = | (0+0) dz+ [(=-0)и+ | (2+9) az 
0 0 0 


= ту + 12 + 42. 
To summarize, 
9 (2, y, 2) = zy + zz + уг + С, 
where C is an arbitrary constant. 
2nd method. By definition, the potential ф (z, у, 2) 
is a scalar function for which grad ф = a. This vector 
equality is equivalent to three scalar equalities: 


I y+, (4) 
Bartz, (5) 
arty. (6) 
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Integrating (4) with respect to z, we obtain 
Ф (т, у, 2)=\(y+a)de=zy+22+fly,2), (7) 
0 


where f (у, 2) is an arbitrary differentiable function of 
y and z. Differentiating both sides of (7) with respect to 
y and taking into account (5), we obtain a relation for 
finding the as yet undetermined function f (y, 2). We 
have 


ap — 
wt» 
or 
zii. 
whence 
of 
== (8) 


Integrating (8) with respect to у, we have 


у 
fy. 2) = |2dy=zy+F (а), (9) 
ó 
where Р (2) is ап as yet undetermined function of 2. 
Substituting (9) into (7), we get 


9 (z, y, 2) = zy + zz + zy + F (2). 


Differentiating this equation with respect to z and taking 
into account (6), we obtain an equation for finding F (z): 


z4y-ityi, 


whence dF/dz = 0 so that F (2) = C = constant. 
Thus we have 


ф (z, у, 2) = ху + yz + zr + С. 


3rd method. By the definition of the complete differen- 
tial of the function ọ (z, y, 2) we have 


do— 59 dz 4 2 dy + Baz. 
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Substituting in place of the partial derivatives дф/дх, 
дф/ду, Op/dz their expressions taken from (4), (5), (6), 
we obtain 


dp = (y + 2) dz + (z + 2) dy + (z + y) dz 
or, after some simple algebra, 
dp = (y dz + z dy) + (z dz + z dz) + (y dz + z dy) 


= d (zy) + d (zz) + d (yz) = d (zy + zz + yz). 
Thus 
de = d (zy + yz + zz), 
whence it follows that 


Фф (z, у, 2) = ху yz + zz + C. 


In the following problems, establish the potentiality of 
the given vector fields a (M) and find their potentials 


Ф (М): 
226. a = 2zyzi + z?zj + Tyk. 
227. a = (yz + 1) i + zzj + zyk. 
228. a = (2zy + 2) i + (22 — 2y) j + zk. 
_ АК 
229. a= TII 
i k 
230. а= nca 


231. а 


r 
T 
232. a=. 


233. a=r-r. 
When the region Q is a star with centre at the coordi- 
nate origin О (0, 0, 0)*, the potential ф (М) ої some 


* The region Q is called a star-shaped то relative to the 
point O belonging to Q if any ray emanating from that point cuts 
the boundary of the region at one poirit at most. For example, star- 
shaped regions in the plane are the plane itself, a parallelogram, 
a circle; in three-dimensional space, the space itself, a parallele- 
piped, a sphere. 
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vector field a = a (M) at the point M (z, у, 2) may 
be found from the formula 


1 
Ф(М)= fa (M^), r(M)) dt--C, C=constant, (10) 
П] 


where r (M) = zi + yj + zk is the radius vector of the 
point M (z, у, 2), and the point M' (tz, ty, tz) forO < t 
< 1runs over the segment OM of the straight line passing 
through the points O and M. 

Example 3. Find the potential of the vector field 


а = yzi + zzj + zyk. 


Solution. It is readily seen that curl a = 0, which 
means the given vector field is a potential field. This 
field is defined throughout three-dimensional space and 
is star-shaped with centre at the coordinate origin 
O (0, 0, 0), and therefore to find the potential we take 
advantage of formula (10). Since in this case 


a (M') = a (tz, ty, tz) = t?yzi + t?zzj + t?zyk, 
it follows that the scalar product of the vectors a(M’) 
and r (M) is equal to 
(a (M^), r (M)) = t? (туз + хуз + zyz) = 3t?zyz. 
The desired potential is 
1 1 
Ф(М) = [с (M^), т(М)) dt = туз Var dt + C=zryz+C. 
0 0 


Ф (M) = zyz + C. 


Using formula (10), find the potentials of the following 
vector fields: 

234. а = ai + Bj + yk, where a, В, y are constants. 
235. а = (y + 2)i + (z + 2) ј + (у + 2) К. 

236. а = уі + zj + ek. 

237. а = e* sin y-i + е cos y: j + К. 


Thus, 
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CHAPTER V 


THE HAMILTONIAN OPERATOR. 


SECOND-ORDER DIFFERENTIAL 
OPERATIONS. 


THE LAPLACE OPERATOR 


Sec. 20. The Hamiltonian operator del 


Many operations of vector analysis may be written in 
abbreviated form and in a form convenient for calcula- 
tions; this is done through the use of a symbolic operator 
called the Hamiltonian operator del: 


+i eke (1) 


This operator combines both differential and vectorial 
properties. We will regard the formal multiplication of 
ólóz by the function f (т, у, z) as the partial differen- 
tiation д} /óz. 

Within the framework of vector algebra, we will per- 
form the formal operations involving del as if it were 
a vector. Using this formalism, we obtain the following. 

1. Ifu — u (z, у, 2) is a scalar differentiable function, 
then, by the rule of multiplying a vector by a scalar, we 
have 


• д . д д 
Уи = |1; К) и 


ди 


. Ou , ды 
=i g ti gy tE = габи. (2) 


2. If a = P (z, y, 2) і + Q (z, y, z) j + R (z, у, 2) К, 
where P, Q, R are differentiable functions, then, by 
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the familiar formula for a scalar product, ме have 
à 1.0 д ; ; 
(у, а) = (13r tig tk Pit Qi+ Rk) 
Фр, àQ 8. e 
=; + зу t = У 8 (3) 
in particular, (V, с) = 0, where c is a constant vector. 


If a = P (z, у, z)i+ Q(z, y, z)j+ R (т, У, z)k, 
then 


[V,a] =|= —— —— |— curl a; (4) 


in particular (V, с) = 0, where е is a constant vector. 

Continuing the formalism of operations with V as 
a vector, we obtain the following from the distributive 
property of scalar and vector products: 


(V, a + b) = (V, a) + (V, b), (5) 
that is, div (a + b) = div a + div b 
(V, а + b) = [V, al + [V, Ы (6) 


or curl (a + b) = curl a + curl b. 

Formulas (5) and (6) may also be interpreted as an 
exhibition of differential properties of the del operator 
(V is a linear differential operator). 

When using the formalism of operations involving the 
del operator regarded as a vector, one must bear in mind 
that del is not a vector, for it has neither magnitude 
nor direction, so that, for example, the vector [V, al 
will not, generally speaking, be perpendicular to the 
vector a (however, for the plane field a — P (z, y) i 4- 
+ Q (z, y) j the vector 


[V, a] 2 curl a= (3-25) k 
will be perpendicular to the zy-plane and, hence, to the 
vector a). In the same way, the concept of collinearity is 


meaningless with respect to the symbolic vector V. For 
example, the expression [Уф, Уф], where Ф and : are 


9* 
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scalar functions, formally resembles a vector product of 
two collinear vectors, which product is always equal to 
zero. But this is not true in the general case. Indeed, the 
vector Уф = grad ф is directed along the normal to the 
level surface Фф = constant, while the vector Vip = grad p 
defines the normal to the level surface ф = constant, 


Fig. 35 


and in the general case these normals are not necessarily 
collinear (Fig. 35). On the other hand, in any differentiable 
scalar field p we have [Уф, УФ] = 0. These examples 
show that the del operator must be handled with care. 

Besides its vectorial nature, the Hamiltonian operator 
del has a differential aspect to it. Taking into account 
the differential aspect of V, we will agree that the operator 
V acts on all quantities that follow it. In this sense, 
(V, а) == (а, V). Indeed, 

(V, a) = div a, 

whereas 


д д д 


is a scalar differential operator. 
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When applying the del operator to a product of any 
quantities, one must bear in mind the rule for differenti- 
ating a product: 


д ди ди 
эх (0) 0 52 Ни аг 


From this it follows that the del operator must be applied 
in succession to each factor, leaving the other factors 
unchanged, and then the sum of the resulting expressions 
taken. In this procedure we are guided by the following 
rules. 

1°. If the del operator acts on some product, first 
take into account its differential nature and only then 
its vectorial properties. 

In order to note the fact that the del operator 
does not act on some quantity involved in a complicated 
formula, that quantity is labelled with the subscript с 
(constant), which may be removed in the final result. 

3°. All quantities not acted upon by the del operator 
are placed in front of the operator (that is, to the left) 
in the final result. 

Example 1. Show that 


div (ua) = и div a + (a, grad и). 
Here, u is a scalar function and a is a vector function. 
Solution. In symbolic notation, 
div (ua) = (V, ua). 
Taking into account first the differential nature of V, 
we have to write 
(V, wa) = (У, и.а) + (V, иа,). 
Regarding the expression (V, u.a), we сап take the con- 


stant factor u, outside the del sign and, as a scalar, outside 
the sign of the scalar product; this yields 


(V, ига) = (u.V, 8) — ш (V, а) =u (V, a). 


(The last step is to drop the subscript c). 
In the expression (V, wa,), the del operator acts only 
on the scalar function u, and so we can write 


(Vi иа.) = (Vu, а.) = (а,, Vu) = (a, Vu) 


134 Second-Order Differential Operations (Ch. V 


to get the formula 
(V, ua) = и (V, a) + (a, Vu) 
or 
div (ua) = u div a + (a, grad и). 
Example 2. Show that 
curl (ua) = и curl a — [a, grad ul. 
Solution. In symbolic notation, 
curl (ua) = (V, ua). 
Taking into account the differential properties of V, we 
first write 
[V, ua] = [V, u.a] + (V, иа,]. (7) 
Then in the first term on the right we take the scalar 
factor u, outside the del sign and outside the sign of the 
vector product, which yields 
[V, ua] = и, [V, а] = u [V, al. 
In the second term in (7) we refer u to the operator V and 
change the order of the factors so that the vector a,, which 
del does not act on, is in front of V. This yields 
[V, ua,] = [Уи, а„] = —[а, Vul. 
Thus 
[v, ua] = u [v, а] — la, Уи] 
or 
curl (иа] = и curl a — [a, grad и]. 
Example 3. Use the symbolic method to find div [a, b] 
Solution. We have 
div [a, b] = (V, [а, bl) = (V, (а, b.]) 
+ (V, la., b)). (8) 
Using the property of cyclic permutation of factors in 
a mixed product, we transform the expression on the 
right of (8) so that all constant quantities are in front 
of the del operator and all variable quantities follow it. 
This yields 
div[a, Ы = ([У, а], b.) — (V. (Ы, a.l) 
= ([У, а], b.)—([V. b]. а.) &(b. [V, a]) —(а, [V, Ы) 
or div{a, b] = (b, сиг а) — (a, curl b). 
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Remark. The use of the symbolic method enables us to 
avoid cumbersome analytical transformations and obtain 
the final results very quickly. On the other hand, the 
various formal transformations involving the del opera- 
tor must be performed with extreme caution in order to 
avoid serious mistakes. For this reason, if there is any 
doubt about the final result, it is wise to verify it by 
the analytical method. 

238. Show that 

(a) у (2) = 22"; 

(b) Vf (u) = f' (u) vu. 

239. Prove that the vector [Vu, Vv] is solenoidal if 
u and v are differentiable scalar functions. 

Use the Hamiltonian operator del (V) to prove the 
following equations: 

240. (a) grad (uv) = v grad и + и grad v; 


(b) сог! (а, b) = (b, У) a (а, V) b + a div b — 


— b div a. 

241. curl [а, г] = 2a, where a is a constant vector. 
242. Prove that the vector a = и grad v is orthogonal to 
curl a. 


Sec. 21. Second-order differential 
operations. The Laplace operator 


Second-order differential operations are obtained as 
a result of a twofold application of the operator V to 
fields. 

Suppose we have a scalar field и = и (M). In this 
field, the operator V generates a vector field Vu — grad u. 

In the vector field Vu, the operator V, when applied 
a second time to Vu, yields the scalar field 


(V, Vu) — div grad u (1) 
and the vector field 
(V, Vu) = curl grad и. (2) 


If a vector field a = a(M) is given, the operator V 
generates in it a scalar field (V, a) = div a. In the scalar 
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ield div a the operator V generates a vector field 
V (V, 8) — grad div a. (3) 


In the vector field a — a (M), the operator V also gener- 
ates a vector field [V, a] — curl a. If the operator V is 
again applied to this field, we obtain the scalar field 


(V, [У, а]) = div curl a (4) 
and the vector field 
[v, [v, al] = curl curl a. (5) 


The formulas (1) to (5) define what are called differen- 
tial operations of the second order. 

Example 1. Suppose a function u — u (z, y, z) has 
continuous partial derivatives up to second order inclu- 
sive. Prove that 


cur] grad и = 0. 
Solution. 1st method. Operating formally, we obtain 


curl grad и = (у, Vul = [v, Vu = 0 


since [V, V] = 0 being the vector product of two identical 
"vectors". 

2nd method. Using the expressions for the gradient and 
curl in Cartesian coordinates and taking into account the 
given conditions, we have 


i j k 
д д д 
Б т ш д [ди д [ди 
curlgradu= [Эт Эу э |а (э) — 9% (9) |! 
дш ди дш 
дт ду дг 


а) не (2-3 501: 
= ( ди au )i дщ Әд?ц i ( д3и ôu )k=0, 


ду д2: az dy 9:05 Oz az дтду dyor 
since the mixed derivatives are equal in this case, 
au д?и u au д3ц tu 


буд: дтду' 0:02 920: дтду дудт' 
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In similar fashion proof is given that for the vector 
field 


а = Р (=, у, 2) і + О (z, y, 2j+ R (z, y, 2) К, 


the coordinates Р, О, R of which have continuous partial 
derivatives of the second order, we obtain div curl a = 0. 

Note particularly the second-order differential opera- 
tion div grad u = (V, Vu). Assuming that the function 
и (z, у, 2) has second partial derivatives with respect to 
т, у, 2, We obtain 


9 d А a 210 


0) ннем 


Thus, (V, Vu) = Au, where the symbol 
_ ёзи , ди, ди 

А = 3+ LT RET 
is termed the Laplace operator (or Laplacian). It may be 
represented as a scalar product of the Hamiltonian oper- 
ator V into itself, that is, 

a no 
A = (V, У) = У? = дз tas ая. 

This operator plays an important role in mathematical 
physics. 

Let us examine another second-order operation curl 
curl a. We have curl curl а = [V, (у, al]. Let us take 
advantage of the formula for a vector triple product writ- 
ten as (A, [B, Cl) = B (A, С) — (A, B) C. Replacing 
A by V, B by V, and C by a, we obtain 


(v, [v, al] = v (v, а) — (V, V) а=У (V, а) — Ла (6) 


or 
curl (curl a) = grad div a — Aa, where Aa = 


= ДР.1 + AQ-j + AR-k. 


The following table is a pictorial display of second- 
order differential operations: 
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Sealar f!eld u Vector field a 


grad div curl 


E grad div a 

ЕЯ 

> div grad u=Au |  — — ^| div curl a=0 

T du=0 curi curl a 

5 шый enn =grad div а—Ла 


Example 2. The laws of the classical theorem of electro- 
magnetism are postulated in the form of a system of 
Maxwell equation. 

In the most elementary case of a nonconducting homo- 
geneous and isotropic medium and in the absence of 
charges and currents, this system is of the form 


IO =[у, Н). (7) 
-+ % у, E], (8) 
(V, Е) = 0, (9) 
(V, Н) =0. (10) 


Неге, Е апа Н are vectors of the electric-field and mag- 
netic-field intensity; e and p are coefficients of the permit- 
tivity and permeability (our assumptions are that e 
and p are constants); c is the velocity of light in empty 
space. 

Since the spatial and time derivatives commute, that is, 


sv. Bi-[v. 52], 
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it follows that by differentiating (7) with respect to t, 
we obtain 
e ó*E -[ oH ] 


v of L'' at 
Replacing A taken from (8), we get + ZE = 
= - £19, Iv, Ell or 
4 SF гу, IV. Ell. (41) 
By virtue of formula (6), [V, |v, Е] = v(V, E) — AE. 
Since (V, E) = 0, it follows from (11) that AT = — AE. 


To summarize, for the vector field E we obtain the 
equation 
OE c? 
ET — en 
This is one of the basic equations of mathematical physics 
and is called the wave equation. 
It is easy to see (check this!) that the vector field H 
satisfies the same kind of wave equation 


aH e 
Ten. mr AH. 


Thus, under our conditions, each of the coordinates 
E, Ey, E, and H,, Ну, H, of the vectors E and Н 
satisfies the equation 


Au Bu u 
— = о? | — чырай rum. 
ot? a (tat ar) 


Неге, а = с/у ep is the velocity of propagation of the 
process. In a vacuum, where = = р = 1, we have а = c, 
that is to say, in a vacuum, electromagnetic processes are 
propagated with the velocity of light. 

243. Show that any solution of the equation [V, [V, Al] — 
— РА = 0 that satisfies the solenoidal condition satisfies 
the Helmholtz vector equation 


VA + РА = 0. 
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Definition. The scalar field и = и (т, у, 2) that satis- 
fies the condition Au = 0 is termed a Laplace (or karmon- 
ic) field. 

Example 3. An important instance of a harmonic field 
is the scalar field u=k/r, k—constant, r= V z? + y? + 22. 
This function is the potential of a gravitational 
field generated by a point mass placed at the coordinate 
origin. It is*easy to verify that the function u = k/r is 
harmonic everywhere except at the coordinate origin, 
where it is not defined. Indeed, 


= 2 (m, г)—-% (V, г) = 
for all г 520 since 
1 1 
(v, M=(V,)=(V—, г)+-—(У,г) 


пазы сайыш: 


Example 4. Prove that in the potential field of a vec- 
tor a its potential function u (т, у, 2) satisties the Pois- 
son equation 


au dtu, at 
Аи == Fat Буз ac Pn 2), (12) 
where р (т, у, 2) is the divergence of the vector a. 
Solution. It is given that 
div a — p. (13) 


Since the field of the vector a is a potential field, 
it follows that a — grad u, where u is the potential of 
the field. Substituting a — grad u into 13), we obtain 
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div grad и = р or, since div grad и = Au, we have 
Au = р. 

In the special case of points of the field where the diver- 
gence is zero, equation (12) turns into the Laplace equa- 
tion Au = 0. The Laplace-Poisson equation permits find- 
ing the potential function и by integrating the partial 


В № 


Ух | +++ +++ + + 


РРР * La aaa ahhh ЛУ! 


4 0 Ay 
Fig. 36 


Wy 


differential equation. This turns out to be more convenient 
in some problems. 

In electrostatics, preference is often given to finding 
the function v = —u instead of и. Then a = —grad v. 
Accordingly, in the theory of the electrostatic field, the 
Poisson equation is of the form 


д% 0% ду 4лр 
att oat ae — е a) 

Let us consider an elementary case where the Poisson 
equation is employed. 

Example 5. Suppose two infinite parallel plates AA, 
and BB, with opposite charges have potentials v, and v,; 
for the sake of definiteness, v, > v,. Find the field E 
between them (Fig. 36). 

Solution. Send the z-axis at right angles to the plates 
in the direction of decreasing potential, and bring the 
yz-plane to coincidence with the positively charged plate 
AA,. We now seek the potential function from the Pois- 
son equation. By virtue of the symmetry of the problem 
about the z-axis and due to the infinity of the plates, 
we can conclude that the equipotential surfaces are planes 
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parallel to the plates, and the function v depends solely 
on the variable z. Equation (14) takes the form 

a? 

30 (45) 
Since space charges are absent throughout the space be- 
tween the plates. Integrating (15), we find 

v = Cz + С», (16) 

where С, апа C, are arbitrary constants. 

We require that for z = 0 the function v take the value 
v, and for z = d, where d is the distance between’ the 
plates, that it take the value v,. This yields C,=v,, v= 
= Cyd + С», whence C, = Vi» C, = (v, — v,)/d. Substi- 
tuting these values of C, and C, into (16), we obtain 

Ug — 01 


v= T. 


Vy— и. 
І+0 = 0, — = i 


The vector E is found from the formula E = —grad v, 
which yields 


Е= 2-2% 1 


so that the field is homogeneous and in the direction of 
the z-axis. The magnitude of E is equal at every point to 
| E | = (v, — v,)/d, that is, it is equal to the potential 
drop per unit length of the shortest distance between the 
plates. 

244. Suppose a scalar function q (M) satisfies the Laplace 
equation. Show that the vector уф is solenoidal and irrota- 
tional. 

245. Show that A (uv) = u Av + v Au + 2 (Vu, Vv). 
246. Prove that if r is the radius vector, then 


E in space, 
Ar= 4 
. in the plane. 


247. Check to see whether the following scalar fields 
are barmonic or not: 

(a) u = T? +4 22у — у?, 

(b) u = zy + ys + 222, 

(c) u = 1? — }%. 
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248. Show that the scalar field 
u=in+, where г= И 22 у? (r0), 


is harmonic. 
249. Find all harmonic fields that depend solely on z. 
250. Find the general form of a homogeneous harmonic 
polynomial of second degree in z and y. 
251. Find all solutions of the Poisson equation Au = 
= 1"? that depend solely on z. 

Example 6. Green's formulas. Let Фф, p be two scalar 
functions of a point. Set up the vector a = q grad v. 
Then 


div а = div (Ф grad p) = ф div grad p + (grad Ф, grad yp) 
= Фф. Ap + (grad Ф, grad $). 
Now apply the Gauss-Ostrogradsky formula 


9 (a, по) do = }| | div adv. 


Note that in our case 


(a, n*) = (o grad р, n°) = ọ (grad p, п) =p, 
We thus obtain Green's first formula: 


| | f irav (де, grad LE LE Фе de, (47) 
which for ф = tp turns into 
1] Ie Ap +lerad oP = ө 99 до. (18) 
If we put p = 1 in (17) we get 
jj | Av do = № ао. 


In (17), interchange Фф and їр and subtract the resulting 
formula 


1j [p Ap + (grad sp, grad g)) йр = 099 99 do 
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from (17). This then yields Green’s second formula: 


Jf J ay eame ( p 2* — yp 2) до. 
/ r 


It is assumed here that all functions that we have to deal 

with and also all their derivatives that occur in the for- 

mulas are continuous in the region under consideration. 
Example 7. Find the surface integral 


I= Dez 


taken over the sphere i: 2 + у + 22 = 1 for ф = 
= 4? + у? and p = у? + 22. 

Solution. By Green S first formula, the desired integral 
is 


I = | | | ФАФ+ (grad o, grad $) dv, 
V 
where the region of integration V is a sphere: z? + y? + 
+ 22 < 1. 
We have Ap = 4, grad ф = 2zi + 2yj, grad p = 
= 2yj + 22k, (grad Ф, grad p) = 4y? and therefore 


I= j || (422 + Ay? + Ay?) dv — 4 11 (22+ 2y?) dv. 


Passing to the spherical coordinates z — r cos q sin 0, 
у = rsin ọsin Ө, 2 = rcos Ө, we obtain 


[=4 | | | (г? cos? ф sin? 6 + 2r? sin? p sin? 0) r2 зіп Ө dr 20 дф 
у 
2л л 1 
=4 | (cos? ф + 2 sin? y) dp | 31020 20 | r*dr 
0 0 0 
2л л 
== f (1+ sin? ф) dp | (1 — cos? 6) sin 0 40 
0 0 


Е f 1 за) |" 16 
=> л ( cos 8+ — cos ө) „^5 7. 
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Example 8. Find the surface integral 


I= $( (o %—% 59) do 


taken over the NS У: г? у —-HBH,.uz-0,;—H, 
Н > 0, provided ф = 2? + y? dicta zy + 
+ 22 + =. 

Solution. By Green's second formula, the desired inte- 


gral is 
1= | | (ФАф—ФА9) dv. 


For the given functions ф and we have Ag = 4, Ay = 4 


and therefore 
I= —4 | № 2 dv. 


Passing to the cylindrical coordinates z = р cos ф, у = 
=p sin Фф, 2 = 2, we obtain 


==. аа-а io | pap | zaz 


| = — 2n R?H?. 
Example 9. Find the surface integral 


I= Goz -2Ф do 


over a closed surface 2 pounded by the planes z + у + z= 
=1, z = 0, y = 0, 2 = 0, provided ф =e sin y + 1. 

Solution. The given function is harmonic since Ag = 
= е“ sin y — ех sin y = 0. Therefore, by (18) we have 


ja Jj | grad ф [2 dv. 
We find the modulus of the gradient of the function q: 
grad ф = е sin y-i + e cos y-j, | grad p | = © 
The desired integral is “ше to 
1 -х- 
I= e? dv= | e? dz ax а | dz= — 9). 
ее | 


10-910 
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252. Compute the surface integral J = 9 p A do over 


a closed surface 2: саа" у=0, y>0}, 
provided ф = 22, ф = 22 + y? — 22. 


253. Compute the surface integral I = 6 (Ф x — y I ) da 


E 
taken over the entire surface of a closed cylinder 2: 
(z?--y?—1, z=0, z—1), provided ф= 222,  — 22 + 22. 


254. Compute the surface integral J = Dem dO, pro- 


vided ф= (z + y + z)/ V Запа Х isa sphere: Poder А?. 
255. Find the surface integral J = 9x do, provided 


p=e*siny+e"sinz+z and 2 is а triaxial ellipsoid 
(22/а?) + (y?/b2) + (z2/c?) = 1. 


Sec. 22. Vector potential 

Let а vector field 
a = a (M) = P (z, y, 2) i + Q (z, y, 2) j + R (z, y, 2) k 
be solenoidal in the region G, that is, div a (M) = 0 in G. 

Definition. The vector potential of a vector field a = 
= a (M) isa vector b (M)=P, (z, y, 2) i+ Q, (т, y, 2) j + 
+ В, (z, y, 2) k that satisfies in С the condition 

curl b (M) = a (M) (1) 

or, in coordinate form, 
дв, _ д0, _ ôP, OR, _ 20 2P, 
ри Uu Uo uua TRE 0) 

For a solenoidal vector field a (M) the SUM potential 
b (M) is not defined uniquely: the condition (1) is also 
satisfied by the vector B (M) = b (M) + grad } (M), 
where f (M) is an arbitrary differentiable scalar function, 
since curl (grad f (M)) = 0. 

Thus, two vector potentials of the solenoidal field a (M) 
differ by the gradient of the scalar field. 
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Finding the vector potential b (М) of the solenoidal 
field a (М) reduces to finding some particular solution of 
system (2) of three partial differential equations for the 
three unknown functions P, (т, y, 2), О, (т, у, 2), 
R, (z, y, 2). 

The vector potential b (М) may be constructed in the 
following manner: Taking advantage of the arbitrariness 
of choice of the vector b (M), we will simplify matters by 
setting P, (z, у, z) = 0, that is, the vector b (M) will 
be sought in the form b (М) = Q, (z, у, 2)ј + 
+ R, (z, y, z)k. Then the system of differential equa- 
tions (2) for finding the unknown functions Q, (z, y, 2) 
and В, (z, y, 2) takes the form 


D gQ, (3) 


From the second and third equations of this system we 
find 


R, (z, у, 2) — | Q (z, у, 2) 2С, (у, 2), 
Qi (z, у, а)= | R (z, у, 2) dz Civ, 2) 


where С, (y, z) and C, (y, 2) are any differentiable func- 
tions of y and z. For the sake of simplicity, set C, (y, 2) = 
= 0 and choose the function C, (y, 2) so that the first 
equation of system (3) is satisfied. To do this, we sub- 
stitute into the first equation the expressions that were 
зач for О, and Ry: 

ac, Әә 


x | Ole, y, дату, | A(z, у, gd 
=P (z, y, 2). 
From this we get 


S | Q(z, у, 2) dz F7 | R (z, y, 2) ат P(z, y, 2). 


10* 
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It is easy to verify that the right-hand side of this 
equation does not depend on z; this is because div a (M) — 
= 0 in G. 

Integrating the last equation with respect to y, we find 


Су, в) = f [5 | Oley 2) 42 
+ | RG, y, 2942 Р(х, у, 2) ]dy4 Сз(). (4) 


Setting С. (z) = 0 in (4) and substituting (4) into the 
expression for А, (=, y, 2), we get a particular solution 
of system (3): 


P,20, (5) 
0,= | R(z, у, 2) ae, 6) 
= | [э ee у, adeo | Ri, y, adz 


*P( y, )]4у— Q(z. y, ааг. (7) 

The vector b(M), whose coordinates P, (z, y, 2), 
О, (z, y, 2), В, (=, у, 2) are defined by formulas (5), 
(6), (7), is the vector potential since it satisfies the condi- 
tion curl b = a. -- 

Example 1. Find the vector potential b = b (z, y, z) 
for a solenoidal field given by the vector 

a = 2yi — zj 22. 
Solution. We seek the potential b in the form 
b = b (z, y, 2) zT о, (z, у, 2) j + R,(z, у, 2) k, 

where Q, (z, y, 2) and R, (z, у, 2) are found from (7) 
and (8). Since in the given case Р = 2y, Q = —2, R= 
= 2z, we obtain 


О, (z, у, 2) = j 25 ат = 1?, 


В, (=, у, 2) = | 212+ j 2y dy = z2 + уг. 


Thus 
b (z, у, 2) = z^j + (zz + y?) k. 
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It is clear, by direct verification, that curl b = a, 
and, hence, this vector is the vector potential of the given 
field. 

Remark. Due to the arbitrariness in choosing the vector 
b, we could require О, (z, y, 2) = 0 or R, (z, у, z = 0 
instead of P, (т, y, 2) = 0. The system of equations (2) 
and formulas (5), (6), (7) would naturally change. 

Find the vector potentials of the following solenoidal 
fields: 

256. a=i+j+k. 

257. a = 2yi + 2zj. 

258. а = (ех — еу) К. 

259. a = 6,21 + 6zj + 62k. 

260. а = 3y*i — 3z^j — (у? + 2z) К. 

261. a = yei + 2yzj — (2zxyze?* + 23) К. 

If the vector field a — a (M) is solenoidal in the re- 
gion G, which is star-shaped (see Sec. 19, Chapter IV) with 
centre at the coordinate origin O (0, 0, 0) [the field a (M) 
may not be defined at the point O], then one of the vector 
potentials b — b (M) may be found from the formula 


1 
b (M) = | а (М7), v (tdt, (8) 
0 


where r (M) = zi + yj + zk is the radius vector of the 
point M (т, у, 2) and the point М’ (tz, ty, tz) runs 
over the line segment OM as the parameter £ varies 
from 0 to 1. 

Example 2. Use formula (8) to find the vector poten- 
tial of the solenoidal field 


а = 2yl — zj + 2%. 


Solution. The given vector field is defined throughout 
three-dimensional space, which is a star-shaped region 
with centre at the coordinate origin, and so we can use (8) 
to find the vector potential. At the point M' (tz, ty, tz) 
we have 


a (M’) = 2tyi — tzj + 2tzk. 
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We find the vector product 


i j k 
[а(М’), r(M)]2|2ty —tz 2tz 
z у z 


= [ — (2zy + z?) i + (22? — 2yz) j + (2y? + хз) k] t. 
Using (8), we obtain 
1 
b(M) = È | — (22у + 22)i + (222 — 2yz) j + (242+ 22) k] dt 
0 
= — 4 Qus (2? ~y2)j++ (2y? + zz) К. 


It is easy to establish that curl b (M) = a (M). 

Remark. In examples 1 and 2 we obtain different vector 
potentials for one and the same solenoidal field a — 
= Qyi — zj + 2zk: 


by (М) = 2%) + (22 + 4?) k, 
be (М) = —3- (ay + 2%) +5 Gf n) bey Qt 22) Е. 


They differ by a term equal to the gradient of some 
scalar field f (M). This term plays the role of an arbitrary 
constant (when acted on by the curl). It may be repre- 
sented as the gradient of some scalar function f (M). 
Let us find this function in our example. We have 


grad f (M) = b, (M) — b; (M) 
= 4- (zy + 22) 1+4 (2+ 2y2) j++ (222+ y?) k. 
To find the scalar field f (M), use formula (3) of Sec. 19, 


in which we take the coordinate origin O (0, 0, 0) for 
the point (хо, Yo, zo). This yields - 


x V z 
(М) = | 0-dz4- | 922 4+ | + zz у) dz+C 
0 0 0 


= (xy + yz +22) +С, 


where C is an arbitrary constant. 
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Example 3. Find the vector potential b of a magnetic 
field H set up by an electric charge e that is moving with 
a constant velocity v. 

Solution. By the Biot-Savart law, the intensity of the 
magnetic field is 


H (M) =EL (9) 


where r is the distance of point M from the charge e. 

Since H is a solenoidal vector, that is, div H = 0, 
there exists for it a vector potential b such that H = 
= curl b or, taking into account formula (9), 


(ev, т] _ e [v, r] 
curl = 4n án в 


Let ыы rewrite this formula as 


curl = 7— + [e St] +[ev, 4 = se) 
= [0 5] ыз D 
= {lt &(F)] HG) + [№ + (5 


Employing the readily verifiable equation 


curla= [t 2] [5 27] [К 2]. 


we obtain 


x Ч 
—— 


)])- 


1 v 
curl b = —— curl —_, 
án r 
whence 


1 ev 
b = Gn b Up * 

Using formula (8), find the vector potentials of the 
following solenoidal fields defined in star-shaped regions: 
262. a = i. 
263. a = 621 — 15yj + 9zk. 
264. а = 5z?yi — 10zyzk. 
265. а = 2 cos zz j. 


—yi-- zi 
266. a— pr. 22+ у? >0. 


CHAPTER VI 


CURVILINEAR COORDINATES. 
BASIC OPERATIONS 
OF VECTOR ANALYSIS 
IN CURVILINEAR COORDINATES 


Sec. 23. Curvilinear coordinates 


In many problems it is more convenient to define the 
position of a point M in space by three numbers (41, qs. 9з) 
instead of the three Cartesian coordinates (z, y, z). These 
numbers often prove to be more suitable to the problem at 
hand. 

Let every point M be associated with a definite number 
triple (qj, 42, 9з) and, conversely, let every number triple 
correspond to a unique point M. Then the quantities 
41, 9, 9з are termed the curvilinear coordinates of the 
point 

The coordinate surfaces in a system of curvilinear 
coordinates 41, q2, ĝa ате the surfaces 


9, = С, (1) 
q2 = Ca, (2) 
Qs = Сз, (3) 


on which one of the coordinates remains constant. 

The line of intersection of two coordinate surfaces is 
called a coordinate line (azis). 

The coordinates q, and q, maintain constant values 
along the line of intersection of the coordinate surfaces 
(2) and (3); it is only the coordinate g, that varies. Simi- 
larly, on the lines of intersection of the surfaces (1) and (3) 
and (1) and (2), it is 9, and qs that vary respectively. 

Let us introduce the unit vectors ej, e,, es directed 
along tangents to the coordinate axes (91), (q2), (9з) at 
the point M in the direction of increasing variables g}, qa, 
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9з respectively (Fig. 37). Let us agree to take the unit 
vectors e,, e,, е, always in that order so that, taken 
together, they constitute a right-handed trihedral. 
The basic difference between curvilinear coordinates and 
Cartesian coordinates is this. In the Cartesian system, the 
vectors e,, е,, ез are constant for all points of space 
and are equal, respectively, to i, j, k. In any other 


Fig. 37 


system, they will, generally speaking, change their direc- 
tions when passing from one point, M, to other points. 
` Аз examples of curvilinear coordinates we consider 
cylindrical and spherical coordinates. 
1°. Cylindrical coordinates. The position of a point M 

in space is defined in cylindrical coordinates by three 
coordinates: 

91 = р, 0 < p < +оо, 

% = Ф DO Ф < 2л, (4) 

qs = 2, —oo < 2 < +æ. 


The coordinate surfaces аге: 
p = constant: circular cylinders with the z-axis; 
Ф = constant: half-planes adjoining the z-axis; 
2 = constant: planes perpendicular to the z-axis. 
The coordinate lines (or axes) are: 
(p): rays perpendicular to the z-axis and having their 
origin on that axis; 
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(Фф): circles with centre on the z-axis and lying in 
planes perpendicular to that axis; 
(z): straight lines parallel to the z-axis (Fig. 38). 
Cartesian coordinates are 
related to cylindrical coor- 
dinates by the following 
formulas: 


т = р соз ф, 
y=psing, (5) 
2 = 2. 


2°. Spherical coordinates. 
In spherical coordinates, 
the position of a point М 
іп space is defined by the 
following coordinates: 


9 = г, 0 < r« +o, 
qa = Ө, 0<6<xa, (6) 
= Ф 0О<еФ<2л. 


The coordinate surfaces 
are (Fig. 39); 
Fig. 38 г = constant: spheres 
centred at O: 
Ө = constant: circular half-angle cones with the 
z-axis; 
Ф = constant: half-planes adjoining the z-axis. 
The coordinate lines are: 
(r): rays emanating from the point O; 
(9): meridians on a sphere; 
(9): parallels on a sphere. 
Cartesian coordinates are related to spherical coordi- 
nates via the following formulas: 


т = г соз ọsin Ө, 
у =r sin ọsin Ө, (7) 
2 = г соз Ө. 
А system of curvilinear coordinates is said to be or- 
thogonal if at every point M the unit vectors е, е,, 65 
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are pairwise orthogonal. In such a system, the coordinate 
lines and the coordinate surfaces are also orthogonal. 
Systems of cylindrical and spherical coordinates are instan- 
ces of orthogonal curvilinear systems of coordinates. 


Fig. 39 


Henceforth we consider only orthogonal systems of coor- 
dinates. 
Suppose г = г (41, 42, 9з) is the radius vector of 
a point M. Then 
dr = H, dq, e, + Н, dq, e, + Нз dg; ез. (8) 


Here 
mey (a) +(ж) (m). i712 


are the Lamé coefficients of the given curvilinear system 
of coordinates. 
In cylindrical coordinates, 


Ч = 0, 9: = Ф, Яз = 2. 
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By virtue of (5) we an 


nome Ue ey E 3 =p, 


m-n- V EEICT- 


In hoien coordinates, 
Ф = г, в» = 0, Is = Ф. 
By virtue of (7) we have 


пену С 7087 


99 
Щи _ дз |? ду \2. (02 Y? 
Н. = Не = и (=) + (=) + (з) — r sin. 
The quantities 
dl, = H, dq, і = 1, 2,3 

that appear in formula (8) are differentials of the arc 
lengths of the coordinates lines, In a number of cases, this 
consideration permits of a more simple computation of 
the Lamé coefficients. For instance, in the case of cylin- 
drical coordinates (4) (see Fig. 38), the differentials of 
the arc lengths of the coordinates lines (р), (p). (z) are 

а (р) = 1:dp, whence Н, = 1; 

а (ф) =р-4ф, whence Н, = p; 

а (2) = 1.02, whence Н» = 1. 
It is just as easy to obtain expressions for the Lamé 
coefficients in the case of spherical coordinates (6). 


Sec. 24. Basic operations 
of vector analysis in curvilinear 
coordinates 


1°. Differential equations of vector lines. Suppose we 
have a vector field 


а = а; (91, 9, 9з) е + а, (9, 9, 9з) е, + аз (ду, 9, 03) ёз. 
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The vector-line equations in curvilinear coordinates 
91, 92, 9з are of the form 
Н, dq, = Hs dq = Нз дз 
а! (91, 9з. 9з) а (41, Gar Ҹә) as (91, qa» 93) ` 
In particular, in cylindrical coordinates (4, = р, 9 = Фф, 
Qs = 2): 
—e — R41 X __. (1) 
a, (0, Ф, 2) а, (0, Ф, 2) аз (р, Ф, 2) ' 
in spherical coordinates (4, = г, 9, = 9, 9з = Ф): 
аг = г 40 _ rSin0dg 
в, (г, 8, $) а, (г, 0, $) ^ аз(г, Ө, Ф) * 
. Example 1. A vector field is given in cylindrical coor- 


dinates 
a (M) = e, + ges. 
Find the vector lines of the field. 
Solution. It is given that a, = 1, а, = q, ag = 0. By 
virtue of formula (1) we have 


whence 


which are Archimedean spirals lying in planes parallel 
to the zy-plane. 

2°. The gradient in orthogonal coordinates. Suppose we 
have a scalar field 


и = и (41, Qs, 03). 
Тһеп 


1 д д ди 
gradu ез, 1+ H, э, t H; э, © 
In particular, in cylindrical coordinates (9 = p, 9, = 
= Ф, qs = 2): 
1 ди 


тай и = эе e+ ete e; (2) 


in spherical coordinates (q, = r, 4, = Ө, qs = 9): 


1 ди ди 
gradu A e, T 0 eot rang во (3) 
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Example 2. Compute the gradient of the following sca- 
lar field specified in cylindrical coordinates (р, ф, 2): 


и = р + 2 соѕ Ф. 
Solution. Using formula (2), we obtain 
1 Я 
gradu=1 "€p— -5 2 Sin Q- e, + COS фе. 

Example 3. Find the gradient of the following scalar 
field given in spherical coordinates (r, 9, q): 
sin 6 

г 


u-r4 — зїп cos ф. 


Solution. Using formula (3), we have 
in Ө егі і 
gradu= (1-225) e + = (-- — соз ) eo + 229 ey. 
3°. The curl in orthogonal coordinates. Suppose 
а=а, (9, 9, 93) € + аз (91, 9, 9) е + аз (ду, 9, 9) €s- 
Then di | 


1 04 4 
TA, © HH, © Н.Н, °з 
curla=| 9 2 EU . 
991 дая 09 


a,H, aH: азИз 


In particular, in cylindrical coordinates (4, = p, 9, = 
= Ф, 9s = 2): 


1 1 
pe es pe 

curla=] 2 ә 2 |; (4) 
др дф д2 


а ра аз 
in spherical coordinates (9 = г, 4» = 8, дз = Фф): 


P MET ы 37 
risin8 ^" rsino 9 r ^9 
cula=| 2 2 a |. w 
ar 80 — 39 


а, ra, rsin®@-a, 
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Example 4. Compute the curl of the following vector 
field specified in cylindrical coordinates: 


eo — рге,. 


Solution. Taking advantage of formula (4), we get 


1 1 
ry ер eo ri ez 
eula-|? ә ә |--Le(0—0)—e,(—2—0) 
р õp д: e 


sinp соѕф — pz 


tle (0 — cos ф) = ze; — E» 


e,. 
4°. The divergence in orthogonal coordinates. Suppose we 
have a vector field 


a—4, (9, 9, 9з) €i + а» (Чу, 9, @з) ез + аз (9, 9, 93) €s 
Then 


div a= 


1 O(a,HsHs) , д(азН.Нз) | ô (a4H,H,) 
H,H,H, [ да\ + ддз LE rA ддз 1. 
In particular, in cylindrical coordinates (q, = р, qa = 9, 
4з = 2): 
1 4 (pa,) 1 да, | даз. 
Р ur tas * 
in spherical coordinates (4, = г, qa = 9, qs = Ф): 
eae | IOS. 1 д(зїпӨ-а„) 1 das 
div a = " or * + rsinü 00 T rsin$ дф ' (5) 
Example 5. ae that the vector field 
= ieee ine e 


ec 


is solenoidal. 
Solution. Using formula (5), we have 


1 д 2cos 0 1 ant 
div a— тт эг (7 гЗ EET я 90 (sin 877 ) +9 
ыш [жш ыы 0 


wherever г =£0, which means that the field of the vector 
a is solenoidal at all points with the exception of r = 0. 
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267. Find the equations of the following vector fields: 
(a) а=е,+--е,--е„ 


(b) a= pep + 9e, + 2е,; 

жоон e, 4- аз 
Find the gradients of the following scalar fields. 
(a) In cylindrical coordinates: 

268. и = р? + 2р cos ọ — е sin $. 

269. и = p cos ф + z sin? ф — 3°, 
(b) In spherical coordinates: 

270. и = г? cos Ө. 

271. и = 3r? sin Ө + е" cos ф — г. 

272. u2u ке» ‚ p-constant. 
Compute the divergence of the following vectors. 
(a) In cylindrical coordinates: 

273. а = pe, + z sin q-e, + e? cos z-e,. 

274. а = arctan p-e, + 2e, — 22ее,. 
(b) In spherical coordinates: 


(с) a= ee, о, = constant. 


275. a = r'e, — 2 cos? p -eo + ат êv- 


Compute the curl of the following vector fields: 
276. а = (2r + a cos ф) e, — а sin 0.eg + r cos Ө-е„, 
о = constant. 

277. а = ге, + 2 cos 0. — pey. 


278. a = cos q. e; — c €, + pre,. 
279. Show that the vector field 
2 0 in 0 
a= = e, +e, 


is a potential field. 
280. Show that the vector field 


a =f (г) e, 


where f is any differentiable function, is a potential field. 
5°. Computing the flux in curvilinear coordinates. Let 
5 be a part of the coordinate surface q, = C, where С = 
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= constant, bounded by the coordinate lines 


Ф = ©, 9 = а (a, « a4); 
ds = В,, ga = В, (В, < Вз). 
Then the flux of the vector 


а = а, (91, 9, 9з) €i + а, (9, 42, 93) ез + аз (9, 9, 9з) ез 
through the surface S in the direction of the vector е, is 
computed from the formula 

а, Bs 


n= { ) a, (C, gar 9) Ha (C, Ч» 43) Ha (C, ga, 93) dgs dga. 


i (6) 


The calculation is similar for the flux of a vector 
through a part of the surface q, — C or through a part of 
the surface qa = C, where С = constant. 

Example 6. Compute the flux of the vector field, speci- 
fied in cylindrical coordinates, 

а = pe, + ze, 
through the outer part of the lateral surface of the cylin- 
der p = 1, bounded by the planes 2 = 0 and z = 1. 

Solution. The cylinder is the coordinate surface p = 

= C = constant and so the desired flux 


2n 1 
n- f j C? dz dg = 2nC?, 
0 0 
whence for the surface p — 1 we obtain 
П = 2n. 


Example 7. Find the flux of the vector field, specified 

in spherical coordinates, 
а = где, + ге Beg 

through the outer side of an upper hemisphere S of ra- 
dius R with centre at the coordinate origin. 

Solution. The hemisphere $ is part of the coordinate 
surface г = constant, namely, г = R. On the surface 5 
we have 


qı=r=R; 9 = 0, 0<0< + ; 
з=Ф, 0<Ф<2л. 
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Taking into account that in spherical coordinates 
H, = Н, = 1, Н, = Н = г, Н; = Но =rsin8, 
we find, via (6), that 
n/2 2л л/2 
П = j ae | НӨ sin Ө dg = 2л А | O sin Ө 29 = 2x R*, 
e 0 0 
Compute the flux of the vector field, specified in 
cylindrical coordinates, through the given surface S. 
281. а = pe, — cos q:e, + 2е,; 5 is a closed surface 
formed by the cylinder p = 2 and by the planes z = 0 
and 2 = 2. 
282. а = pep + pqe, — 2ze,; S is a closed surface 
formed by the cylinder p —1, the half-planes ф = 0 and 
ф = л/2, and by the planes 2 = —1 and z = 1. 
283. Find the flux of the vector field a — (1/r?) e, through 
a sphere of radius R with centre at the coordinate origin. 
284. Find the flux of the vector, specified in spherical 
coordinates, 
а = ге, + r sin 0-e9 — 3rq sin 0-e, 


through an upper. hemisphere of radius R. 
285. Find the flux of the vector, specified in spherical 
coordinates, 

а = r'e, + А? cos p'ep 
through the sphere r — R. 
286. Find the flux of the vector, specified in spherical 
coordinates, 

а = re, — r sin 0-е, 


through a semicircle of radius R located in the half-plane 
ф = л/4 (the flux is taken in the direction of the vector 


ey). 
287, Find the flux of the vector, specified in spherical 
coordinates, 


a=rsin $ eor sin Ө cos p.e, 
through the outer side of part of the half-angle cone 


y 32 = z’ + у, bounded from above by the plane z = 
=V30<z< ИЗ). 
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6°. Finding the potential in curvilinear coordinates. Giv- 
en, in curvilinear coordinates q,, 93, qs, the vector field 


а (М) = a, (9, 9», 9з) € + а» (41, 9, 05) ез + 
+ аз (9, 9, 4з) ез; 
this is a potential field in some region © over which 
the variables q}, 42, 9з range, that is, curl a = 0 in ©. 

To find the potential и = u (ду, 42, 9з) of this field, 
write the equation a ui = ge u (M) as 

1 du 1 ди 

аце, + але, + азе; = F, On ^ e, Zz 7 Be e, ur ез. 

From this it follows that 

ди ди ди 

en ‘да, = Hs ‘да; = @3Нз. (7) 
This is a system of partial differential equations whose 
integration yields the desired potential u=u (41, ga» 9з) + 
+ C, where C is an arbitrary constant. 

The system (7) of differential equations is solved in 
the same way as in finding the potential in Cartesian 
coordinates 

The system (7) of differential equations is of the fol- 
lowing form: 

(1) In cylindrical coordinates (4, = р, 4, = Ф, 93 = 2), 

и и П 
TE =p, Fe = Pg, Sr = Gy. (7) 
(2) In spherical coordinates (4, = г, 4» = 0, дз = Фф), 
ди 0, Bara, S=raind- (7°) 
ar @т, "30 9 "ag "09: 

Example 8. Find the potential of the following vector 

field specified in cylindrical coordinates: 


= a,H,, 


a= (E соз 9) ер — зір P-e t 7-7 iere 
Solution. By formula (4) we find 
1 1 
"p e? е rA 
д a д 
curla= p ЭФ os =0 (р2> 0). 
m: 2 А Іар 


———-Feosp —psin@ +s 
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This is a potential field. The desired potential и = 
= и (p, Ф, 2) is the solution of the following system of 
differential equations: 


ди — arctan 2 

e + cos Фф, | 
E sin 

дф I р Ф, 

Ou p 

ôz 14-22 


From the first equation, integrating with respect to p 
we obtain 


и = In p-arctan 2 + p cos Фф + С (Ф, 2). (8) 
Differentiating (8) with respect to s we get 
ди 
99 ^ —psin e +E 


and since ди/дф = —р sin Ф, it follows that дС/дф = 0, 
that is, С = С, (2). Thus 
и = In p-arctan z + p cos o + С, (2), 
whence 
x 1 , 
= Tor 60. 
By virtue of the а equation of the system we have 
1 1 " 
TES pps +С @ 
or С; (2) = 0, whence С, (2) = С = constant. 
То summarize, the potential of the given field is 
и (p, 9, 2) = Inp-arctanz + p cos ọ + С. 
In the following problems, verify that the vector fields 
given in cylindrical coordinates are potential fields and 
find their potentials. 


288. a= e+ e+e. 
289. a = ре,---#-е„--ге,. 
290. а = pze,+ Ze + рфе,. 


291. а = e?sing-e,+ E eP cos Q- e, + 22е,. 


Sec. 24] Basic operations of vector analysis 165 


292. а = 9 cosz-e, + cos ze, + pp sin 2.е,. 
Example 9. Find the potential of the following vector 
field given in spherical coordinates: 


a= 5. e99e, + Sar e99e, + mr peeg. 
Solution. Using (4), we find that 
е, ree r sin Ө.еф 
д д д 
сапа = 95| or 28 дф |=0. 


= e olnr-e Ө Inr.e9e 


This is a potential field in the region where г > 0,0 nn 
(n = 0, +1,...). 

The system (7) of differential equations for finding the 
potential и = u(r, Ө, q) із of the form 


2t — pe Inr, (9) 


Integrating the first equation of system (9), we obtain 
и = e lnr + С (Фф, 9). (10) 

Differentiating (10) with respect to Ө and taking into 

account the second equation of the system, we have 
pe 1n r = pee Inr I. 


That is, 0C/00 = 0, whence C (Ф, Ө) = С, (Ф) and 


therefore 
и =e ln r + С, ($). (11) 


Differentiating (11) with respect to q and taking into 
account the third equation of system (9), we obtain 
099 In г = Be Inr+ C, (9) 
or С,(ф)==0, whence С, (ф) = С = constant. The 

desired potential is 
u (r, 9, Фф) = е8 Inr + C. 
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Establish the potentiality of the following vector 
fields, specified in spherical coordinates, and find their 
potentials. 


293. a = Өе, + ee. 


1 1 
294. а = 2re, + бе b ee. 
= Ф 0 
295. а = > фе, + ang eo + — ёо. 
296. а = cos фѕіп 0-е, + cos pcos Ө .ев — sin 9-e,. 


A 2 
297. а =e" sin 0-е, +- e" cos 0-eo + Ga ae es. 


7°. Computing the line integral and the circulation of 
a vector field in curvilinear coordinates. Suppose a vector 
field 
a (M) = a, (ау, ds. 9s) е + аз (41, 92, 43) €z 
+ аз (9, 9, 9з) ез 


is defined and is continuous in a region Q over which the 
orthogonal curvilinear coordinates 41, 92, qa range. 

As we know [see Sec. 23, (8)], the differential dr of 
the radius vector г of any point M (qis 42, 93) EQ is 
equal to 

dr == H, dq, e, + Н dq, ey + Нз dq; es. 


Therefore the line integral of the vector a (M) over an 
oriented smooth or piecewise smooth curve L с © is 


| ка, а= | aH, dg, + аан. ад. + азН44. (12) 
L L 


In particular, for the cylindrical coordinates q, = p» 
92 = Ф, q3 = 2 we have 


а = а, (p, P, 2) eo + ay (p, Ф, 2) е + а. (p, P, 2) ez» 
dr = dp-e, + раф-е, + dz-e,, 
and therefore 


| в, а) = | а, dp + asp do + a, dz. (13) 
L L 


Sec. 24) Basic operations of vector analysis 167 


For the spherical coordinates q, = г, qa = 0, дз = Фф ме 
have 
а = а, (г, Ө, 9) е, + ав (г, Ө, Ф) eo + а (7, 0, 9) e,, 
dr = dr.e, + г 49-е + r sin 0 dq ‘ey, 
and consequently 


| (a, dr) = | a, dr + тав 00 + ra, sin 0 dp. (14) 


The circulation C of the vector field a (M) is computed 
in the curvilinear coordinates g}, 4», qa via formula (12) 
in the general case; for cylindrical or spherical coordi- 
nates it is computed from (13) and (14) respectively. 
Example 10. Compute the line integral in the vector 
field, given in cylindrical coordinates, 
а = 4p sin p-ep + zee, + (p + 9) e, 


along the straight line 


= л . 
L: Mig , 
2=0, 


from the point О (0, л/4, 0) to the point A (1, x/4, 0). 
Solution. In the given example, 
a, = 4р sing, а» = ze?, a, = р + Ф. 
By formula (13), the desired line integral is 


f (a, dr) — | 4p sin фар + pze? dp + (p + Ф) dz. 
L L 


On the straight line L we have 
p=}, dp=0; 2=0, dz=0; Oxp«. 


Therefore 
1 


f (a, dr) = | 2V2pdp=V2 f 2pdp— V2. 
L L 0 
Example 11. Compute the line integral in the vector 
field, given in spherical coordinates, 


а = e sin 0.e, + 30° sin p-eẹ + robe, 
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along the line 


in the direction from the point М, (1, 0, x/2) to the 
point M, (1, л/2, x/2) (Fig. 40). 


Fig. 40 


Solution. The curve L is an arc of a circle with centre 
at the coordinate origin and radius r — 1 located in the 
yz-plane. The coordinates of this vector are 


a, = e зіп Ө, ae = 30° sin ф, a, = гфӨ. 
By virtue of (14) the line integral is of the form 


f (a, dr) = f е” sin Ө dr + 302г sin ф 20 + г2ф0 sin 0 dọ. 
L L 


Taking into account that the following conditions hold 
on L, 
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we obtain 
1/2 $ 
л 
| (а, dr) = | 382 40 = ( 30246 = T. 
L L 0 


Example 12. Compute the circulation of the vector 
field, given in cylindrical coordinates, 
а =p sin ф-е„ + р 2е + pe, 
over the curve 
= Sin Ф, 
L: {Pro 0<е<л, 
directly and via the Stokes theorem. 
Solution. The coordinates of the vector are 
ар =psing, a, = pz, a, = р?. 
The contour L is a closed curve located in the plane z = 0 
(Fig. 41). 
(1) Direct calculation of the circulation. 


Fig. 41 


Substituting the coordinates of the vector into (13), 
we obtain 


C= $ p sin фар + p?z de + p? dz. 
L 


On the curve L we have 
z= 0, dz = 0; р = зіп ф, dp = cos q dg, 
О<Ф< л: 
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Therefore the desired circulation is equal to 
л 


C =Q psin e dp | sin? o cos o do =0. 
L 0 


(2) Computing the circulation via the Stokes theorem. 
By the Stokes theorem, the desired circulation is equal 
to 


C= $ (a, dr) = | | (curla, n°) dS, 
L S 


where $ is the surface spanning the contour L. 
We find the curl of the given field: 
е рер е, 
1 д à д 
curla = 5- Эр 3e 3: = — pe, — Зр?еф 
psing р22 p° 
+ (2z — cos 9) e,. 
At points where p = 0 we redefine the value of curl a 
with respect to continuity, setting 
curl a (0, p, 2) = (22 — cos 9) e.. 
Thus, curl a is defined throughout three-dimensional 
space. Since the curve L lies in the plane z — O, for 
the surface S spanning this curve we take that portion 
of the plane z = 0 that is bounded by the curve L. Then 
we can take the unit vector e, for the unit vector of the 
normal n? to the surface S, that is, n? = e,. We find 
the scalar product: 
(curl a, n?) = (—pe, — Зр?е + (22 — cos 9) ez, е,) 
= 22 — cos ф 
because by virtue of the orthonormality of the basis 
ер, €y, е, we have 
(ep, e) = (es, е.) = 0, (е., e) = 1. 
The desired circulation is 


C= Jj (22 — cos g) dS. 
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Taking into account that 2 = 0 on 5 and the element of 
area 4$ of the coordinate surface 2 = 0 is equal to 


dS = p dp dg, 
we finally get 


C-— | cos p d$ = — № cos фр dp dp 
л sing л 


= — | cos фаф j фас = -4 f 3112 pcos 9 dọ = 0. 
0 0 0 
Example 13. Compute the circulation of the vector, 
specified in spherical coordinates, 
а = re, + (В + г) sin 8.e, 
around the circle 


in the direction of increasing values of the angle ф, 
directly and via the Stokes theorem. 
Solution. In this example, 


аа =r, а = 0, a, = (В + г) зіп Ө. 


(1) Direct computation of the circulation. 
By formula (14) the desired circulation is equal to 


C - brár- (R^-r)sinr sin Odp= 
L 
=rdr+r(R+r) sin? Ө dọ. 
L 


On the given circle L, the centre of which lies at the 
coordinate origin, we have 
г= В, dr=0; 0-7; 0<ф< 2л, 


and, consequently, 


C -2Rm6 ау =? | do indt 
L 0 
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(2) Computation of the circulation via the Stokes 


theorem. 
The desired circulation is, by the Stokes theorem, equal 
to 


C= $ (a, dr) = | | (curl a, n?) 45, 
L 8 


where S is the surface spanning the circle L. 
We find the curl of the given vector: 


e, reo rsin@-e, 
а= КД Кд EN 
~ risin8|ór 00 дф 
г 0 (Ar+r?)sin?6 


=2 (Я +r) cos @-e,—+ (R+ 27) зіп 0-e9, 


For the surface S spanning the circle L we take, for 
example, the upper hemisphere of radius А: г = В, 
0 < 6 < 2/2, 0< q < 2n. The unit vector of the normal 
n? to the outer side of the hemisphere S is directed along 
the vector e,, and so we take n? —e,. We find the scalar 
product 
R+2r 


r 


2G t7) cos 0.e, — 
T r 


(curl a, n°) = ( sin 0-е, e.) 


= AREN cos 0 
since (e,, e,) = 1, (e,, eo) = 0. 

Taking into account that г = R on the surface S, we 
obtain the following expression for the desired flux: 


n= || HREN cosdds = 4 | 08645. 
"s S 
In spherical coordinates, the element of area dS of the 


coordinate surface r = R, that is, the hemisphere S, is 
equal to 


45 = В? sin Ө 29 dq 
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n=4{ | cos 0 А? sin Ө 20 dp 


8 
1/2 2л 
а 4m cos 0 sin Ө ae аф = 4x R2, 
0 0 


By the Stokes theorem we get С = 4nR?. 

For the surface S spanning the circle L we could take 
the lower hemisphere, the unit vector of the normal of 
which, n? = —e,, and the result would be the same: 
С = 4n R*. 

Note that it is not desirable to take the circle bounded 
by L for the surface S spanning the circle L because the 
circle has a point r — O (the centre of the circle) in which 
the curl of the given vector has a discontinuity. 

Compute the line integral along the given curves L in 
the following vector fields specified in cylindrical coordi- 
nates. 

298. a = ze, + ppe, + cos g-e,; L is a segment of the 
straight line: (p = a, p = 0, 0 xz « 1). 

299. а = pe, + 2pqe, + zez; L is the semicircle: (p = 1, 
z=0,0< 9 x n). 

300. а = е? cos g-e, + p sin g-e, + ре,; L із a turn 
of the helical curve: {р = R, 2 = Фф, 0x: 9 « 2x). 

Compute the line integral over the given curve L 
in the following vector fields given in spherical coordi- 
nates. 

301. a = e" cos 6-e, + 20 cos ф-еь + ge,; L is the semi- 
circle: {г = 1, 9 20,0 0 < n). 


302. а = 473 tan + e, + Opes + cos? p-ey; Lis a segment 
of the straight line: (9-3: d=, 0<r<i}. 
303. a= sin? 0-e,+sin 0-eg 4-rq0es; L is a segment of 


the straight line: (9-5. к=. ixl). 


Compute the circulation of the following vector fields, 
specified in cylindrical coordinates, about the given con- 
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tours, directly and with the aid of the Stokes theorem. 
304. а = ze, + рте, + ре,; L is the circle: {p = 1, 
z = 0}. 
305. a = p sin q-e, — p?ze, + р?е,; Г is the circle: 
{р = R, z = R}. 
306. a = 2 cos p-ep + рез + ф?е,; L is the loop: {р = 
= sin g, 2 = 1). 

Compute the circulation of the following vectors, given 
in spherical coordinates, along the given contours L, 
directly and with the aid of the Stokes theorem. 


307. а = ге, +r sin 6-e,; L is the circle: {r =1, 0 = 
= 7} 

=F}. 

308. а = r sin 0.е, + Өеёеө; L is the loop: {r = sin Ф, 
0 = 5, 0 < Ф< л} à 

309. a—rq0e,; L is а contour bounded by the semi- 
circle: {r=R, ф= 7, 0<9<л) and its vertical 


diameter {o = i 6 = 0} Я 


Sec. 25. The Laplace operator 
in orthogonal coordinates 


If u = и (9, 4», 9з) is a scalar function, then 


1 ди ди ди 
rad u =-~- —— ез. 1 
g H, àq €, etr T; E е? +- 7, gg. 3 ( ) 


If 
а=а, (41, 92, 9з) е + а, (41, 9, 93) ез + аз (91, 4», 93) ез, 
then 


2 (a4H5Hs) + 2 (агНзН\) 
1 92 


div a= LEA 
H,H,H, | да 


+з НН]. Q) 
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Using formulas (1) and (2), we obtain the following expres- 
sion for the Laplace operator Au: 


i 1 Н.Н, ди 
Au = div grad и= и lee 2з 
Ё НаНзНз м | Ну дау 
+2 (НН, ди 2 IT = 
даз Н, ôq 943 Н, ддз 


In cylindrical coordinates, 


1 д ди д 1 ди 
А s (096) tay o ag) +3; (07) | 
1 ôu , ди 
) ra 023 * 
Ín spherical coordinates, 
4 a ee, 
Au— Ass [2 (resino 2) +3 (sin ө) 
д 1 ди 1 д ди 
tay ay (mar =) === (" a) 
1 Я ди \ , 1 д?и 
+ ғаш я (sin ө29) T asino 99: 

Example. Find all the solutions of the Laplace equation 
Au = 0 that depend solely on the distance r. 

Solution. Writing the Laplace equation in spherical 
coordinates and taking into account the spherical sym- 
metry of the solution (it must not depend on Ө or $), 
we have 


1 8 д 
Аи= у (rt)=0 (=u), 
whence 
ди 
rz = Са 
so that 
um A С, 


where C, and C, are constants. 
310. Given: a scalar field и = и (M) in cylindrical coor- 
dinates 
и (p, Ф, 2) = p*o + zig? — poz. 
Find Au. 
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311. Given: a scalar field u = и (M) in spherical coor- 


dinates 
и (г, Ө, 9) = r*q0 + Pg? + 9 + 07. 
Find Au. 
312. Are the following functions harmonic? 
(1) и = p? cos 29. 
(2) и = rcos 28. 
313. Find all possible harmonic functions that 
(a) depend on Ө alone, 
(b) depend on q alone 
(in the spherical system of coordinates). 
314. Find all solutions of the Poisson equation 


Au = rh 


in the spherical system of coordinates, provided u = и (г). 


ANSWERS 


| z22, 
1. (a) The half-line | у= — 2,0220, 2<0 is traversed 


twice when— оо < t< + oo. (b) и А "s 


U(—1, +œ) the point r(t) = 1+ 


war en j 
twice traverses the half-line z -- y = 1, sot, ,<->. 
з 23 
ка| (d) у=, =. 
(c) z=1;)° 
(e) z?+-y2+22=1, r—y=0. 
7. i--k. 8. i+k. 9. —j+ k 
10. —i+k. 11. ei—j+2k. 12. No. 14. No. 
dr . 2 д?г ү, dir 
17. (a) 2 (3 г) +(r, а); (с) E ar |. 


21. Circles lying in planes perpendicular to the vector a. 

22. The hodograph of velocity is a helical curve: z = 

= a cos t, у = a sin ї, 2 = 2bt; the hodograph of accele- 

ration is a circle: z = —a sin t, у = a cos t, 2 = 2b. 
da da du , да _ d?a (du da d*u 

26. а= 97 ae ae = aur (ae) tae a 


28. (2—1) e'i-- + (t£ — sin 2t) j — arctan t- К е. 
29. jn ЧР 4 4- + ej -- sin t-k4- c. 


30. eint. 4 — 3. sin 2.jtitk+e. 
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E] t 
31. $ i+ (£cost — sin t)j+ 125 К+е. 
2 


32. zitak. 33. (1—e- 2) i+ (e!/2—1)j+(e—1)k. 


, | 2 
34. —In2-j+k. 35. 274 --лј + л2К. 36. R= iu ; 


37. R= 2 |11(14-9/2)2/, 38. R=6. 39. А =--ал. 


40. В = 2а соѕһ2 1. 41. х +у = 0. 


5 1 1 
42. z—-y—V/ 2z=0. 43. T^3- 44. 


d scc cr. 
T  2acosh?t ^ 


45. Рича -C, a family of triaxial ellipsoids. 


46. z? + у? —z = C, a family of paraboloids. 

47. z? + y? = Cz, a family of paraboloids. 

48. 2y? + 922 — C, a family of elliptical cylinders. 
49. z + 2y — z = C, a family of parallel planes. 

50. A family of planes resulting from the sheaf of planes 
az + а,у + аз = C (bız + bay + 6:2) passing through 
the straight line 


ат + ау + аз2 = 0, | 
biz + Бу + bz = 0. 

via the elimination of the straight line itself. Here, a;, а», 

аз are the coordinates of the vector a; b,, bg, bs are the 


coordinates of the vector b. 
51. 22 + y? + 2 = C*, a family of concentric spheres. 


т y2 

52. (a, b, г) =С or |a, a, аз| =C, a family of parallel 
b, b, b; 

planes. 

53. 22 — y = C, a family of parallel straight lines. 

54. у = Cz, С> 0, z #0, a family of rays. 

95. у? = Cz is a family of parabolas with the vertex 

O (0, 0) deleted. 

56. z? — у? = C, a family of hyperbolas. 


57. 
98. 
62. 
67. 


70. 
73. 


75. 


81. 


87. 
88. 


90. 
91. 
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у = —х1п C — С, C > 0, a family of straight lines. 


У15 59 —Y? 60, УЗ o. бол, 


3 y3. 63. 1. 64. 0. 65. 2У3 (y3+3). 66. 0. 


—2. 68. м. 69. $ (i--j— k). 


К. 71. ф=л. 72. 9 —0. 
Фф = 0. 74. у = —z + 2пл, n = 0, +1, +2,.... 


22+ у2+22=1. 78.—-. 79. а. 80. a (b, r)+b (а, r). 


2|а|?г — 2 (a, r)a. 86. LUN 


дг тг’ 
a Cb. a 
и cos (г, 1) . и _ 
A^ xi э=0 for r | l. 
ди 1 ди 
аг я. 89. a= 


ди _ (grad u, grad v), ди 0. 
ЭТ” Татай vj iF =0 if grad u | grad v. 


(a) 1 in the direction of the y-axis; (b) 3 in the 


direction of the vector a= —i—2j+ 2k. 


92. 


94. 


у=е,т; z=e,2. 93. у= 2+ С; z= zd C. 
1 


1 1 1 1 
22+ y2+ = С}, | 95. soya тъз 
+ у+=бС;. 


96. 12=Сцу, 2=0,. 97. 2=Cyz, у=С.. 98. zy—C, z=}. 


99. 


101. 
103. 


z=C,, 9у%—%=С,. 100. +-у=бь 2=C. 
y2+22=C,, х=С.. 102. х= Слу, z—Cs. 
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where Бо, Bos, bog are cooordinates of the vector by. 
104. П = —3. 105. П = xR*y. 106. П = x Rh. 


107. П = 4x Ref (В). 108. II— 5-. 109. НЕХ. 
110. П= лл. 111. Пл. 

81 л 
112. п= S x, из. T=, 


144. П=0. 115. П= 2, 116. n=. 417. IT— 40k, 


2 

118. (a) M=>—2; (b) П= — 7; © П=л. 119. П=0. 
120. II—6zR. 121. П=0. 

У? 
122. П-л. 123. П=0. 124. II- 2л (1— 7). 
125. П=8 f. 
126. П-И? a. 127. П=45л. 128. п= 22g. 
129. I1— 0. 


130. П=л. 131.  (r) = С-. 132. 7r*. 133. 0. 134. 0. 


135. їр (2) = С — 2, С = constant. 136. П = 4x Аз. 
137. div E = 0 (г 520). 


143. 16x. 144. nH3*. 145. X n. 146. 0. 147. +. 
148. 4n. 149. Dr. 150. 2 д. 151. 2R3. 


152. S. s. 153. —1. 154. — 


155. Solenoidal field. 
156. Nonsolenoidal field. 
197. м field. 


159. ф(г)=—, г=20, C —constant. 
161. ал 162. In, 163. ++, 164. 0. 
1 
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4 4 14 2 
166. — я. R. 167. ^L. 168. (а) — =; (b). 169. 0. 


170. >. 171. 3V3. 172. 
173. — ла?. 174. 1. 

175. — 2л. 176. — ЭЖ“. 177.4. 

179. — 2 (zi 4- zj + yk). 

180. 3 (12—22) j. 181. (z--y) К. 191. o — curly = —4j, 


193. f(z, 2) = 12 - zx -z4 C, C —constant. 
195. 4л. 196. —4л. 197. +. 198. —2n. 199. 18. 
200. 729л. 201. 0. 202. —V2n. 


203. 20ла2. Hint: v=[o, г]. 
204. и. = 1. 205. р, = 3. 206. Dependent. 


207. Independent. 208. Dependent. 209. — 1. 210. $. 
2 л 1 
2и. 0. 212. 2. 213. 1. 244. +. 


1 
35 ° 


216. >. Hint: Supplement the path of integration L 


with the line segment OA of the z-axis. 

217. No. 218. Yes. 219. No. 220. Yes. 221. No. 222. No. 
223. Yes. 226. ф = z?yz. 227. р = z + хуг. 228. 9 = 
= 22у — у? + zz. 229. ф = | т + y +z]. 230. ф = 
= arctan (zyz). 231. ф = г. 232. 9 = Inr. 


233. 9-4. 234. o — az -- By -F yz 4- C, C = constant. 


235. ф = zy + yz + 22 + C. 236. ф = zy - e + С. 
237. ф = e* siny + 2 + С. 247. (a) Yes (b) No, 
(c) Yes. 

249. и = Сл + C,. 

250. и = Az? + Bry + Ay?, where А and В are arbi- 
trary. 


тау +С if n1, 


251. T 
z ln |z| - Coz 4- C; if n—1(23*0). 
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252. T= n. 253. [= —2. 254. [= i nR. 


3 
255. І = 0. 
256. b = zj + (y — х) К. 257. b = (у: — 2zz) К. 
258. b = (e* — ze") j. 259. b = 3212} + (2j? — 622) К. 
260. b = —= (z + у?) ј + (2° + у) К. 


261. b = — (xz? + у2ех°) j — 2zyzk. 
262. b= (— zi - yk). 

263. b = —8yzi + zzj + 7zyk. 

264. b = 2zy*zi — 3z?yzj + дук. 
265. —Lsin 21.11 sin zz-k. 


266. b— AB z—k. 267. (a) р=ф-+ Сү, p=z+Cy; 


(b p= — ug. р Са; (0) ф= C» г= C, sin. 


268. grad и= 2 (р + соз ф)е, — (2sin 9+ te cos o) es 
—e’ sin ф:е,. 
269. grad и = (соѕф — 3? 1n 3)e,+ (= sin 2ф— зіп ә) е, 
+ sin? ф.е,. 
270. grad и = 2г соз Ө .е, — г sin Ө.еө. 
271. grad и = (бг sin Ө + e" cos p — 1) e, + 3r cos 0 -eo 


ет sin @ 
rsinð T 


272. grad u=—p (22е è+ sin’ e). 


273. div E 
ivaz-X. ОЕ. t d 
274. div a— р arctan p+ T+? (224+ 22) e. 
É 2 
275. div a=4r—Ž cosp cot 0-- cr rT: 


_ соз 20 5їпф asin Ө 
276. сипа =- È ‚— (2co50.- 557 5- rein ITT) ео Cy. 


277. 
278. 
281. 


285. 
287. 
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curl а= — cot Ө. е, + T. ео + 20088 о ер. 


curl a= —2pe, + ——- sin oe 

24л. 282. 5л. 283. 4л. 284. 2 лн. 

án Ң*. 286. -$ ns. 

48. Hint: Write the equations of the surfaces in 


spherical coordinates. 
288. u =p 4-9 + 2 + С. 


289. 


и=+ (p^ 2)4C. 290. и=рф2 + С. 


‚ и = е? sin ф + z? + C. 292. и = рф соз 2 + С. 


293. u = гд + C. 294. и = г + ф + 9 + С. 
в p? + 0?) + С. 296. и = r cos o sin 0 + C. 
297. u = p rU nr 

298. 1. 299. n*. 300. 2лВ. 301. n*. 302. 1. 

303. eia 

304. 0. 305. —2nR*. 306. л. 


314. 


. (1) u (8) C, In 


. л. 308. 0. 309. 0. 
. Аи = 4g — $ om +298, 


. Ди = = 698 + 1279 + = Z p. cot 6+ 29 cot 0. S 
. (4) Yes, (2) No. 


а (2) «(9) — C9 - C. 
гп 

went a rcc 

и (г) = In r4. — C, n=—1, (г 5 0) 


| рас, п= — 2. 


APPENDIX | 
BASIC OPERATIONS OF VECTOR 
ANALYSIS IN ORTHOGONAL 


CURVILINEAR COORDINATES 


1. The scalar field is given in orthogonal curvilinear 
coordinates, и = и (41, 9, 43). Then we have 
1 du 1 ди 1 ди 
ad u = 4— — ее — 
gradu Н, да, (++, дд» e Нз ддз x 
The Laplace operator is 
mold [2 (++ u) 0 (HH, ôu 
aw HHH, E ( Н, н ы ( Н, 095 
д | A,H, ди 
+, ( H, 2%, } | 


Special cases: (a) The scalar field is given in cylindrical 
coordinates, и = и (p, ф, 2). Then ^ "m 
ди 1 ди 

grad и = + p 9 e + 2% are 
The Laplace operator is 
ey fou ди. 
Au а (р жук ap 22 
(b) The scalar field is pom in spherical coordinates, 


u = u (г, Ө, $). Then we have 
gradu=e,++ : A ө +1 ee t 

The Laplace operator is 
Ty (7-9, ) + mae ae (99 sin) 


ди = ar 
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2. The vector field is given in orthogonal curvilinear 
coordinates: 


а = а, (d, 93, 93) €x + аз (41, 9, 9з) е» 


+ аз (41, д», qs) ез. 
Then we have 


А 1 д (a, H4H4) д (аз НН) | ô (аз Н.Н) Нз) 
div a= H,H,H; [ 941 + 94з 943 | 
Н.Н, ! H,H, 2 Н.Н, 3 
curl а = 2 2 2. 
дау 94 94з 


aH, а,Н, asl; 


Special cases: (a) The vector field is given in cylindrical 
Coordinates, 


а = а: (р, Ф, 2) е, + G4 (p, Ф, 2) е + аз (р, Ф, 2)е.. 
Then we have 


1 ô (pa) 1 даз Dr 
V3.5 т , 
diva tp ap T A. 
1 
р ё е — е; 
curl а= |2. oy 
EE" ap 02 
a, Pa, аз 


(b) The vector field is given in spherical coordinates: 


a = a, (г, 0, фе, + a, (г, Ө, Ф) ев + аз (r, 0, Ф) ep. 
Then we have 


_ 1 @(a,r?) 4 д (ag sin Ө) 1 даз 
div ам 3— trane a0 ^ Tsino o9" 
1 1 ё á 
resin © rsinà ^? г? 
айса: | © 2 EU 
4 or 00 дф 


а; га, agr sin Ө 
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AREA ELEMENTS 
OF COORDINATE SURFACES 


Coordl- 
nates 


Coordinate surfaces 


Area elements 


General | ^*^ С = constant |d$,—H (С, qs. 9з) Н» (С, 9з, 43) дз баз 
qs = C = constant |45.=Н, (ду, C. ga) Н, (ау, С, 93) 49, 443 
91, Yar 9з яз = С =constant |d$5— Н | (91, 92) С) Н» (91, 9», С) dq; йау 


Cylind- 
rical 
gı =P | p= С —constant 
q= 9| p=C = constant 
93 =2 | 2= С = constant 
Spherical 


qı =r | r=C = constant 
9: = 9 | = С = constant 
gs— 9| Q—C = constant 


dS —C dq dz 
dS — dp dz 
dS =p dp dp 


dS = С? sin Ө 40 dp 
dS =r sin C dr dp 
dS = г dr 40 
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